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Uber eine Verallgemeinerung der Purcellschen Methode 
zur Auflésung linearer Gleichungssysteme 


Von E. Egervary, Budapest 


Zusammenfassung. Es wird gezeigt, wie eine vom Verfasser eingefiihrte ,,rangvermindernde 
Operation“: bei Matrizen unter anderem zur Auflésung homogener linearer Gleichungssysteme 
verwendet werden kann. Das dargelegte finite Iterationsverfahren enthalt als Spezialfall die 
bekannte vektorielle Methode von KE. W. Purcell. 


E. W. Purcell hat kiirzlich ein Iterationsverfahren zur Auflésung eines Systems 
von linearen Gleichungen mitgeteilt!, das nach endlich vielen Schritten’ zum Ziel 
fihrt und die folgende geometrische Deutung gestattet?: Sei 


n 

aga = 0,0 Hei we 

j=l 
das aufzulodsende homogenisierte Gleichungssystem. Betrachtet man n zu einer 
Spaltenmatrix y zusammengefaBte GroBen x,, x,..., %, als homogene Koordinaten 
eines Punktes in einem (n — 1)-dimensionalen Raum und wird die Transposition 
durch einen Stern angezeigt, so kann die i-te Systemgleichung durch a; r = 0 wieder- 
gegeben und als Bedingung dafiir angesehen werden, daf der Punkt ry in der Hyper- 
ebene a; liegt. Man nimmt nun den Ausgang von den n Punkten 


] 0 0 

0 1 0 

t1° = ed |: Le? = Salle ? ee a é 

0 0 il 
Aus 1,° als Zentrum werden die iibrigen Punkte zunachst auf a, projiziert, womit 
man n — 1 Punkte y,}, 1,',..., 2,1 erhalt. Von diesen werden aus y,! die tibrigen 
auf a, projiziert, womit man n — 2 Punkte 7,2, 1,?,..., Yn” in a, und a, erhalt. Auf 


diese Weise fortfahrend, gelangt man schlieBlich zu einem Punkt r%—* in ay, ay,..., 
Qn—, der eine Lisung des Gleichungssystems darstellt. 


Jede Projektion verlangt gemaf 


* k—1 S 
yd a get wet, Oth go} ibe 2, ee ey n k, (1) 
kk +h >: i i Ea Ny a eat 


lediglich die Berechnung von Skalarprodukten. 


1 BE. W. Purcell: The vector method of solving simultaneous linear equations. J. Math. 
Physics 32, 180—183 (1953). ; 

2 Wir folgen der Interpretation, die G. E. Forsythe in Math. Reviews 15, 471 (1954) ge- 
geben hat. : 
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Allerdings la8t sich gegen diese Methode ebenso wie gegen die Triangulairfaktori- 
sation von Banachiewicz ein Einwand erheben: Beide Methoden versagen, wenn 
irgendeiner der linken oberen Hauptminoren der Koeffizientenmatrix 2{ verschwindet. 


Es ist der Zweck der vorliegenden kurzen Note, ein allgemeines Verfahren dar- 
zulegen, welches in der Wiederholung einer rangvermindernden Operation besteht?, 
ohne Ausnahme fiir jedes homogene System anwendbar ist und die Purcellsche 
Methode als Sonderfall einschlieBt. 


Setzen wir vorliufig voraus, daB die gegebenen Gleichungen 


I 
S 
S&S 

I 
S 
> 

I 
= 
ee 


at eee ery, (2) 


j=l 
linear unabhangig sind. Der erste Schritt unseres Iterationsprozesses besteht in 
der Berechnung von ” — 1 linear unabhangigen Vektoren, welche die erste Gleichung 
a; ¢ = 0 befriedigen. Zu diesem Zwecke wenden wir auf die Kinheitsmatrix n-ter 
Ordnung & = (¢, ¢:,...-,n) die folgende rangvermindernde Operation aus: 

e, ay 


¥,=€-——-—; 


* 
ay er, 


hierbei ist aje,, irgendein nichtverschwindendes Element von aj. X, hat den 
Rang »— 1 und geniigt der Beziehung a; X, = 0, wie man sofort sieht. Ferner 
gilt X,e,, = 0, demnach verschwindet die »,-te Spalte von X, und die verbleibenden 
Spalten sind linear unabhingige Losungsvektoren der ersten Gl. (2). 


Der allgemeine Iterationsschritt lautet: 


* 
Ee—1 ey 0, Xp—4 


Ba eas Sale a 


(3) 


ake 

) ibaa ea! a 
Die lineare Unabhangigkeit der Vektoren a; (6 1, 25.05, 7)> mlehbt age eee 
nach sich und aj, %;,-1 ey bezeichnet irgendein nichtverschwindendes Element 
von Os X;,-1 (das stets existiert). X; hat den Rang n — k und befriedigt die ersten 
k Gleichungen a; &, = 0, tes a, X, = 0. Ferner gelten die Gleichungen X, ¢,, = 0, 
Peo es = 0, so daB die »,-te, v,-te,...,-te Spalte von X, verschwinden, 


wahrend die restlichen Spalten linear unabhingige Lésungsvektoren d 
k Gln. (2) darstellen. é 8 er ersten 


In dieser Weise fortfahrend, gelangt man mit dem r-ten Schritt zu einer Matrix 
(n — r)-ten Ranges 


* 
eee yO, ky 


Sep eR 


> 


* 
a et ey, 


deren v,-te, Ye-BOy ic, " y,-te ‘Spalte verschwinden. Die n — r nichtverschwindenden 
Spalten sind notwendigerweise linear unabhangig und stellen mithin ein vollstandiges 
System von Lésungsvektoren des gegebenen Gleichungssystems (2) dar. 


3 Vgl. des Verfassers Abhandlung ,,On rank-diminishing o i i i 
é | mr : eratio d i 
to the solution of linear equations’. Z. angew. Math. Phys. ee Deasaees pai ga 


. 
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Sind nun aber die gegebenen Gleichungen nicht linear unabhangig, dann eliminiert 
unser Verfahren automatisch jene Gleichungen, welche Folgen der vorangehenden 
sind. In der Tat, ist ein Vektor a, X= 0 fir p >q, so ist ajr=0 eine Folge 
_ der ersten q Gleichungen und kann daher weggelassen werden. In dieser automatischen 


Elimination linear abhingiger Gleichungen liegt ein besonderer Vorzug unserer 
Methode. 


Um den Zusammenhang zwischen Purcells Algorithmus und unserem Verfahren 


zu zeigen, nehmen wir an, daB die ersten k — 1 Spalten von X,_, verschwinden, 
das hei®t, daB X,_, die Gestalt 


Tea Osa OAT ia rd Be) 
hat. Es gilt somit 


ee ill * = * = —_ 
Xe 7 te=% ? Oe Repay = (020%, 0h ape eer an y. as Oy Le 1 Ce = 0; LE - 


und vorausgesetzt, da a; x'—1=- 0, liefert unsere Rekursionsformel (3) 


ano ote? 

any a cae nN 

ee Ory pea aig be eke (0... HE, HS 
%, Uy Oy Uy 


Xx, enthalt infolgedessen die nachstehenden nichtverschwindenden Spalten 


* k—1 
ae 1 Khetk th 
een eae le eae h=1,2,...,.0—k, 
kee 
was mit der Purcellschen Rekursionsformel (1) tibereinstimmt. Im iibrigen la8Bt sich 
auch unsere allgemeinere Methode geometrisch als wiederholter Projektionsvorgang 
interpretieren®. 


Das Verfahren 148t sich tiberdies so ausgestalten, daB es die gleichzeitige Auflosung 
von zwei Gleichungssystemen leistet. Schreibt man diese Gleichungssysteme in 
der Gestalt 

Gf Ot te Oe ar 


y* b; = 0, Po WQARH 15 


so hat man zunachst auf die Einheitsmatrix n-ter Ordnung € die rangvermindernde 
Operation 


beds 
Be ee 
31 © b, 
anzuwenden. Hierauf folgt 
‘ Bi be a5 Bi 
aaa) ee usf., 
A> B1 be 


bis man schlieBlich nach r Schritten zu einer Matrix 8, vom Range n — r gelangt, 
deren Spalten die Lésungsvektoren y und deren Zeilen die Lésungsvektoren )* an- 
geben. Natiirlich ist darauf zu achten, da die in den Nennern auftretenden Skalare 
nicht verschwinden’. Dies wurde im oben beschriebenen Verfahren dadurch erreicht, 
daB fiir die Vektoren 6, jeweils geeignete Hinheitsvektoren e,, gewahlt wurden. 


(Hingegangen am 24. September 1957) 
. 18* 
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Die diinne Kreisringplatte mit grofer Ausbiegung 
Von K. Federhofer, Graz 


Zusammenfassung. Durch unmittelbare numerische Integration der beiden Grundgleichungen _ 
der diinnen Ringplatte mit kleiner Innenlochzahl g,; und kleiner Belastungszahl #* wird bewiesen, 
da& die in einer friitheren Arbeit erfolgte genaherte Integration dieser Grundgleichungen mit 
Hilfe des Verfahrens von Galerkin schon bei Verwendung nur eines Freiwertes dann zu praktisch 
vollkommen ausreichenden Ergebnissen fiihrt, wenn die nichtlinearen Glieder in den Grund- 
gleichungen durch jene Werte ersetzt werden, die sich aus dem Grenzfall der biegungssteifen 
Ringplatte ohne Membranwirkung (lineare Theorie) ergeben (Fall II a). Verwendet man hingegen 
bei der Berechnung dieser Glieder die Werte des Grenzfalles der reinen Ringmembran (Fall IT b), 
so ergibt sich die groBte Einsenkung um 14% zu niedrig, die gré8te Beanspruchung um 40% 
zu ungiinstig. SchlieBlich zeigt die durchgefiihrte Untersuchung, da8 die bisher tubliche Ver- 
nachlassigung des Hinflusses der Biegungssteifigkeit, das hei8t die Anwendung der reinen Ring- 
membrantheorie bei kleinen Kennzahlen zu véllig unbrauchbaren Ergebnissen fihrt. 


Einleitung 


Wenn eine diinne, am AufSenrande festgehaltene Kreisringplatte unter dem 
Kinflusse 4uBerer Lasten Durchbiegungen erfahrt, die mit der Plattendicke vergleich- 
bar oder noch gréBer als diese werden, so ist es iiblich, solche diinne Ringplatten 
mit groBer Ausbiegung nach der Ringmembrantheorie zu berechnen, das heiBt 
auf die Biegungssteifigkeit der Ringplatte keine Riicksicht zu nehmen. Fir zwei 
wichtige Lastannahmen, nimlich fiir den gleichfoérmigen Oberflachendruck und fiir 
eine nur am nicht gehaltenen Innenrand wirkende Belastung, wurde nach dieser 
Theorie die vollstandige Berechnung von E. Schwerin! durchgeftihrt, und es konnte 
fiir den zweitgenannten Belastungsfall die Losung der hierher gehorigen Grundgleichung 
— einer nichtlinearen Differentialgleichung 2. Ordnung — sogar streng in geschlossener 
Form dargestellt werden. Um die Zulassigkeit der Anwendung dieser Ringmembran- 
theorie bei der Berechnung von diinnen Ringplatten mit groBer Ausbiegung beurteilen 
und ihren Giiltigkeitsbereich abgrenzen zu konnen, ist aber die Kenntnis des Verhaltens 
solcher Ringplatten bei Beriicksichtigung ihrer Biegungssteifigkeit notwendig, dies 
um so mehr, als aus friiheren Untersuchungen von St. Way? und des Verfassers? 
bekannt ist, daB sich bei der diinnen, vollen Kreisplatte mit groBer Durchbiegung 
recht betrachtliche Unterschiede in der Formanderung und in den Spannungen er- 
geben, wenn anstatt der genauen Theorie die einfachere, von H. Hencky‘ entwickelte 
Membrantheorie benutzt wird. Fiir den nach der Membrantheorie bereits streng ge- 
losten Fall der nur am Innenrande belasteten Kreisringplatte habe ich daher in einer 
friiheren Arbeit® die beiden Grundgleichungen einschlieBlich des von der Biegungs- 
steifigkeit herriithrenden Beitrages aufgestellt und diese zwei nichtlinearen Differential- 
gleichungen 2. Ordnung sodann genahert mit Hilfe des Verfahrens von Galerkin 
integriert. Die numerische Auswertung der allgemeinen Ergebnisse bestitigte vollauf 
die urspriingliche Vermutung, daB sich so wie bei der vollen Kreisplatte auch bei der 
diinnen Ringplatte betrachtliche Unterschiede in den Spannungen ergeben, wenn bei 
deren Berechnung die Biegungssteifigkeit vernachlassigt, das heiBt anstatt der genauen 
Theorie die Ringmembrantheorie benutzt wird. Die Abweichungen wachsen mit 
abnehmender Lochverhaltniszahl (Verhaltnis des Innen- zum AufSenhalbmesser des 
Ringes) und mit abnehmender Belastungszahl; fiir Lochverhialtnisse gréBer als 0-7 
ist wenigstens die Formanderung durch die Ringmembrantheorie zutreffend dargestellt. 


1 KE. Schwerin: Z. techn. Physik 10, 651 (1929). 

2 St. Way: Trans. Amer. Soc. mechan. Engr. 56, 627 (1934). 

3 K. Federhofer: Luftfahrt-Forsch. 21, 1 (1944). 

4H. Hencky: Z. Math. Physik 68, 311 (1915). 

° K. Federhofer: Z. angew. Math. Mechan. 25, 5 (1945); Osterr. Ing.-Arch. 1, 21 (1946). 
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Da der Nachweis fiir das Versagen der Ringmembrantheorie bei kleinen Lochverhiltnis- 
zahlen und kleinen Belastungszahlen auf Grund einer nur genaherten Integration 
der beiden vollstéindigen Grundgleichungen erbracht worden ist (und zwar bei 
Beschrankung auf einen einzigen Freiwert in dem bei Verwendung des Galerkinschen 
Verfahrens benutzten Naherungsansatz), so war es noch wiinschenswert, iiber den 
Grad der Genauigkeit dieser geniherten Integration Aufschlu8 zu bekommen. 
Diesem Zweck dient die vorliegende Untersuchung, in der fiir Sonderwerte, die 
im Bereiche kleiner Kennzahlen liegen — und zwar fiir eine Lochverhiltniszahl 0-1 
und fiir eine Belastungszahl 4 —, die vollstindigen Grundgleichungen unmittelbar 
numerisch integriert worden sind. Der Vergleich mit den nach dem Verfahren von 
Galerkin und bei Verwendung nur eines Freiwertes erhaltenen Ergebnissen zeigte, 
daB diese praktisch bereits vollkommen geniigen. 


I. Die Grundgleichungen 


Im folgenden werden die schon in der friiheren Mitteilung gebrauchten Be- 
zeichnungen und die mafgebenden Grundgleichungen angegeben. 


Tis a Innen- und AuSenhalbmesser der Ringplatte, 
h deren Dicke, 


P/2r;a Belastung je Lingeneinheit des Umfanges am Innenrande, 
YD artes Biegungssteifigkeit, 
g(r) Neigung der Meridiantangente am Orte r, 
w(r) Einsenkung am Orte 7, 
EH Elastizitatszahl des Plattenstoffes, 


u Poissonsche Zahl. 


Mit uw = sin ¢ gelten fiir die Biegungsmomente G,, G, im Meridianschnitte und im 
dazu senkrechten zweiten Hauptschnitte die Formeln 
du U uU du \ 
G=—Dietus) &=-D(F +4): 


r 


Die Membrankrafte 7',, 7’, in den genannten Schnitten lassen sich mit Einfiithrung 
einer dimensionslosen Funktion y(o) und mit @ = —. darstellen durch 
a 


DY ae ace 
Gag | tt ae de 


fare 
Mit Einfiihrung einer zweiten dimensionslosen Funktion ®(¢) gemaé8 


P(0) = 2/3 u(o) 


und mit dem homogenen Differentialoperator 
a ad 1 
EAs do? 5 0 do. o@ 
lauten die vollstindigen Grundgleichungen fiir die nur am Innenrande belastete 
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worin die dimensionslose Belastungskennzahl (Plattenparameter) on in ihrer Ab- 
hingigkeit von der Belastung P und von den Abmessungen der Ringplatte durch 


le = T, \3 
oe = 4 12 780 — 2") (48) (3) 
bestimmt ist, wihrend m eine numerische Konstante ist, die gemaB 
1— 
eee ) (4) 


nur vom Plattenstoffe (uw) abhangig ist. 


Il. Lésung der Grundgleichungen 1 und 2 durch numerische Integration 


Zunichst fiihren wir die beiden Grundgleichungen mit Hilfe der Substitutionen 


Say She wes ne peas 
nek. eae pa 
in die von den ersten Ableitungen freien Gleichungen 
aC 
4 & dé Spe (5) 
2 dn * 
en og Oe” (6) 


tiber, die mit den Randbedingungen 


fiir £ = 0,2 (Innenrand) und = 1 (AuBenrand) zu integrieren sind; die beiden ersten 
Bedingungen bringen das Verschwinden des Spannungsmomentes G, am gelenkig 
ausgebildeten Innen- und AuBenrand zum Ausdruck, die dritte und vierte Bedingung 
die Festhaltung der Rander (ausgedriickt durch Nullsetzen der tangentialen 
Dehnung «¢,). 


Die numerische Integration wird durchgefiihrt mit der Annahme uw = 1/; (m = 4/5) 
fiir den Sonderfall eo, = 0:1 und #* = 4. Dann lauten die beiden Grundgleichungen 


d2 
96 Fe +1 =0, 


dn ss 

$e 6 +4§=90, 

und die Randbedingungen 

dn 
dé 

ao oe 

fiir €, = 0-0Y und 4,9; 

Nach numerischer Ermittlung von 9 Naherungsfolgen, wobei mit dem aus der 


friiheren Naherungsrechnung entnommenen Ausgangswertepaar 7; = 0°133, ¢; = 


= 0°0235 begonnen wurde, ergab sich fiir den Innenrand bis auf die vierte Dezimale 
genau 


nH, = 0°1340977, ¢C; = 0°01973649 
und fiir den Au8enrand bis auf die dritte Dezimale genau 


Nq — 19348502, €, = 0°6453855. 
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9 : d 
Die ersten Ableitungen a und = am Innenrande wurden unmittelbar aus den 


Randbedingungen fiir den Innenrand entnommen, jene fiir den AuBenrand wurden 
beim schrittweisen Vordringen vom Innen- zum AuBenrande zu 


"i = 0°4489819, es = 0°4301891 
berechnet. 


Hiemit liefern obige Randbedingungen am AuSenrande 


3En, — Yq = — 0:0015559 (anstatt Null), 
3&1 — 26, = — 0:0002037 (anstatt Null). 
Die Einsenkung w(e) an der Stelle 9 berechnet sich aus 
h 
— d _—— ee \ 


bei Beachtung der am AuSenrand 0 = | zu erfiillenden Randbedingung w(1) = 0 zu 


h 
wleh= sig | Clg) da. 

Die Tabelle 1 enthalt die hienach berechneten Ordinaten der Meridianlinie fiir 
eine Reihe der o-Werte zwischen 071 und 1; in ihr sind auch die nach dem Galerkin- 
schen Verfahren berechneten Ordinaten zum Vergleiche eingetragen. Dabei sind zwei 
Falle (bezeichnet mit Fall IL a und IIb) zu unterscheiden, entsprechend den beiden 
Grenzfallen der biegungsfesten Ringplatte ohne Membranwirkung und der reinen 
Ringmembran. Fiir beide Grenzfalle sind die zugehdrigen Grundgleichungen ge- 
schlossen integrierbar, so dai: Naherungslésungen der in strenger Form nicht integrier- 
baren Gl. (1) und (2) dadurch zu gewinnen sind, daB die darin vorkommenden nicht- 
linearen Glieder durch jene Werte ersetzt werden, die sich aus einem der beiden 
Grenzfalle berechnen lassen’. 


SchlieBlich enthalt die Tabelle 1 auch jene Ordinaten, die sich unmittelbar aus 
der Ringmembrantheorie ohne Berticksichtigung der Biegungssteifigkeit ergeben. 


Aus den Angaben der Tabelle 1 ist zu entnehmen, daf die unter Zugrundelegung 
des Falles IIa nach dem Galerkinschen Verfahren mit Benutzung nur eines Frei- 
wertes berechnete Forminderung der Ringplatte praktisch bereits vollkommen mit 
der durch unmittelbare numerische Integration ermittelten Formanderung iberein- 
stimmt; nicht so giinstig liegen die Verhaltnisse im Falle IIb, in welchem sich bei 


Tabelle 1. Hinsenkungsverhdltnisse w/h 


umeneché Galerkinsches Verfahren Heind Wiisbetea 

; pase Fall IIa Fall IIb Hall Auk 
Ol 0°35103 0°35005 0°30079 | 0°68100 
0°22361 0°30658 0°30721 0°24214 0°55 
0°34641 0°26168 0°26328 0°19428 0°445 
0°4 0°24124 0°24310 0°17525 0°405 
0°56569 0°17554 0°17759 O2TTS 0°278 
0°69282 0°12371 0°12541 0°08321 07193 
0°77460 ~ 0:09026 0°09158 0:06003 0°14 
0°91652 0°03291 0°03345 0°02165 0°053 
0°95917 0°01600 0°01625 0°01051 0°025 
1:0 0 0 0 | 0 
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Anwendung des gleichen Naherungsverfahrens Werte fir die Formanderung ergeben, 
die gegeniiber den genauen Werten um rund 14% zu niedrig liegen. 


Das Versagen der reinen Ringmembrantheorie im Bereiche kleiner Lochverhaltnis- 
zahlen und kleiner Belastungskennzahlen ist aber erwiesen durch die Zahlenangaben 
der Spalte 5 in Tabelle 1, die um rund 94% hoher liegen als die nach der strengen 
Theorie berechneten Einsenkungswerte. 


Die in den Randfasern + ue durch die Spannungsmomente G4, G, erzeugten 


Biegungsspannungen betragen 


h2 , h2 
EH 7 V1 bzw. EH Te Vo» 


worin die Zahlenwerte y, und y, bestimmt sind durch 


1 d® @ 
V3 dat) beam tuys 

1 © d® 
ema er aar) ic. +a) 


Analog lassen sich die von den Membrankraften 7, 7’, herriithrenden Spannungen 


2 2 
darstellen in der Form # se ei Heese mit 
‘a la 


1 on bike 1 dV 
Ise eS tes Too py ao. 


Uy = 


Tabelle 2. Die Beiwerte y,, y, der Biegungsspannungen sowie 7), t, der Membran- 
spannungen ftir den Innen- und AuBenrand der Ringplatte 


at 5 
OF = 45 1 = oa V1 Ye | Gl Ts; 
Numer. Integration. . 0 1°9355 0°18503 0°061677 
imnenrand Balla cope rent 0 | 1863 07182 0:061 
6; OL SRalbsTT vis hs sete cites 0 1 227535 Ora 0°057 
Ringmembrane ..... 0 0 0°870 0°261 
Numer. Integration. . | — 0°00016 (anstatt 0) 0°19458 0060505 0°020156 | 
AuBenrand Pall Il ai... 20... 0 0°198 0:061 0-020 
| Fall IEE: ong Rahim es 0 0°128 0°037 0°012 
Ringmembrane ..... 0 0 0-180 0°054 


Die Tabelle 2 enthalt eine Zusammenstellung der fiir den Innen- und AuBenrand 
berechneten Spannungswerte y und t. Auch hier stehen die fiir den Fall II a genahert 
berechneten Spannungsbeiwerte fiir die Biegungs- und Membranspannungen in sehr 
guter Ubereinstimmung mit den durch numerische Integration der Grundgleichungen 
erhaltenen Werten. Die fiir den Fall II b durchgefiihrte Niherungsrechnung beurteilt 
die maximalen Spannungen (tangentiale Spannungen am Innenrande) um rund 40% 
zu ungiinstig, wahrend die nach der reinen Ringmembrantheorie gewonnenen Ergebnisse 
zeigen, dali die gebrauchliche Vernachlassigung der Biegungssteifigkeit hinsichtlich 
der Beanspruchungen ein vollkommen unrichtiges Bild gibt, denn die am Innenrand 
auftretenden tangentialen Biegungsspannungen betragen mehr als das Doppelte der 
maximalen reinen Membranspannungen. 


(Hingegangen am 25. September 1957) 


H. Jung: Bestimmung der Warmespannungen in ungleichférmig erwairmten Kontaktéfen 257 


Uber die Bestimmung der Warmespannungen in ungleichférmig 
erwdrmten Kontaktéfen* 


Von H. Jung, Stuttgart 
Mit 6 Textabbildungen 


Zusammenfassung. Es werden das Temperaturfeld und die dadurch auftretenden Warme- 
spannungen untersucht, wenn an einem Kontaktofen der Ammoniaksynthese eine Betriebsstérung 
auftritt. Wir zeigen, daB ein nicht rotationssymmetrisches Temperaturfeld die Betriebsspannungen 
herabsetzt, so da8 durch eine Betriebsstérung der Ofen nicht gefaihrdet wird. 


I. Einleitung und Problemstellung 


Betriebsstérungen in Kontaktéfen zeigen sich durch eine ungleichmaBige Er- 
warmung des Mantels. An der Oberflache des Mantels lassen sich durch Messungen 
mit Thermocolorfarben und Thermoelementen die Temperaturfelder einfach ermitteln. 


An den Betriebsingenieur tritt nach Bestimmung der Temperaturen die Frage 
heran, ob die von ihm ermittelten Temperaturen im Kontaktofen nicht auf unzulassige 
zusatzliche Warmespannungen fiihren. 


Erste Untersuchungen tiber nicht rotationssymmetrische Warmespannungen 
wurden von F. Schau! und A. Stodola durchgefiihrt. Hierbei wurde angenommen, 
daB die Temperaturverteilung sich mit der Rohrlinge nicht andert. Bei dem vor- 
liegenden Problem andert sich die Temperatur sehr stark mit der Rohrlinge, so daB 
wir diese vereinfachende Voraussetzung im vorliegenden Fall nicht anwenden kénnen. 


Wir wollen dieses fiir den Hochdruckapparatebau wichtige Problem erneut an- 
greifen. In Kap. II geben wir ein allgemeines Verfahren zur Bestimmung der 
Temperaturverteilung in Rohren an, wenn die Temperaturverteilung des AuSen- 
mantels und die Warmetibergangszahl am AuBenmantel bekannt ist. Das in Kap. II 
angegebene Verfahren vereinfachen wir in Kap. III. Die in Kap. III erhaltenen 
Temperaturverteilungen vereinfachen wir in Kap. IV noch weiter so, daB wir das 
vorliegende Problem in eine Aufgabe fiir den ebenen Spannungszustand und eine 
Aufgabe der Plattenbiegung aufspalten konnen. Diese beiden Teilaufgaben behandeln 
wir in Kap. V. Ein Zahlenbeispiel soll dann noch in Kap. VI durchgerechnet werden. 


II. Bestimmung der Temperaturverteilung 


Der Kontaktofen (kreiszylindrisches Rohr mit AuSenhalbmesser r, und Innen- 
halbmesser 7;) bestehe aus einem isotropen homogenen Werkstoff mit einer von der 
Temperatur unabhangigen Warmeleitzah] A. Das den Ofen umgebende Medium habe 
die konstante Temperatur 7. Die Warmeiibergangszahl an der Aufenwand des 
Ofens sei «,. 

Die stationare Temperaturverteilung 7’ im Ofen geniigt dann bei Verwendung 
eines Polarkoordinatensystems (7, gy, z) der Differentialgleichung 


gh I gh: ES eXeh eg e) eT 


or? +e or a Cl | a ee 0: (1) 
Die Randbedingungen lauten 
deel AG's), 
a ea fir r= 14. (2) 
Ln ray aes 


* Vorliegende Untersuchung wurde mit Unterstiitzung der BASF Ludwigshafen durchgefihrt. 
1 BE. Melan und H. Parkus: Warmespannungen, S. 32. Wien: Springer-Verlag. 1953. 
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Die Lésung der Differentialgleichung (1) setzen wir in der Form an: 
T=T, +a +a,mr+ >’ T, cosm@. (3) 
(3) in (1) eingesetzt ergibt 


Er eh ac De 
or r or | @2 


ney a == (4) 


Mit (3) gehen die Randbedingungen (2) tiber in 
T,=T7, +a, taq,lnr, +2) T alta 2) 008 me; 


T (5) 
0 ee, PSE 4 Ty) _ 20a th @. 
Xe i or he a 


Xe = 


Ist Q die bei einem bestimmten Durchsatz des Kontaktofens frei werdende Warme- 
menge, so konnen wir a) und a, so bestimmen, da diese Warmemenge gerade durch 
ein rotationssymmetrisches Feld abgefiihrt wird. Wir erhalten dann fiir die Temperatur- 
verteilung, wenn h die Lange des Ofens ist, 


on : 
(Macs Wide fal eee - 4 In 2) + >’ T (7, 2) COS Mm @. (6) 
(6) in (5) eingesetzt, fihrt auf 
Q ise) 
Dl hee aera a Sai Mira ai?) C08 
a a (7) 
i en 
0 ce i or T’,, Rati 
Mit der Reihenentwicklung der Oberflachentemperatur 
ao 3 —- Q 
Dt Tale eonatiie: To Les a Rt (8) 
erhalten wir 
P nl as 2) = Tnl(2). | (9) 


Durch die zweite Gleichung aus (7) und durch (9) sind uns die Randbedingungen fiir 
die Funktion 7',,(r, z) gegeben. 


Die Lésung 1a4Bt sich mit Hilfe der Fourier-Transformation durchfiihren. Mit 
den Zylinderfunktionen J,,(ip7) und N,,(ipr) und den Abkiirzungen 


: ‘ Lp ih 
Z,(i pt.) =I mi pte) +5 


[Jim _1(t_p Pa) aa Spe -1(t_p a ih 


Zip 1) = Nip ra) +2 [Nn sli Pe) — Nn (69 70)] 


erhalten wir nach einiger Zwischenrechnung 


T=) 4+557(++4+m), 


Q2ahi See ies r 
ive | Zeer) T ep Repea ee ee oe 
dale s al? P Tq) Y mit Pr) — CPT, mt p Pr) aah ip(e— 
+ Ba = come \ Tipe, Tali Dre) = Dpto eD | T(t) rc eae 


@ 


Die Auswertung der in (10) auftretenden Integrale wird nur in den wenigsten Fallen 
moglich sein. Wir sind daher gezwungen, um praktisch brauchbare Ergebnisse zu 
erhalten, eine speziellere Lésung, die den im Betrieb auftretenden Temperatur- 
verteilungen méglichst gut angepaBt werden kann, zu ermitteln. 
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Il. Ableitung einer speziellen Lisung 


Die meisten im Betrieb auftretenden Temperaturverteilungen lassen sich durch 
Deg) = 0 26 57 A, err ge" 003 32 — a) (11) 
annéhern. Fir viele Falle kénnen wir uns auf die Temperaturverteilung 


LAS OC) = 2" A eo? cos (11a) 
beschranken. 


Gehen wir mit der Temperaturverteilung (lla) in (10) ein, dann wird 


wth _ +o | 
T = 1+ 4 7940 (pe — 7) inte + coe FV FV Fetter. —n + | a 


2 
x poi a2 


(posi) eine.) be Rap 


Das in (12) auftretende Integral 14Bt sich mit Hilfe der GauB-Transformation aus- 
werten. Wir erhalten 


Ta %a 


gee te ee = Fe) In “2 + A cos $y/ Te erm cos 2u2z(r7g—1T) + 


(13) 


ay ed) 
+ ( at : ee e#* cos 2u2C ae 


2a 
0 


Beachten wir, daB u(r, — r)? <1 ist, so k6nnen wir die Exponentialfunktionen in (13) 
durch ihre Potenzreihen ersetzen. Das in (13) noch auftretende bestimmte Integral 
148t sich dann auch auswerten. Es wird damit die Temperaturverteilung 


im 1 
fe he 4 788 (T* — T',)In =2- +A cost / *e fern cos 2u2(%a—1) + | 
1 wd 1 : 
-( 2 au 7 l(a (ra — 1) — gr) sin 2 uz (re — 7) ao (F4) 
= “" cos 212 (ra —1)|}. 


Fir z = 0 und r = 7, erhalten wir das fiir den Betriebsleiter wichtige Temperatur- 
maximum. 


ca Tepe eactigh 
Te Ee Tin + A \/ 2\e ated Gee 2, ) (15) 


Durch (14) ist ein unseres Erachtens fiir den Betrieb brauchbarer Ausdruck gegeben. 


IV. Umformung und Vereinfachung der Temperaturverteilung 


Zur Bestimmung der Warmespannungen miissen wir die Temperaturverteilung 
noch vereinfachen. Wir fiihren zu diesem Zweck das durch Abb. 1 gegebene kartesische 


Koordinatensystem ein. Es ist, damit 


Li TOs Yi ils 2%. 
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Mit diesen Verianderlichen wird der nicht rotations- 
symmetrische Anteil der Temperaturverteilung (13) 


fice # fy) dee 
T = A cos VE Ween cos 2uzy + 
to 1 ? oe 
+( ai 7,3) oF | et cos 22 y dy. 
1 
Ist 1 >> und 1 y? <1, so wird, wenn wir Vrs 


und die Exponentialfunktionen in Potenzreihen ent- 
wickeln; wobei von diesen Potenzreihen nur die linearen 
Glieder beriicksichtigt werden, 


Abb. 1. Koordinatensystem 


T= A cos 5—(1 +> ce Cae 1 + (4472 —~+)4|. (17) 


Ta 


Es ist zweckmaBig, die dimensionslosen GroBen 
g a & 


Oe = 
ee te Ir, 2r 


= 2s Ds ae 


a 


n (17) einzufiihren. Beriicksichtigen wir nur die linearen Glieder in y, dann wird 


es _ 


T = A Gos per? 1p (18) 


Die Temperaturverteilung (18) spalten wir in einen symmetrischen und asymmetrischen 
Bestandteil auf. 


Mit 
= 6 6 1 
Rae Nae ay Diss Aap (ta 2t3) 
erhalten wir 
T=7, 473 

~ ee 2 Xt, 6 
T,, = A.cos GE-#* |\T=- “at I> (19) 
T, = A*%a%e £ G9 % eK, (20) 


Die Vereinfachungen, die auf (19) und (20) fiihren, bilden die Rohroberfliche auf 
eine ebene Scheibe und eine Platte ab. 


Bezeichnen wir mit 7 die Querdehnungszahl und » die lineare Warmeausdehnungs- 
zahl, so kénnen wir die durch die Temperaturverteilung (19) bedingten Spannungen 
aus dem Verschiebungspotential bestimmen. Das Verschiebungspotential geniigt der 
Differentialgleichung 


ary ay ee 2 6 
mae tos =f es =) yA cos Be-##(] + “Sate 2] (21) 
Bezeichnen wir mit G den Schubmodul, dann sind die Spannungen 
a Oy 2 
Ong = — 24 Cn = 2G oe = 1¢— (22) 
Ox Oz 


Die Spannungen (22) miissen den Grenzbedingungen geniigen 


z—> + ©, Cxg = Cyy = Oz > 0. (23) 
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Die Temperaturverteilung (20) fiihrt auf Spannungen, die sich mit Hilfe der Platten- 
durchbiegung w berechnen lassen. Die Durchbiegung w muB der Differentialgleichung 


gentgen: 2 
Ow ew +o(1+—)A 


Ox? 622 


Aus der Durchbiegung ergeben sich die Spannungen 


2 Las = 
a Cos Pere? 10), (24) 


Hé 1 ew C?w 1 2 r Se ee ’ 
Rieke SHES} Prot ga t(he  )A og cee mee, 
q 
HE 1 ew ew 1 Oa, Sateen 
eo (eh oat oat ts A rt ase ae “*, ¢ (25) 
| 
HE C2w 
Lt) eras 


J 


Die durch (25) gegebenen Spannungen miissen ebenfalls den Grenzbedingungen (23) 
genugen. 
V. Bestimmung der Warmespannungen 


Zur Lésung der Differentialgleichung (21) machen wir den Ansatz 


VW = cos D(z). (26) 
(26) in (21) eingefiihrt, ergibt 
a2@ = ~ 1 2 a‘’‘a 0) i 2 
~~ @) =(14+—) vA (1+ 25) em : (27) 
Die den Grenzbedingungen (23) gentigende Lésung der Differentialgleichung (27) lautet 
am) 3 
a Cae 1 a \ -(S +ivz 
O@) = —A(1 +8726) (1 +4)Vz\e (a lap: (28) 
Damit ist die Verschiebungsfunktion fi 


P(E, 2) = — Acosz(1 + “478 6)(1 + ; )» Velie lap (29) 


Nach einiger Umformung wird 1 ) 
1 9 2Vae 
Se os a a’a Ue a ee eB = 
(@, 2) = —Acosz(1 + 297 §\(1 +7) peu ve \e dl Go} z — 
> (30) 
z 2 4 
= | _ G+ 2uo? ; | _ (2nd)? 
—// Zing +6 \e ae 6 d= ere We Veet aa 
ag 0 0 H 
In entsprechender Weise erhalten wir die Durchbiegung w zu 
1 
: why Vie 
w(@, 2) = A cos # "272 (1 + | s-e4# | \e® do Gofz — 
he (31) 


2 ‘ 2 Vu 
se oY a wae 4 ee eae? 
—~//2Ginz+e ets idee te Ne 4a a 
7 0 


0 


ES 
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Aus (30) und (31) bestimmen wir die durch das nicht rotationssymmetrische 
Temperaturfeld (18) hervorgerufenen Spannungen. Fiihren wir die urspriinglichen 
Koordinaten 7, y, z ein, dann ist 


2 i) 
1 Gaiam = @ mile 
Sop =2GA (1 | deat fae ul cos £6 (gf) — be (a) |, | 
1 Wives A g 
On =—2GA(14 [14 o ite r) cos £5(5—), 
1 
pry » Be cree Nias (32) 
On =2GA (1 + | = f came r)| sin $e47 vere Dy Je dg + 
227 | 227, x 
— Z _ (+ 20d Hoy \ Q=2008 gp 
+ |/ = Gof 5 +e?" Je 4u 6 AC +e 7% Se 4h : 
lb Ta 0 0 


In (32) haben wir zur Abkirzung gesetzt 
1 


} 

rr 2 a 

a 2 = ae C2 7s Z Ht ee a "a a+2ngy 
Slat )=| flere at [Baier teal a 
0 0 \ 

+ (33) 
2/27 
tt Nal BED) 
—e 7% je ae ae 


Die in (32) und (33) auftretenden Integrale lassen sich einfach mit Hilfe der Tabellen 
von Jahnke-Emde? berechnen. 


Fir die AuBenwand des Rohres wird 


Opp = 2GA(I =) ”_ cos r/o ( aexamaal 


On = —2G-A (1 + -) 5 COs eS 


47 Vu 
1 
Tie) geese 1 2s : Zz (34) 
= | 16 & seme yee 
Oe = 2GA (1 7 )ysinte TM maT dl Sm 
0 
2/297 2/27 


Aus (32) ergibt sich, da im Rahmen dieser Anniherung die Radialspannungen durch 
die nicht rotationssymmetrischen Temperaturverteilung nicht beeinflu8 werden. 
Aus (34) entnehmen wir, das das Temperaturfeld die Spannungen an der AuBenwand 
herabsetzt, so das durch die ungleichformige Temperaturverteilung der Kontaktofen 
nicht gefaéhrdet werden kann. Wir wollen an einem Zahlenbeispiel untersuchen, wie 
die im Lautfe der Rechnung eingefiihrten Vereinfachungen sich auswirken. 


* Jahnke-Emde: Tafeln héherer Funktionen, 8. 23. Leipzig. 1952. 
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VI. Zahlenbeispiel 
Zur Durchrechnung des Zahlenbeispieles wihlen wir 
rg=O07m; r;,=05m, «, = 10 keal/m?h °C; 4 = 50 keal/mh °C. 
y= 1210-8 sa; G = 0'85- 10% kg/em?; 4 = 3. 

Lo 20° CP Tt = 110° Cr A= 90°C; pp = 1,0996. 

Mit diesen Zahlenwerten wird (14) 
T(r, p, 2) = es | 
= 110 + 126 n-* + 90 / 2 cos -$ fe-t0006 [2-(,—') cog 2°1992 z (7, — 7) — 


sta Sa 35 
— 9-234 2" [(1 + 1.0906 (r — 7)? — SP) sin 219922 (7, — 7) + oe 
Zz 
+2 a cos 2°1992 2 (r, — r) \ 
n Abb. 2 gibt die Temperaturverteilung 


als Funktion von z und g fiir r = 1;. 
In Abb. 3 haben wir die Temperatur- 
verteilung fiir g = 0 als Funktion von 
z mit dem Radius als Parameter auf- 
getragen. Da wir die Temperatur- 


ik jets verteilung (14) zur Berechnung der 
g = Warmespannungen durch (18) an- 
740 : ° : 

[_ nihern, ist in Abb. 4 die Temperatur- 


differenz (7’ — 7’) als Funktion von z 
mit r als Parameter gegeben. Aus Abb. 4 
entnehmen wir, daf fiir das vorliegende 
Zahlenbeispiel die durch (18) gegebene 
Annéherung ausreichend ist. 


025 


Ho 


Abb. 2. Temperaturverteilung als Funktion 
von z und @ 


140 


120 


100 15 2 


hy 
Abb. 3. Temperaturverteilung, y = 0 Abb. 4. Temperaturdifferenz, 7 — 7' 


Wir bestimmen nun noch die Warmespannungen. Abb. 5 gibt die Warmespannun- 
gen durch den rotationssymmetrischen Anteil der Temperaturverteilung 


“fe 
8 == 1628 In are 
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Abb. 6 gibt die aus (32) ermittelten Warmespannungen fiir gy = 0 als Funktion von z 
fiir r—r,;. Die Zahlenrechnung zeigt, da sich die Spannungen fiir r = 7, nur un- 
wesentlich andern, so daB® wir auf ihre Angabe verzichten konnen. 

Den Warmespannungen  iiberlagern 
sich noch die Betriebsspannungen durch 


2,0 
1y=05 lm =06 14°07 [in] 


Abb. 5. Rotationssymmetrische Spannungen Abb. 6. Zusatzspannungen 


den Innendruck p;. Bekanntlich ruft der Innendruck p; eine tangentiale Zugspannung 
hervor. Diese tangentiale Zugspannung wird im vorliegenden Beispiel des Kontakt- 
ofens durch das nicht rotationssymmetrische Temperaturfeld herabgesetzt. 

Dieses im ersten Augenblick tiberraschende Ergebnis ist unseres Erachtens fiir 
den Betrieb von Kontaktéfen von groBer Bedeutung, da es dem Betriebsingenieur 
erhohte Sicherheit bei Betriebsst6rungen gibt. 

(Hingegangen am 4. Oktober 1957) 


Transverse Bending of Thin Shallow Shells of Translation 
By Gunhard-Aestius Oravas, Venezuela 
With 4 Figures 


Summary. Thin translational shells are investigated in this treatise for transverse boundary 
perturbations. Mushtari-Vlasov thin shallow shell equations were employed in the derivation 
of the translational thin shell’s solution in terms of normal displacement function. This solution 
is implemental in the satisfaction of higher order boundary conditions (displacement, stress 
couple and rotation) of the thin shell of translation. The remaining boundary conditions are 
satisfied by the momentless state of stress. 

An illustrative example demonstrates that the transverse bending region is actually a narrow 
zone around the edges of the thin translational shell. It also seems evident that along most part 
of the boundary the transverse bending is very little dependent upon the coordinate parallel 
to the boundary. In the vicinity of the shell’s corners the transverse bending is heavily dependent 
upon both coordinates. 

Many boundary supporting arches may be of such proportions, which can offer considerable 
rotational restraint to the edges of translational shells. In such cases it is advisable to design the 
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thin translational shell both for simple and fully fixed support. Full restraint condition exists 
at the common boundary of the twin thin shells shown in Fig. 3 


Synopsis 
A number of papers* have been written on the momentless theory of thin 


translational shells. The transverse bending of translational shells on the other hand 
has received relatively little attention. 


In this paper an approximate solution is developed for boundary perturbation of 
shallow shells of translation which enclose a square base. Some general conclusions 
are drawn. 

Introduction 

Thin shallow shells of translation, with positive Gaussian curvature = 2 = > 0) 
and enclosing an almost square base, are investigated in this paper for Geisaucta 
bending perturbations propagated by normal surface loading. 

The basic shallow thin shell theory, that is utilized in the following development, 
was derived in terms of generalized orthogonal curvilinear coordinates by Mushtari?® 
and Vlasov”. 

An equivalent thin shallow shell theory was developed by Marguerre® through 
energetical considerations in “almost’’ Cartesian coordinates. Nazarov arrived 
at the same set of linearized equations by poignant analysis of shallow shell’s 
differential geometry. 

According to Vlasov the radii of the curved generatrices, which describe the middle 
surface of the shell, can be assumed to be constant for shallow shells (ratio of shell 
rise to span << 1/5). This approximation in practical design has negligible influence 
on the final solution. 

The two fundamental shallow shell equations admit an elegant reduction to one 
single complex differential equation which yields to an elementary solution for a 
translational shell over a square base. 

To the writers knowledge the true complex variable method was first used in 
problems of elasticity by Russian mathematician G. V. Kolosov in his doctorate 
dissertation’ at the University of Tartu (Dorpat), Estonia in 1909. 

The ‘‘complex”’ method of analysis introduced in this paper is not a true complex 
variable method as outlined by Kolosov in his dissertation. The “complex”? method 
used here was introduced into the theory of shells by Goldenveizer?* in 1939. Since 
then it has been further developed and amplified by Russian mathematicians 
Novoshilov# 1%, Goldenveizer, Lure, Georgian mathematician Vekua and 
Estonian mathematician Alumae (Alumyae). 

On the classical methods of shell analysis (neglecting the transverse shear 
deformation and excluding any possible finite displacements in the behavior of thin 
shells) the book by Goldenveizer® gives an extensive and rigorous treatment on the 
theory of shells with applications of ‘‘complex’’ methods. 

In Western Europe the German elastician W. Zerna and English mathematician 
A. E. Green developed the complex methods of analysis in connection with classical 
and non-classical theory of shells®. 

In America E. Reissner! 16 and P. M. Naghdi’ have contributed significantly 
to this aspect of thin shell analysis. The contributions by Green, Zerna and Naghdi 
do not neglect the transverse stress deformation upon the general behavior of the , 
shell and therefore, are theoretically more accurate than other available expositions. 


* See 1, 2, 6, 14, 17, 18, 19 and 21. Numbers refer to Bibliography. 
Ingenieur-Archiv XI/ 4 19 
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The practical application of the theory of shells in the field of structural engineering, 
however, makes many simplifications possible. The transverse shear deformation 
has very small influence on the outcome of the solution for such structures and is 
always neglected under normal conditions of loading. Also the shallowness of the 
shell’s geometry enables still further simplifications to be introduced in the rigorous 
classical theory of shells. 


The final outcome of this judicious reduction is expressed in Mushtari-Vlasov 
equations, that are used as a starting point of analysis in this paper. A solution of 
the said differential equations is demonstrated in the following work for a thin shell 
of translation spanning a square base. Even though the derivation of the simplified 
equations is very difficult and a task worthy of the most talented in the profession, 
the solution obtained from them in this paper is elementary and easy to apply. All 
the steps required in the execution of this solution should be easily understandable 
for engineers with the usual background in college mathematics. 


For the consolation of structural engineers, whose calculus has become rusty 
over years of inactivity in theoretical work, it can be ascertained that the knowledge 
of higher mathematics is not necessary for the sole purpose of applying the ultimate 
solution to practical problems. One such representative example is included to illustrate 
the behavior of the translational shell under the action of transverse bending. 


Finally readers’ attention is called to the effort made in this paper, following 
leading authorities in this field, to propagate a practice of more precise terminology 
by referring to internal resisting moments and forces as stress couples and stress 
resultants respectively. 


Fundamental Differential Equation for Shallow Shells 
Formulation of Complex Differential Equation 


The generalized differential equations of Mushtari-Vlasov for thin shallow 
shells are 
V4d— EnV? w= — (1 —7)V?T— EhV? (y T), 
y 1 1 (1) 


Di Raetiers 


lo Pn 
V4w + As) VO = D 

where 
I’ = external load intensity potential function (see p, and Ppp below), 


Pn = load intensity normal to the middle surface of the shell, 
or’ 


Pa ~ ~~ “Aa | Intensity of tractional loading acting in direction of « and f coordinates 


oe ar respectively, 
St eas Bos) 


w = Displacement function (w is normal to shell’s middle surface), 
® — Stress function, 

h = Constant shell thickness, 

D = EARF/12 (1 — +), 

y = Coefficient of thermolinear expansion, 


Dc] (or) a Temperature distribution function that describes the differential 
thermal increase or decrease from a stress and strain free temperature level 
of the shell bringing about an extension or contraction of shell’s middle surface, 
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E = Young’s modulus, 


v = Poisson’s ratio. 


The differential operators convey the following meanings: 
Sor le v[ye hee a{A a]\ 
f tala sal + ale al) 


— AB\é@a|A d«|" OB|B a 
Ra Tie a ee se Ofsted 
Va an \ealep aol tole sol 
and 
Y4=/2p2, 


Na 


® 


Na 


Fig. 1. Differential Element of Shell’s Middle Surface with Stress Resultants, Stress Couples, 
Surface Load Intensities and Embedded Curvilinear Coordinates 


The coefficients A? = [A(a, 8B)? and B* = [B(«, B)}? are the first order funda- 
mental magnitudes of the middle surface of the thin shell. The radii 1/R, and 1/R, 
are the curvatures of the curvilinear coordinates 6 = constant and « = constant 
respectively (Fig. 1). The parametric coordinates « and f are assumed to coincide 
with the principal lines of curvatures on the middle surface of the shell*. 


The sectional quantities, stress resultants and stress couples, are expressed 
conveniently through stress function ® and displacement function w? (Fig. 1): 


N le 4) (= e) iF de oe 

ae —~|B OB\B Op A? B @a 0%\? 

N Ld 1 o@ 1 oA o®@ 

w=|q ale a) A B? ap a 
N N 1 ( 2D 1 oB e@ 1 0A OD 
OD erase BO A B |\ dx op B oa op A Op da’ 
Q.=—Dl| az Vv), 


=—Dlza Vw), 


* For more detail on the differential geometry of the surface, see ‘‘Differential Geometry 
of Three Dimensions” by C. E. Weatherburn, Cambridge at The University Press, 1955. 
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aia te) + ame op (ae) +) le apa ae) + ae ae (ae lb 
Bala ap) + aw te (aa)| +7 [a ae (aap) + ame ae (oe) 


Mug = My =—DA—” (ap [seap — Bee (oe) ~ a oF (oe Ih 


It seems convenient to introduce a complex function of the form 


V=w+i09, 


where w designates as yet unknown constant and ¢ = //— 1. 


Multiplying the first equation of (1) by iw and adding the resulting expression 
to the second equation of (1) yields 


V4(w +in®) —V*(ioBhw— ®) = 


ite (re: : a) A (1—»)toFl—EhioV@y 7). 
On the left side of this equation the first term already has the desired dependent 
variable, but the second term’s dependent variable is not the desirable one. As w 
in that relation is yet a free constant, the desired form of dependent variable is 
imposed, hence 
2 (ioKhw— + %)=iobh?(w +a, | 
D Si LhDo i 
Now 


wtappe P=wtio®. 
In order that the two complex quantities can be equalities, both their real and 
imaginary parts have to be equalities. 


The real parts are already equalities. The imaginary parts have to yield to the 
equation 

a : 

hep oon 


Rearranging 


eet 
ay BA 


/ To PI 
ey ie Dae eee 


The complex shell equation can finally be expressed by 


or 


Th eee ee c=! )r r= 

ae 2) 
Di Ian 1 -V12 (1 — v ( 
= it — ay (ae te i eee 


For shallow shells the first order fundamental magnitudes of its middle surface 
behave almost as constants’, i. e. 


0 ) 
Fa (4 v) eA, 
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Consequently it is feasible to propose absolute coordinates 


ArBwrl 
and 
OS g; 
B=. 
The differential operators become simply 
_ @ e 
y2 = AD “On? 
and 
2s Gi: e 1 
= ela + orl me 


A tacit assumption was made that the radius 


Ro Re (i + + - .. .), 
where 
|u| <1. 


Then it follows that for all practical purposes 


1 1 1 


n, Re Re 


In addition for shallow shells the tractional surface loads p; and p, are insigni- 
ficant under normal conditions in structural practice and, therefore, can be neglected 
in the analysis. 


When all these arguments are taken into consideration equation (2) simplifies to 


V2(V2 12) V = iB (yT), (3) 
where 
rie V 12 (19) 

Rh 


The solution of (3) can be constructed from three part solutions 


Vere Vip ar es 


where 
V, represents the particular solution, V, and V, are extractable from equations 
(V2? —1 #7) V, = 90, (4) 
V2 V, — 0. (5) 


In conformity with the assumption that the thin translational shell persists 
essentially in a momentless state of stress and that its transverse bending is effectively 
restricted to a limited narrow zone lining the shell boundaries, the perturbation zone 
can be considered as a local disturbance of the primitive membrane state of stress. 
Tf there exists such a membrane state of stress, provided the transverse bending 
does not affect the equilibrium equations to a noticeable degree, then the solution 
of equation (5) may be dropped and the membrane solution used instead. 

Hence the boundary perturbation is resigned to depend upon solution of equa- 
tion (4) and the particular function. 
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Complementary Function 
Operational form of equation (4) 
(V2 —7 1?) V, = 0 
is the representative of the expression 


o oO : 
(ee On? Pat 22) V; cae. Oe (6) 
In analogy with ordinary differential equations the introduction of solution 


V,= et tin 
(a and 6 are arbitrary constants), in equation (6) gives the condition 
e+bt—172=0 
to be satisfied by a and b, hence 
a ere Ss 


In view of the fact that the stress and strain patterns for the thin shells under 
consideration have 4-ply symmetry, the analysis henceforth is restricted to one side 
only (7 = 0). From a study of boundary conditions it can be determined that 
function V, has to have periodicity with coordinate &. 


Consequently — 
Gat; 


where c is a new arbitrary constant. 


Then 
b= V+iP=+[44+16], 
where 
A =75V Ve +A +e? 
and 


Then the solution of (6) is 
V, =A dlick+ 41+ 16) a + Belicé—(4+1)n), (7) 


Particular Solution 


In the process of extracting the particular solution, equations (1) can be inte- 
grated to 
ih 
V20, = - wv, — Ehy T’) 
and 


V4 wy + VIO, = Fe. 
A simple form of particular stress function that satisfies both equations above is 
®,=p,k n?/2 


and the particular normal displacement function becomes 


Wy = Pr R/EK + Rly T). 
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General Solution 


In the spirit of inspiration that the transverse bending of thin translational shells 
is restricted to a relatively narrow zone around its boundary, only that part of 
solution (7) is retained which decreases with increasing coordinate normal to the 
boundary under consideration. 

Enforcement of this concept in solution (7) yields an approximate function 


Vi=w, +i ®, =e“4"{(C, + iC,) (cos On — isin@y)} sinc €. 


After decomposing this solution into real and imaginary parts and then equating 
the real and imaginary parts on both sides of the equation, both the complementary 


stress function 
E h2 


V12 (1 — 22) 
and the complementary normal displacement function 
w, =e “"[C, cos On + C,sin Oy] sinc é 


c 


e-47(C,cosO 7 — C,sin@y]}sincé 
q] 4, 


are obtained. 
The general solution is the sum of complementary function and particular solution, 
hence 


E he ae haps Pn Fn 8 
: ee at [C, cos 0 yn — C, sin O yn] sinc € + 5 (8) 
and 
2 
w =e “"[C, cos On + Cz sin Oy] since € + ae +R(yT). (9) 


Stress Resultants and Stress Couples 


Finally all sectional quantities can be expressed in absolute coordinates £ and 7: 


Nes 

oe 

Qe = —D lew), 

Q, =— Dz Vrw)), 

Me = Dae vel (10) 
My —D [age + 9 ae 

My, = My = —D(L—») [see , 

Few Q: + a 

PF, ~Q, + “ae 


F; and F, designate approximate reactions acting along the edges € = 0 and 
7 = 0 respectively. 
Boundary Conditions 
As dictated by the limitations of the imposed conditions on the behavior of the 
thin translational shell all boundary conditions involving stress resultants, that 
act in the middle surface of the shell, are satisfied by the membrane solution not 
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covered in this paper*. The boundary conditions which involve stress couples, edge 
rotations and edge displacements are to be satisfied by the transverse bending solution 
presented in this paper. The assumption is made that in general the effect of membrane 
stresses on transverse bending and the effect of transverse bending on the membrane 
state of stress is insignificant. 

When it comes to the imposition of proper mathematical boundary conditions 
on shells of translation the stress analyst finds himself in a rather uneasy quandary. 
The standard practice in the design of supporting members for such shells has been 
to provide them with boundary arches that posses relatively small lateral restraining 
capacity. Yet in the proximity of the corners 
the capacity of the boundary restraint 
increases considerably. 


Fig. 2. Geometry of a Translational Surface Fig. 3. Twin Translational Thin Shallow Shell’s 
over an Almost Square Base 


A consequence of this edge behavior manifests itself in boundary conditions that 
are functions of the arch coordinate. 

An attempt to specify and satisfy such boundary conditions would be extremely 
complicated and practically untenable. Observations on the behavior of actual 
translational shell structures have indicated that the boundary conditions exhibit 
an almost constant characteristic all along the length of boundary arches and that 
marked variations appear only in the immediate vicinity of the corner of the shell. 
Therefore for practical purposes constant boundary conditions are considered adequate 
to describe the average condition of the shell’s edges. 

The following complete boundary conditions seem reasonable for the thin shell 
of translation shown in Fig. 2: 


At = 0, Wh: w=), (&) 


Meg=0,  (b) 
Ne re 0, (c) 
eae (d) 
At 7 = 0: w= 0; (6) Or 
M sq Tees (f) 
Nin=9, (g) 
& — 0, (h) 


where ¢, and ¢; designate strains along the coordinates 7 and & respectively. 


* For membrane (momentless) solution see 1, 2, 6, 14, 17, 18 and 19. 


Transverse Bending of Thin Shallow Shells of Translation 273 


In case of shells shown in Fig. 3 the boundary conditions have to be modified 
as follows: 


arié=0;L° vw 0; (a) 
Mz=0,  (b) 
Neg =0, (¢) 
é, = 0. ~ (d) 
At y=0: w=0, (e) (12) 
dw 
‘dn apes (f) 
Nin = 9, (g) 
ee = 0. (h) 


Consistent with the basic assumptions, made at the outset of this treatment, 
boundary conditions (12c) and (12g) are handled by the membrane 
solution? ?, ® 14, 17, 18, 19, 21, Conditions (12d) and (12h) are satisfied rigorously only 
when exact solutions, corresponding to functional satisfaction of equations (4) and (5), 
are applied to the shell problem. 


Hence in the present treatment of edge restraint only the boundary conditions (11 a), 
(11 b), (lle), (11 f), (12a), (12 b), (12 e) and (12f) are pertinent. 


A consideration of boundary conditions (11 a), (11 b), (12a) and (12 b) requires 
c= (2k+1)a/L. 


For the further analysis particular solutions are expanded into Fourier sine series 
in variable €. Solutions (8) and (9) can be expressed by 


Eh <a ae 
y12(1—v) =a 4,7[C,? cos 0, — C,) sin 0,4] + 


k=0 
2R 540 Bike tos) 
+ Pn Qk+ hx “| ae aed 2) 
and 
w= Sle [(C,1 cos 0,7 + C2 sin 0,7] + 
k=0 
4 py RB |. @kt la 
tees HF + Re D)|psin iF i (14) 
where 


Boundary conditions (11 e) and (11 f) are satisfied by 


4 , R? 
Cf = ok be: | Riy 1) ae 


and - 2k+1 R2 
og = — ACY Pe + Ry TD). (16) 
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Boundary conditions (12e) and (12f) are satisfied by 


4 Foi LE 
Che ee ee) fer + Ry) a 
and 
4 A eee 


For any particular translational shell the integration Onetaae CA, C,? are 
calculated from (15), (16) or (15), (17) depending upon the nature of its edge support. 
These values are substituted in (14), which represents the normal displacement of 
the translational shell brought about by the boundary restraints. 


As stipulated by the basic assumptions no use is made here of the stress 
function (13). The stress resultants that act in the plane of the shells surface are 
all given by the membrane solution. The limited effect perpetrated by the transverse 
bending upon these stress resultants is neglected. 


All the sectional quantities in (10), that depend on normal displacement function w, 
are easily computed by means of solution (14). 


Application 


For a thin reinforced concrete shell of translation (Fig. 2), enclosing a square 
base with sides B, the following design data is given: 


R= 100 ft,, 
B= 105i 
h — 0.5 tt., 
y = 0.16, 
Pn = — 9.15 k. p. ft.2 (k. = kilo pound = 1000 lbs.), 
B == 100%. 
If the normal displacement function (14) is limited to four term series, then 
k=3 
w= a, jen" [C3 cos 0,7 + C,2 sin O,7] 4 Qk + We oe oe E= 


a4 { e— 4001 cos Oy n + C,? sin Oy 4] + — 


ae We GE= 
{f —anga @) Oh say 4 Pn fr ™ 

ae ye [C,* cos On + Cy’ sin 0; 4] + 3 

+] e “"[C,! cos 0,7 + C,? sin O,y] + oe Be 


ei {e-* "IC; 1 cos O,7-+ C,? sin O, n] += — Te = [Pea |! sin &. 


The stress couple acting normal to the shell’s boundary is 


b=3 
2k + 1)? 7? 
M,,= — D S {e+ l02 a v) — ~~ cos 0,7 + #sin O, n| “e 
k=0 


LL? 
4(2k+1 ae Be 
(2k +1) % Pn \ sin (2k + 1)2 


+¢2 |— 2? cos O;, + (1 — 2) (2% + 1)? 2? ain xn eae 


L Eh oR g. (18) 
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The boundary reaction is given by 


k=3 
[ 
F,=D S) e4:"\C2|(4,3—3.4, O2 — (2 — ») A, SA) cos 0,7 — 
k=0 


—(62—3 426, (ge 5) 0, = ™ ) sin 6, 7] a3 


2k+ 1)? 72 
Oa gir 


2 Deira : 
++ (4! — 3A,0, —(2—») A, + Le ) sin Ox] } sin ee (19) 
The supporting arches of the thin shell are assumed to offer no resistance against 
lateral rotation. 
The proper integration constants are obtained from (15) and (16), hence 


C2 = 3820/((2k +1)H and C2 = (50.5) (2k + 1)/E. 


Substitution of these constants 7 in fF. 
in (18) yields the stress couple 0 
M,,. The calculated stress couple 
distribution along the coordinate 
Hhaebe fox Li 2.0ee = L/4: and: 
&=L/8 is depicted in Fig. 4. 

If more terms of the series 
were used in the solution, the 
variation of stress couples with 
the coordinate € would have been 
further reduced. 

This example demonstrates 
that the boundary perturbation 
can be treated approximately as 
being independent of the coordi- 
nate parallel to the edge of the 
shell under consideration, except 
in the immediate proximity of 
the intersecting arches. It is Fig. 4. Diagram of Stress Couple My, 
also evident from the stress 
couple diagram that the initial assumption of a narrow transverse bending 
zone around the translational shell’s boundary was a reasonable one. 

Furthermore, the boundary reaction F’,, obtainable from (19), gives that part 
of the surface load which is transmitted to the boundary supports by transverse 
stress resultant Q,. 


+ C2 |(O2 — 3.420, + (2—») 0 


Stress Couple May in ft. kips 
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Uber Stabilititsuntersuchungen mit Hilfe des Verfahrens der 
,yharmonischen Balance‘ 


Von W. Peschka, Wien 


Mit 6 Textabbildungen 


Die Stabilitat eines Systems ist durch sein Verhalten gegentiber einer beliebigen 
Stérung definiert. Eine Stérung bewirkt eine Anderung des urspriinglichen Zustandes 
des Systems. Bleibt diese Anderung beschrinkt und verschwindet sie mit wachsender 
Zeit, so spricht man von Stabilitat. Ist die bewirkte Anderung des ungestérten Zu- 
standes jedoch nicht beschrinkt, so spricht man von Instabilitat. 

Wenn nun der Betrag der Storamplitude einer Systemgrd8e sehr klein gegen den 
Wert der ungestérten SystemgréBe ist, wenn es sich also um ,,kleine Storungen“ 
handelt, dann kann die Frage, ob stabil oder instabil, relativ einfach beantwortet. 
werden. Sind jedoch die auBeren Stérungen gro8 gegeniiber den ungestérten Werten, 
handelt es sich also um ,,groBe Stérungen“, dann ist die Frage nach dem Verhalten 
des Systems ,,im groBen“ viel schwerer zu beantworten. Eine Stabilitatsuntersuchung 
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geht bekanntlich folgendermaBen vor sich: Man stellt zuerst die Gleichungen auf, 
die das Verhalten eines Systems beschreiben. Dann untersucht man mit Hilfe dieser 
Systemgleichungen, ob eine Stoérung des urspriinglichen Systemzustandes zeitlich 
beschrankt bleibt oder nicht. Damit ist dann die Frage, ob stabil oder instabil, ent- 
schieden. Beschrinkt man sich nun auf ,,kleine Storungen“, so erhalt man fiir das 
weitere Verhalten der Stérungen lineare Differentialgleichungen mit konstanten 
Koeffizienten. Dies ist eine groBe Vereinfachung gegeniiber dem allgemeinen Fall, 
weil bei einem solchen System von Differentialgleichungen das Verhalten der Losungen 
mit wachsender Zeit angegeben werden kann, ohne da8 man die Lésungen zu kennen 
braucht. Es gibt hier verschiedene Kriterien, wie z. B. das bekannte Hurwitzsche 
Kriterium, das Nyquist-Kriterium, das ,,Aufsuchen der Wurzelorte usw. (vgl.» 2). 

Das Verhalten eines Systems ,,im kleinen“ kann also mit obigen Methoden beurteilt 
werden. Hine Beurteilung des Verhaltens ,,im groBen“ ist viel schwieriger. Ist bei 
,kleinen Stérungen“ das Verhalten des Systems von der Form der Stérung unabhangig, 
so ist dies bei ,,groBen Stérungen“ nicht mehr der Fall. Man kann also beim Verhalten 
eines Systems ,,im grofen‘ streng genommen nur von Stabilitét bzw. Instabilitat 
gegeniiber einem bestimmten Stérungsverlauf sprechen. Fiir einen anderen Stérungs- 
verlauf muff dann die Stabilitaétsuntersuchung von neuem gemacht werden. 

Zu Stabilitétsuntersuchungen ,,im grofen“ ist nun das Verfahren der ,,harmoni- 
schen Balance“ sehr gut geeignet. Es mégen zuniachst die allgemeinen Gesichtspunkte 
kurz behandelt und dann einige Beispiele gebracht werden*. 

Die Gleichungen, die das Verhalten eines Systems beschreiben, kénnen in folgender 
Form geschrieben werden (vgl. * °°): 


n 
d t 
- =e = Ajy Uy + My; hi (as Xo) ee es Le (1) 


v=1 


Durch die Werte x, wird der Zustand des Systems beschrieben. Die x, konnen z. B. 
Schwingungsausschlag, Geschwindigkeit, Druck usw. sein, ¢ bedeutet die Zeit, aj, 
sind Koeffizienten, die Gré8en wu; Parameter und die Groen f; sind Funktionen, 
die mindestens quadratisch sind. (Hat f; noch einen linearen Anteil, so kann dieser 
ZU Gi x, hinzugefiigt werden.) 

Wenn nun das System eine oszillatorische Instabilitaét besitzt, dann miissen an 
der Stabilitatsgrenze anhaltende Schwingungen auftreten. Das Gleichungssystem (1) 
mu also dort eine periodische Lésung besitzen, deren Grundwelle folgendermafen 


lautet: x, = A,sin ®, = A x, sin ®, (2) 


mit ®, =at+q, ow als Frequenz, g, als Phasenwinkel und x,(m) = A,/A als mit 
einer Normamplitude A gebildeten Amplitudenverhiltnis. 
Nun wird zu (1) ein lineares Ersatzsystem gesucht, das in folgender Form ge- 


schrieben werden kann: in A 
= ot ee 1) (3) 
ate wat aii 
(Jal 
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5 K. Magnus: Uber ein Verfahren zur Untersuchung nichtlinearer Schwingungen und 
Regelvorginge. VDI-Forschungsheft 451. . 

= W. Hahn: Nichtlineare Regelvorgiinge. Vortrige der Tagung in Darmstadt 1955. Olden- 


burg. 1956. 
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“Die Koeffizienten aX und a bestimmen sich aus der Bedingung, daf die Grund- 
harmonische einer nach Voraussetzung existierenden periodischen Lésung von (1) 
fiir jede Amplitude A, mit der Lésung von (3) tibereinstimmt. 

Es sei nun vorausgesetzt, daB man die Funktionen /; in Fourier-Reihen entwickeln 
kann. Es sollen nur die Glieder, die die Grundfrequenz enthalten, beriicksichtigt 
werden. 

Die Fourier-Entwicklung von /,; lautet dann: 


f; =a, cos ®, +6, sn®,+... (4) 
mit den Koeffizienten: 
20 
an =~ | fi(Amsin ®, .. ., A x, sin B,) cos ®, d,, 
: (5) 


22 


bu =—\ f(A asin, ..., A x, sin ®,) sin , dy. 
0 


Als Grundharmonische der Funktionen /; kann dabei die Schwingung der Variablen z, 
genommen werden. Hier wurde die erste, also x, gewahlt. 
Setzt man nun (5) und (4) in (1) ein, so erhalt man: 


am 


n 
“E =D) a -A-x,-sin D, + mu; (a cos ®, + 6,, sin ,), (6) 
vy=1 
d nN 
a = >) (ah A+ %- w+ cos ®, + ay A+ x + sin ®,). (7) 


v=] 
Durch Koeffizientenvergleich erhalt man dann: 


Gs One. DNAs Goo Omurery i, 
ferner: 


=o Mi 
Giy = Uy + A uy, bas 


8 
ose it (8) 
ay = A wy Ay 
mit (5) erhalt man daraus: 
27% 
Gy =O, + ie \f:(Amsin®,,..., Ax, sin ®,)sin B, dD, 
20 . (9) 
tity = Paes: | f,(Ax,sin @,,..., Ax, sin ®,) cos DB, dQ,. 
1 


Das Verfahren der harmonischen Balance beruht also, wie schon erwaihnt, auf 
der Einfiihrung eines linearen Ersatzsystems (3) fiir das nichtlineare Ausgangs- 
system (1). 

Es bestehen jedoch zwei wesentliche Unterschiede gegeniiber der Methode der 
kleinen Schwingungen: Die Koeffizienten @ und a* sind nicht mehr konstant, sondern 
Funktionen von A. 

Ferner tritt auf der rechten Seite von Gl. (3) noch die zeitliche Ableitung der 
Verinderlichen x, auf. Dadurch wird die Anwendung des Verfahrens der harmonischen 
Balance auch bei Funktionen f;, die vom Hysteresetyp sind, méglich (vgl. 5). Das 
Verhalten der Lésungen des linearen Ersatzsystems (3) kann nun mit Hilfe der- 
selben Verfahren, die fiir die Stabilitiétsuntersuchung bei kleinen Stérungen ent- 
wickelt wurden, diskutiert werden®, 
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Stabilitatsgrenze: Zuniichst ist die Stabilitaitsgrenze selbst von Bedeutung. Man 
bendtigt bekanntlich hiezu zunichst die Sakulargleichung oder auch charakteristische 
Gleichung genannt, die aus (3) mit Hilfe des Lésungsansatzes e*% folgt: 


las + 0% -2— 6, a] = S'c,_,- 0? = 0. (10) 
»=0 
A sind hiebei die Exponenten des Lésungsansatzes — die im allgemeinen komplex 


sind —, dy, ist das Kroneckersche Symbol. Die Stabilitat der Lésungen von (3) ist 
nun gesichert, wenn die Koeffizienten c der charakteristischen Gleichung unter- 
einander gleiches Vorzeichen haben, und die » Hurwitz-Determinanten H, positiv 
sind. Ist c) = 1, was durch Division mit 4” in (10) stets erreicht werden kann, dann 
sind die Bedingungen c, > 0 und H,_, > 0 die schirfsten. Diese Bedingungen be- 
stimmen also den Stabilitatsbereich®. Dabei besitzt die charakteristische Gl. (10) 
an der Grenze H,_, = 0 ein imaginéres Wurzelpaar, an der Grenze c, = 0 eine 
Nullwurzel. Hs sind also im ersten Fall ungedimpfte Schwingungen méglich, wahrend 
im zweiten Fall das System aperiodisch aus einer gegebenen Lage abwandert. Der 
Fall c, = 0 soll jedoch aus unseren Betrachtungen ausgeschlossen werden, da in 
diesem Fall das Verfahren der harmonischen Balance aus naheliegenden Griinden 
nicht mehr angewendet werden kann. AuSerdem kommt dem Fall c, = 0 keine 
so groBe praktische Bedeutung zu, wie dem Fall H,,_, = 0, wo Schwingungen auftreten 
kénnen. Im folgenden soll nur dieser Fall weiter untersucht werden. Die Frequenz 
der Schwingungen an der Stabilitaétsgrenze kann leicht ermittelt werden? ?: >. 


re: eee n> 3, | (11) 
Ha LG = Ee] 
Fiir » <3 sind die Ergebnisse noch einfacher: 
i ig OS as 
: (12) 
See pe ee 
eu Co” 


Stabilitétsuntersuchung: Die Stabilitaétsuntersuchung linearer Systeme geht nun 
in der Weise vor sich, indem man prift, ob der Bildpunkt, der dem jeweiligen Zustand 
des Systems entspricht, innerhalb des Gebietes H,_, > 0 liegt oder nicht. Da das 
Ersatzsystem (3) im Grenzfall verschwindender Amplituden A = 0 mit dem System 
identisch wird, das man mit der Methode der kleinen Storungen erhalt, so ergibt 
sich also zunadchst derselbe Bildpunkt wie beim Verfahren der kleinen St6rungen. 
Die Koeffizienten a und a* in (3) sind jedoch Funktionen der Amplitude A und der 
Bildpunkt wird mit wachsendem A auf einer Kurve wandern, die A-Kurve bezeichnet 
wird. 

Aus der gegenseitigen Lage dieser A-Kurve und der Kurve (bzw. Flache) H,_; = 0 
ist es nun moglich, das Verhalten eines Systems im grofen zu beurteilen. Man nennt 
dieses Verfahren deswegen auch A-Kurvenverfahren. 

Die Eigenschaften eines solchen Systems, das im Parameterraum durch H,_, = 0 
und eine A-Kurve dargestellt wird, mégen nun angefiihrt werden (entnommen aus °). 
Hiebei wird der Fall A = 0 Verhalten im kleinen und der Fall A — co Verhalten 
im groBen genannt (s. Abb. 1). 


1. Ist die A-Kurve fiir alle A im Stabilitatsbereich H,,_, > 0, dann ist das System 
stabil im groBen wie im kleinen. 

2. Das System ist instabil im groBen wie im kleinen. : 

3. Das System fihrt stationiire Schwingungen mit der Amplitude A,’ aus. 

4. Das System ist im kleinen stabil, im grofen jedoch instabil. 
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5. Das System ist im kleinen stabil; erreicht die Stéramplitude den Wert A,’, 
so wandert das System nach A,” ab und fihrt dort mit dieser Amplitude stationare 
Schwingungen aus. 

6. Das System ist instabil im kleinen, fiihrt an A,’ mit dieser Amplitude unge- 
dimpfte Schwingungen aus und wird mit A =A,” instabil im groBen. : 


=e 
927 
/y-7 0 
8277 Dee eS a eee A=0 
7 
In-1=4, 2 w 
Y VY Yp 
U5 
4-129 
3 
Abb. 1 Abb. 2 
Anwendungen 


I. Rohrengenerator (entnommen aus ®*: ”) 
Die zustindige Differentialgleichung ist die von van der Pol? ’. 
Z—u(l—2)z+o%2=0. (13) 


z ist die Variable, w die Frequenz und mw Parameter. Mit z,=z und z,=2z laBt 
sich (13) in der Form (1) schreiben: 


eae 
; 14 
‘nes eae ce 
mit dem Ansatz 

4, =Asinot, 

t,=Aacoswt 


ergeben sich die Koeffizienten des Ersatzsystems: 


Maced ae * * * e 
Ay, = 9, Mg = 9, Ay = AQ =—Mg=—0, Ay =—a.=—1, 
2% 
= be é : 
Gay =O, + | (1 — A*sin? ®) A w cos Gsin Odo = a,, = — w? 
0 


22 | (15) 
ecto | 


? 


a = >a5 | (1 — A? sin? B) A @ cos ® cos ®- d® = (1 ~+). 
0 


Das lineare Ersatzsystem lautet somit: 


a oe 225 . | (16) 
Ze = — wiz + ad, 2, | 
die charakteristische Gleichung lautet: 
#2 — at, 14 +02? =0 (17) 


7 Van der Pol: On relaxation oscillations. Phil. Mag. 2 (1926). 
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a 


und die Stabilitaétsbedingungen sind: 


tg=o > 0, ce) = By = = ate 0; (18) 
die erste ist hiebei stets erfiillt. 

(18) kann in einer azi, w?-Ebene graphisch dargestellt werden (Abb. 2): Die 
A-Kurve ist eine Gerade und beginnt mit 4 = 0 im Punkte a%, — pf und @ = wy. 
Sie verlauft nach unten und trifft die Stabilitiitsgrenze bei A’ — 2, bei weiter an- 
steigendem A lauft die Kurve in das Innere des Stabilitétsbereiches hinein. Das 
System ist im kleinen instabil, im groBen stabil. Der Arbeitspunkt (@o,.A-= 2) ist 


stabil und das System hat dort einen weichen Schwingungseinsatz (mit kleinen Ampli- 
tuden beginnend) (Genaueres s. 5). 


Il. Treibstoffordersystem eines Fliissigkeitsraketenmotors 
(Monopropellant)® 


Die wichtigsten Bezeichnungen sind in Abb. 3 angefiihrt. 

p; bedeutet den Brennkammerdruck, p, den Druck unmittelbar nach der Pumpe P, 
m die Férdermenge in der Zeiteinheit. F ist der gesamte lichte Querschnitt der 
Treibstoffzuleitungen, / ihre mittlere Linge. Die Gleichungen fiir dieses System lauten 
dann in dimensionsloser Form: 


d 
en? Ets 
d 1 1 (19) 
oo alt Owe 
mit: 
_. Pk=Pr egpllaien tn my eoee | 
Pr ? Na i‘ ? o,2 (20) 
sce ae ele dn 


p, und m bedeuten die Stationirwerte von p;, und m, t ist die Zeit, 9, ist eine beliebige 
Bezugszeit (hier die mittlere Laufzeit eines Gaselements in der Brennkammer). f(m) 
ist die stationire Pumpencharakteristik bei konstanter Drehzahl. 

Das lineare Ersatzsystem lautet: 


dy 


a art aay 
d 
a te Seis, Hae |. | 


f*(1 +) wird nun an der Stelle » = 0 in eine Potenzreihe bis zu Gliedern vierter 
Ordnung entwickelt: 


fl + pu) —L= ayy +a, + agus + Ogu +... (22) 
Der Ersatzkoeffizient a,. lautet dann mit Hilfe von Gl. (9): 
27 ! 
hee dq, = = . — [a, A? sin? ® + a, A?sin? © + a, A*sin*®]sin@d®, (23) 
a 
0 
dag = Lay Got AD (24) 


, G. Heinrich und W. Peschka: Uber Stabilitatsuntersuchungen an flissigkeitsgetriebenen 
Raketenmotoren mit Hilfe des Verfahrens der Harmonischen Balance. Vorgetragen am 
VIII. Intern. Astronautischen KongreB in Barcelona (6. bis 12. Okt. 1957). KongreB8bericht. 
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mit dem Lésungsansatz 


@ = Aiet*) Tea ete (25) 
erhaélt man aus (21): 
A(A+1)—B=0, 
(26) 


A+— + (A—ay) B= 0. f 


Eine Lésung des homogenen Gleichungssystems (26) existiert nur dann, wenn seine 
Koeffizientendeterminante verschwindet. 


i 
M+ (1 — G9) A — ( i =) =a, (27) 
Daraus ergibt sich A. 


Stabilitaét ist nun vorhanden, wenn die Lésung (8) beschrankt ist, wenn also der 
Realteil von A, oder / selbst, negativ ist. Daraus ergibt sich die Stabilitatsbedingung: 


Ao <1. : (28) 


Das System wird aperiodisch, wenn gilt: 


1—a,, \2 , / 1 
Geran) (29) 

4 222 
es 

ope 
2 

= L=44 
OP; : 
Ar lastabif oes 
Tb ae = = 
Cli VU dddg Ze U4 
Ar tabi! 
stabi 

: == Po eee 

A 4 5 6 ere 

pi 
-7 7 
2=0 
2. 
R=O5 

oc) Lat 

= 92=7 
-4 Z 2=4,4 


Abb. 4 


In Abb. 4 sind die Ergebnisse (28) und (29) in einem Ao, d-Diagramm aufgetragen. 
Das Gebiet unterhalb a,, = 1 ist also stabil. Ferner sind die Kurven, auf denen die 
Frequenz* 2 konstant ist, aufgetragen. Q ist einfach der Imaginarteil von A. 


or = (“S%)'+ (ay, — 2), (30) 


* 2 =o %8,; » Kreisfrequenz. 
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In dem durch die Kurven 2 = 0 begrenzten Gebiet kénnen keine Schwingungen 
auftreten, wahrend sonst iiberall Schwingungen mdglich sind. 

Zuletzt mu noch a,, durch a und A ausgedriickt werden, was mit Hilfe von (24) 
bereits getan wurde. 


3 
Oy, =~ 0, +2 A. (24) 


Der Arbeitspunkt des Systems liegt also fiir A = 0 an der Stelle (asa, a1): 
Mit wachsendem A wird er sich auf einer Kurve, der A-Kurve, bewegen. Die A-Kurve 


ist hier eine Gerade parallel zur Achse a,,. 
Damit ist nun eine einfache Diskussion aller vorkommenden Falle méglich (Abb. 4). 


it 0 gO 

a) Der Arbeitspunkt ist P}. 

Bei Zunahme von A geht a. —> — oo, bleibt also im Stabilitatsbereich. Das 
System ist im kleinen wie im groBen stabil. 

b) Der Arbeitspunkt ist Pi oder P3!. 

Das System ist im kleinen stabil. Im groBen fiihrt es Schwingungen an den 
Punkten A} oder A}’ aus. Frequenz und Amplitude kann aus dem Diagramm 
enthnommen werden. 

Bea SOs, = 0: 


Ist a ee dann ist das System im kleinen stabil, im groBen aber instabil, 


ist ~-a, > 1, dann ist das System im kleinen wie im groSen instabil. 
dope 0. a0, ; 
Ist 2 a, <1, dann ist das System im -kleinen wie im grofBen stabil. Ist — a, > 1, 


dann tritt Fall 2 auf. 
4. a, <0, a, > 0: Das System ist im kleinen stabil, im groBen jedoch instabil. 


= 


Das Stabilitaétsdiagramm ist aber auch fiir Untersuchungen im kleinen geeignet, 
also fiir den Fall A — 0. Die Stabilitatsbedingung lautet hier einfach nach (28): 


mI 
By, FO, 
Wird der Parameter //F der Rohrleitung klein, dann 
reduziert sich im Falle 1/F = 0 GI. (31) auf die Be- 
dingung: 


i 
aaa < Le OGOF ». x, <0. = (31) 


eral (32) 


Diese Bedingung ist von den Kreiselpumpen und 
Axialverdichtern her bekannt, wonach ein stabiler 
Betrieb im schraffierten Bereich in Abb. 5 nicht méglich 
ist. Aus Gl. (31) sieht man jedoch, da die Rohrleitung 
stabilisierend wirkt. 

Um also stabilen Betrieb zu erhalten, mu8 der Arbeitspunkt erstens einmal in 
einem Bereich liegen, wo a, <a. Dann ist das System im kleinen stabil (Fall 1 und 4). 
Nun hangt es ganz von a;, also im wesentlichen von f*’” (1), der dritten Ableitung 
der Pumpencharakteristik, ab, ob das System auch im grofen stabil wird. Stabilitat 
im groBen ist vorhanden, wenn a, < 0. Ist jedoch a, > 0, dann ist das System im 
groBen instabil. 

Dieser Fall ist von besonderem Interesse: Wenn die Kennlinien der Treibstoff- 
pumpen Parabeln sind — wie es bei Radialpumpen meist angenommen wird —, 
dann ist /*’’’ = 0 und das System verhilt sich dann im grofen ebenso wie im kleinen.- 

20% 
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Im allgemeinen treten aber stets Abweichungen von diesem Verhalten auf, und be- 
sonders bei Grenzleistungsverdichtern, wie sie gerade im Raketenbau vorkommen, 
ist es méglich, daB die Kennlinien Wendetangenten besitzen, und also zumindest 
die dritte Ableitung von Null verschieden ist. 

Interessant ist iibrigens, daB die geradzahligen Ableitungen der Kennlinie keinen 
direkten Einflu8 auf die Stabilitat besitzen. 

Es kénnen aber auch Abweichungen vom Parabelverlauf der Kennlinie auftreten, 
wenn die dynamische Kennlinie mit der statischen nicht mehr tibereinstimmt. Des- 
wegen wiire es nétig, dem Verlauf der dynamischen Kennlinie mehr Aufmerksamkeit 
zu widmen. 


Ill. Flissigkeitsraketenmotor mit Totzeit in der Verbrennung und 
konstanter Treibstoffo6rderung! ® § 


Die Gleichungen des System reduzieren sich hier auf eine einzige (vgl. °). 
dp 1 + v2) __|” 
mit 6 = = [s. Gl. (20)]. 7 ist hiebei die Totzeit. n ist ein Parameter. 


Wenn der Treibstoff in die Brennkammer eingebracht wird, so bendtigt er noch 
eine gewisse Zeit (t), bis er ein brennfahiges Gemisch bildet und bis zur Zindtemperatur 
erwirmt ist. Dann setzt erst die Verbrennung ein. Diese Totzeit + kann nun Ursache 
einer oszillatorischen Instabilitét eines Raketenmotors sein, wie zuerst L. Crocco 
gezeigt hat®. Diese Untersuchungen wurden von ihm jedoch nur mit der Methode 
der kleinen Schwingungen durchgefiihrt. In*® wurden diese Untersuchungen ver- 
allgemeinert, indem auch grofe Stérungen betrachtet wurden. Das Verfahren der 
harmonischen Balance hat sich hiebei sehr bewahrt. 

Das lineare Ersatzsystem zu (33) lautet: 


d Si: a 
a + (1 — 1) ¢(2) + 2-9 (2 — 6) =a | +4 9 (2), (34) 
die Ersatzkoeffizienten ergeben sich nach langerer Rechnung: 


ete = {(n — 1) (n — 2) — 2(2n? — 1)cosQ26+n(n 4+ 1) (14 bps 


ae (35) 
a1 = Zo (2 ( +n + 1)sinQ 6 — (n + 1) nsin 22 6}. é | 


Da sowohl in (33) als auch in (34) nicht nur die Funktionswerte an der Stelle z, 
sondern auch an der Stelle z — 6 vorkommen, ist dies keine gewdhnliche Differential- 
gleichung. Derartige Gleichungen nennt man Hysterodifferentialgleichungen”. 

Gl. (34) geht in die Gleichung der linearisierten Theorie® itiber, wenn man A und 
damit aj, und a,, gegen Null gehen lat. Mit Hilfe von (34) und (35) kann nun 
wieder ein Stabilitétsdiagramm gefunden werden. Die numerischen Berechnungen 
erfordern jedoch einen enormen Aufwand an Arbeitszeit. Die Berechnung wiirde sich 
aber bei Verwendung einer programmgesteuerten elektronischen Rechenmaschine in 
relativ kurzer Zeit ausfiihren lassen. 

Das Stabilitatsdiagramm fiir den Fall » = 0°75 zeigt Abb. 6. Eine Diskussion 
der vorkommenden Fille ist damit leicht durchfiihrbar. Man sieht, daB z. B. fir 
6 = 0°25 das System instabil ist, und da die Kurve 6 = 0°25 die Stabilitatsgrenze 


* L. Crocco: Aspects on Combustion Instability in Liquid Propellant Rocket Motors. Jet 
Prop. 1951/52, 6/1. 


<i N. Minorsky: Self Excited Oscillations in Dynamical Systems Possessing Retardet 
Actions. J. Appl. Mech. 9, Nr. 2 (1942). 
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schneidet, wird dort das System mit 
der zugehérigen Amplitude A Schwin- 4 
gungen ausftihren, weil diese Kurve 

mit wachsendem A aus dem instabilen 

in den stabilen Bereich eintritt. zs 

Die Frequenz kann an der Q 6-Achse 
abgelesen werden. slaPh 

Interessant ist, daB die Kurven 
6 < 0°25 die 2 6-Achse nicht erreichen, 
das hei8t, daB der Wert A =0 im 
groBen in diesen Fallen nicht méglich 7% 
ist. Dies deutet auf einen harten 
Schwingungseinsatz hin. Es kann nim- 
lich das System im kleinen stabil sein, 
aber bei einer endlichen Stérung setzt 6-01 
plotzlich eine Schwingung von end- 
licher Amplitude ein. 

Die Annahme konstanter Treibstoff- 
forderung ist allerdings eine starke 
Idealisierung. Auch die Elastizitét der 
Rohrleitungen, die hier vernachlissigt 
wurde, kann noch Korrekturen zu 
obigen Ergebnissen liefern. 


ANNAN ARR entcsseseeserenss 


~ 


QAO AA 


AYE 


stabil 
stabil 


NISSEN ESS CO SEES SES CE SE ESE eee 


re 
y 
| 


26 
Q 7 2 
Stabil. Ben fA=0 | 


Nachtrag 


In letzter Zeit ist eine weitere Arbeit Abb. 6 
von K. Magnus in der Z. f. angew. 2m 
Math. Mech. (11/12, 1957) unter dem Titel ,,Uber den Zusammenhang verschiedener 
Naherungsverfahren zur Berechnung nichtlinearer Schwingungen“ erschienen. Darin 
wird gezeigt, daB das Verfahren der Harmonischen Balance sowie auch andere Naherungs- 
methoden zur Berechnung nichtlinearer Schwingungen als Sonderfalle des allgemeinen 
Ritz-Galerkin-Verfahrens aufzufassen sind. Unter diesen Gesichtspunkten sind dann 


sogar noch allgemeinere Aussagen modglich. (Hingegangen am 21. November 1957) 


Mitteilung aus dem Institut fiir Einfiihrung in die Stark- und Schwachstromtechnik, Wien 


Uber das Drehmoment im elektromagnetischen Feld 
Von H. Hofmann 


Erster Teil 
Elektrostatisches Feld 


Mit 6 Textabbildungen 


Zusammenfassung. Die in vielen neueren Arbeiten behandelte Frage des Drehmomentes 
im elektromagnetischen Feld, speziell im Falle anisotroper Materie, wird einer grundlegenden 
Untersuchung unterzogen. Es wird gezeigt, daB wegen der Polaritat der Materie im elektro- 
magnetischen Feld im Gegensatz zur Mechanik des Kontinuums unsymmetrische Spannungs- 
tensoren auftreten. Unter anderem erweist es sich auch, daB die Maxwellschen Spannungen 
bei anisotroper Materie ihre Giltigkeit behalten, wahrend zum Moment der Maxwellschen 
ponderomotorischen Volumskrifte in diesem Fall noch ein Zusatzglied zu treten hat: Ks ist das 
vektorielle Produkt aus Polarisations- und Feldvektor und stellt die Momentendichte der im 
Materieinneren an den polarisierten Volumselementen angreifenden Kraftepaare dar. 
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I. Vorbemerkungen iiber die Schreibweise von Tensoren (Vektoren) und Formelzeichen 


Im Zuge der Angleichung an bestehende, insbesondere an inzwischen neu er- 
schienene Normen erfolgt eine geringfiigige Abiinderung der Schreibweise gegentiber 
derjenigen, die in den bisherigen Veroffentlichungen des Verfassers angewendet wurde. 


In der symbolischen Schreibweise der Tensoren (Vektoren) bedeutet jetzt 


AB statt (2B) das skalare, 
Ax B statt [2B] das vektorielle, 


%:B statt WB das dyadische Produkt der beiden Vektoren Xf und B (s. ONORM 
A 6411, DIN 1303, Ausgabe Juli 1946, und DIN 1303, Neuentwurf August 1957). 
Die Zusammenfassung von Summanden und Produkten erfolgt nur noch mit Hilfe 
der runden oder der geschlungenen Klammern, das Komma wird hierfiir nicht mehr 
verwendet, also 


¢ x (B +C) statt (2, B+], 
A{B(C;D)} statt (2(B,C D)) usw. 


Bei den Formelzeichen wird die Bezeichnung ‘)t fiir den magnetischen Polarisations- 
vektor in § abgeandert und fiir die Stromdichte der GroBbuchstabe & an Stelle des 
Kleinbuchstabens g verwendet (s. ONORM A 6401). . 


II. Der Kraftangriff des elektromagnetischen Feldes an der Materie 


Die Differentialgleichungen fiir das elektromagnetische Feld beinhalten unmittel- 
bar keine Aussage iiber den Kraftangriff des Feldes an der Materie und man muf 
daher zusitzlich irgendeine Annahme oder ein als allgemeingiiltig angesehenes Prinzip 
(Energieprinzip, Impulssatz oder anderes) heranziehen, um zu einer derartigen Aus- 
sage zu gelangen. Je nach der Art dieser Annahme erhalt man verschiedene Aus- 
driicke fiir Kraftdichten, denen jeweils eine ganz spezielle physikalische Bedeutung 
zukommt. Eigenartigerweise wurden diese verschiedenen Bedeutungen bislang nicht 
immer entsprechend gewiirdigt, so da die praktische Anwendung der einzelnen 
Kraftdichten auf sich scheinbar widersprechende Resultate gefiithrt hat. So errechnen 
Sommerfeld und Bopp in einer neueren Arbeit! fiir das Drehmoment auf einen 
Magneten, der sich in einer permeablen Fliissigkeit befindet, zwei verschiedene 
Resultate, die sich wie pu, :1 (u,... relative Permeabilitét der den Magneten um- 
gebenden Fliissigkeit) verhalten. Sie empfehlen, ohne weiteren Versuch einer 
theoretischen Klaérung, eine Entscheidung durch das Experiment. Uber die Kraft 
oder das Drehmoment liegen unter anderem weitere neueste Arbeiten von Diessel- 
horst?»3, Déring*:> und Sommerfeld und Ramberg® vor. Uber Unklarheiten im 
Zusammenhang mit den Kraft-, Energie- und ImpulsgréBen des elektromagnetischen 
Feldes ‘diskutiert ausfihrlich auch v. Laue in der neuesten Auflage seines Buches 
tiber die Relativitatstheorie’ (1. Bd., § 19). 


' A. Sommerfeld und F. Bopp: Zum Problem der Maxwellschen Spannungen. Ann. 
Physik 8, H 1/2 (1950). 

* H. Diesselhorst: Feldwirkungen in permeablen Medien und Definition des magnetischen 
Momentes. Ann. Physik 3 (1948). 

* H. Diesselhorst: Magnetfeld und Drehmoment bei einem magnetischen Ellipsoid in 
permeablem Medium und Fremdfeld. Ann. Physik 9, H. 6/7 (1951). 

a We Doring: Zur Definition des magnetischen Momentes. Ann. Physik 6 (1949). 

* W. Doring: Uber die Kraft und das Drehmoment auf magnetisierte Kérper im Magnet- 
feld. Ann. Physik 9, H 6/7 (1951). 

U A. Sommerfeld und EK. Ramberg: Das Drehmoment eines permanenten Magneten im 
Felde eines permeablen Mediums. Ann. Physik 8, H. 1/2 (1950). 

7 M. v. Laue: Die Relativitétstheorie, Bd. I, 6. Aufl. Vieweg u. Sohn. 1955. 
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Alle erwaihnten Arbeiten sowie auch diejenigen des Verfassers behandeln das 
Problem des Kraftangriffes vom makrophysikalischen Standpunkt aus, das heiBt 
unter Verwendung der GréBen des Makrofeldes, wobei man sich diese allerdings durch 
Mittelwertsbildung aus den GréBen des Mikrobereiches entstanden denken kann. 
Die Ergebnisse setzen also die Annahme voraus, da die Materie als Kontinuum 
autgefaft werden kann und sind damit, in gleicher Weise wie auch die gesamte Maxwell- 
sche Theorie, auf den Bereich der kontinuumsmaBigen Darstellung beschrinkt. Die 
praktische Anwendung insbesondere in der Elektrotechnik mu8 ihr Auslangen im 
Rahmen der Kontinuumstheorie mit den GréBen des Makrofeldes finden und findet 
sie auch, doch sei an dieser Stelle darauf verwiesen, daB fiir Untersuchungen, die 
eine mikrophysikalische Behandlung benétigen, die Rechnung unter Verwendung 
der GréBen des Makrofeldes bekanntlich nur eine mehr oder weniger grobe Naherung 
darstellen kann, mitunter versagen die Aussagen der Kontinuumstheorie auch vollends. 
Es wird im folgenden anlaBlich der Besprechung der Polarisationsvektoren 8 und $ 
nochmals auf die Bedeutung der kontinuumsméBigen Darstellung und ihre Grenzen 
eingegangen. Auch in der Mechanik arbeitet man vorwiegend und mit Erfolg 
kontinuumstheoretisch und die Parallelen zwischen dem Kraftangriff am mechanischen 
Kontinuum, speziell im Rahmen der Elastizitatstheorie, und an der kontinuierlich 
gedachten Materie im elektromagnetischen Feld treten immer wieder in Erscheinung. 
Kin Unterschied ist allerdings evident: In der Mechanik des Kontinuums gibt es 
nur homogen mit Masse erfiillte Volumselemente, die sich allerdings nicht immer 
isotrop verhalten miissen. Im elektromagnetischen Feld erweist sich hingegen ein 
Volumselement kontinuierlicher Materie nicht nur als homogen mit (,,wahrer‘‘) 
Ladung und Leitungsstrom erfiillt, sondern es ist auSerdem noch polarisiert. In 
Zusammenhang mit dieser Polaritét, einer Erscheinung, die eben in der Mechanik 
des Kontinuums nicht auftritt, steht unter anderem die Tatsache, daB im Gegensatz 
etwa zum symmetrischen Spannungstensor der Elastizitatstheorie im _ elektro- 
magnetischen Feld auch unsymmetrische Tensoren auftreten. Das polare Verhalten 
der Materie im elektromagnetischen Feld ist gleichsam der Schliissel fiir das Ver- 
stindnis dafiir, da die Parallelen mit der Mechanik des Kontinuums nicht voll- 
kommen sein kénnen, und es laBt sich zeigen, daB die immer wieder erhobene Forderung 
nach ausschlieBlich symmetrischen Tensoren im elektromagnetischen Feld vollkommen 
unberechtigt ist. 

Die Behandlung des Kraftangriffes des elektromagnetischen Feldes an der Materie 
durch den Verfasser in *: % 1° 11 brachte verschiedene prinzipielle Ergebnisse, die auch 
fiir die Berechnung des Drehmomentes von Bedeutung sind. Sie sollen zunachst 
zusammengefabt dargestellt werden: 


A. Der Unterschied zwischen bewegender (ponderomotorischer) und 
direkter (kérperfester) Kraft 


Die Kraft, mit der ein Feld direkt (korperfest) an dem punktformigen Trager 
einer Menge angreift, ist durch das Produkt ,,Menge mal Feldstarke“ definiert, also 
z. B. im elektrischen Feld durch das Produkt ,,Ladung mal elektrische Feldstairke“: 

& =O. (1) 

8 H. Hofmann: Uber den Kraftangriff des stationaren elektromagnetischen Feldes an der 
Materie. Osterr. Ingenieur-Arch. 10, H. 4 (1956). ; ‘ 

9 H. Hofmann: Uber den Kraftangriff des Magnetfeldes an Elementarstrémen. Osterr. 
Ingenieur-Arch. 11, H. 1 (1957). ae, 

10 H. Hofmann: Die Behandlung von Sprungflachen beim Kraftangriff des stationéren 
elektromagnetischen Feldes an der Materie. Osterr. Ingenieur-Arch. 11, H. 1 (1957). 

1H. Hofmann: Uber den Kraftangriff des allgemeinen, nichtstationaren elektromagneti- 
schen Feldes an der Materie. Osterr. Ingenieur-Arch. 11, H. 2 (1957). 
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Man wird vom réiumlichen Ausmaf des Tragers der Menge unabhangig, wenn man 
zu dem allgemeineren Ausdruck mit Hilfe der DichtegroBen iibergeht, im elektrischen 
Feld somit auf 

hao’ (2) 
(| ... Kraftdichte oder Volumskraft, @ ... Ladungsdichte). Die gesamte am K6rper 
in seinem Volumen vom Felde angreifende Kraft ist dann durch das Integral 


% =| fdV =\eEdV (3) 
v Vv 


gegeben, wobei hier iiber die néhere Bedeutung der GréBen Ladungsdichte @ und 
Feldstairke © zunachst noch nicht diskutiert werden soll. 

Die meBbare, den Kérper bewegende (wigbare, ponderomotorische) Kraft ist 
im allgemeinen Fall jedoch nicht identisch mit dieser am Korper in seinem Volumen 
direkt angreifenden Kraft des Feldes, da der Kérper auch an seiner Oberflache durch 
die umgebende Materie einwirkenden Kriaften unterliegt. Sie werden ebenfalls durch 
das betreffende Feld, namlich durch Einwirkung auf diese umgebende Materie, 
hervorgerufen. Im Gravitationsfeld ist fiir die Resultierende 8 dieser auf die Ober- 
flache A des Korpers ausgeiibten Krafte 


B= pda (4) 
A 


(p ... Flachenkraft) der Ausdruck ,,Auftrieb“ gebrauchlich, der auch hier im elektro- 
magnetischen Feld ohne weiteres im selben Sinn verwendet werden kann und soll. 
Die bewegende, ponderomotorische Kraft %yonq setzt sich also aus der k6rperfesten, 
im Volumen des Korpers direkt angreifenden Kraft 7%} und dem vom Feld indirekt 
durch die umgebende Materie hervorgerufenen, an der Oberflache angreifenden Auf- 
trieb J gemafB 


Spa =F + B=|fdV +opa (5) 

zusammen. us “ 
Ks gelingt nun — sowohl im (stationéren) Gravitationsfeld als auch im (stationaren) 
elektromagnetischen Feld — die ponderomotorische Kraft Gl. (5) in ein Volums- 


integral tiber den betrachteten Korper allein umzuwandeln (s. #2, Abschn. II), dessen 
Integrand f,onqg dann formal gréBenmaBig die Bedeutung einer Kraftdichte besitzt, 
physikalisch jedoch keineswegs als die an jedem Volumselement (je Volumseinheit) 
vom Felde tatsaichlich angreifende Kraft gedeutet werden darf: 


pond = 8 + B= [FAV +h pdd =| fomadV. (6) 
Diese formale ponderomotorische Volumskraft f, nq hat lediglich die Bedeutung, 
daB sie allein tiber das Volumen des K6rpers integriert unmittelbar bereits die gesamte 
bewegende Kraft %,ona, die sich physikalisch gesehen ja aus der direkten, im Volumen 
angreifenden Kraft und dem Auftrieb % gema&B Gl. (6) zusammensetzt, liefert. 
Man kann zur ponderomotorischen Kraftdichte auch tiber Energiebetrachtungen 
kommen: Es gelingt namlich, die bei einer virtuellen Verschiebung 63 verbrauchte 
Feldenergie — dW, die nach dem Energieprinzip gleich der gewonnenen mechanischen 
Arbeit 6A sein mu8, unter gewissen Voraussetzungen in ein Volumsintegral iiber den 
betrachteten Korper allein umzuwandeln: 


~ OW. =:+..04 =5\ fone ser. (7) 
4 
Die auf diesem Wege erhaltene Ableitung zeigt, daB der mit der virtuellen Verschie- 


bung 63 multiplizierte Ausdruck f,.,q eben wiederum nur einen formalen, mit dem 


12 §. FuBnote 8, S. 287. 
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Energieprinzip vertraglichen Ausdruck fiir eine Kraftdichte darstellt, der die Wirkung 
der physikalisch deutbaren, im Volumen wirksamen direkten Kraftdichte f und der 
an der Oberflache angreifenden Flichenkraft p in sich vereinigt. 

In gleicher Weise mu8 zwischen dem Moment der kérperfesten und dem Moment 
der ponderomotorischen Krafte unterschieden werden. ZweckmiBigerweise iiber- 
nimmt man dann die analogen Bezeichnungen ,,kérperfestes (direktes) Drehmoment 
und ,,ponderomotorisches (bewegendes) Drehmoment‘‘. Das kérperfeste Drehmoment 
muB, vermehrt um das Moment der an der Oberfliche des betrachteten Korpers 
angreifenden Auftriebskrafte, wiederum das meBbare ponderomotorische Moment 
liefern. 

Befindet sich ein Ko6rper im Vakuum, so wird mangels an Auftriebskriften die 
Integration der kdrperfesten und der ponderomotorischen Volumskrafte unmittelbar 
auf dasselbe Resultat fiihren, was in  (Abschn. VIII) gezeigt wurde. Dieselbe Uber- 
einstimmung besteht in einem solchen Fall auch fiir das kérperfeste und das pondero- 
motorische Drehmoment, was noch spater im Rahmen dieser Arbeit nachgewiesen 
werden soll. 

Im elektromagnetischen Feld ist demnach streng zwischen den etwa aus dem 
Lorentz-Ansatz fiir die Kraftdichte im Mikrobereich . 


1 
t — 0) : e(2) + ‘5 oy x h® 14 (8) 


((, § ... elektrische bzw. magnetische Feldstarke im Mikrobereich nach der 
Elektronentheorie) abgeleiteten k6rperfesten Volumskraften des Makrofeldes, die 
auf Produkte aus Mengendichten und Feldstarken zuriickgehen, und den Maxwellschen 
ponderomotorischen Kraften zu unterscheiden. Der Lorentz-Ansatz Gl. (8) beschreibt 
ja nur den direkten Kraftangriff des Feldes an der Ladungsdichte 0, die aus dem 
Energieprinzip gewonnenen Maxwellschen ponderomotorischen Volumskrafte ent- 
halten hingegen auch den Kraftangriff der den Korper umgebenden Materie, also 
den Einflu8 der Auftriebskrafte. In (Abschn. VIII, 3) wurde gezeigt, da man 
bei zusatzlicher Bericksichtigung dieser Auftriebskrafte zu den aus dem Lorentz- 
Ansatz fiir das Makrofeld abgeleiteten k6rperfesten Kraften auf dasselbe Resultat 
fiir die bewegende Kraft des Feldes auf einen Korper gelangt wie direkt mit Hilfe 
der Maxwellschen ponderomotorischen Volumskrafte. Es sei an dieser Stelle ver- 
merkt, daB diese Ubereinstimmung der Resultate fiir die resultierende pondero- 
motorische Kraft auch bei ganz allgemeinem Zusammenhang zwischen D und & 
bzw. 8 und §, also auch bei anisotroper Materie, besteht. 


B. Die Bereichsabgrenzung polarisierter Materie 


Die Mittelwertsbildung der fiir den Mikrobereich giiltigen Feldgleichungen fiihrt 
auf die Beziehungen des Makrofeldes (s. etwa 1®, Abschn. 109, und 1’). Bei An- 
wesenheit von Materie bedeutet dann die Verwendung der elektromagnetischen Makro- 
groBen zwingend, daf man sich auf den Standpunkt der Kontinuumstheorie gestellt 
hat, was die Autoren in fast allen Fallen gar nicht mehr eigens betonen (s. die bereits 
erwahnten Arbeiten!’), doch soll diese Tatsache, um allen EKinwendungen vorzu- 
beugen, nochmals ausdriicklich an die Spitze dieser Erérterungen gestellt werden. 


13 §. FuBnote 8, 8. 287. 

14 Es wird das Lorentzsche Mafsystem verwendet. 

is §,. FuBnote 8, 8. 287. 

16 ©]. Schaefer: Einfiihrung in die theoretische Physik, Bd. III/1. Berlin: de Gruyter 
u. Co. 1950. 

17 H. Hofmann: Uber die Deutung der Maxwellschen Gleichungen mit Hilfe elektrischer 
und magnetischer Mengen. Acta Physica Austriaca 11, H. 2 (1957). 

18 §, die FuBnoten 1 bis 7, 8S. 286. 
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Die allzu formale Anwendung der Makrogré8en war nun mit ein Grund, dab tiber 
den Kraftangriff des elektromagnetischen Feldes an der Materie bislang so viele 
Unklarheiten und offene Fragen bestanden haben, und damit auch ein Anlaf fir die 
erwihnten Untersuchungen des Verfassers!®. Es sollen daher tiber die Anwendung 
und Bedeutung der MakrogroBen in der Materie auch jetzt im Rahmen der Unter- 
suchungen iiber das Drehmoment nochmals die folgenden Tatsachen angefthrt werden: 


Der makroskopische Mittelwert der mikroskopischen elektrischen Ladungs- 


(el) pes 
ee ge) — o/@) — div € = div (D — B) =o +o (9) 


stellt die resultierende, als Quellendichte des elektrischen Feldes in Erscheinung 
tretende Ladungsdichte dar und wird allgemein ,,freie“‘ elektrische Ladungsdichte 
genannt. Sie ist durch die algebraische Summe der 

ee ee eal durch Leitungsstrom zugefiihrten ,,wahren“ Ladungs- 
dichte pa und der durch Polarisation der Materie frei 


(el) 


gewordenen Polarisationsladungsdichte eo” gegeben. Hs 


P 
soll nun ein beliebig abgegrenzter Bereich von Materie 


betrachtet und dabei zunachst vorausgesetzt werden, 

da sich die Material- und Feldwerte stetig andern. 

; Spriinge dieser Werte werden also vorlaufig ausge- 

Abb. 1 schlossen oder zumindest durch schmale Schichten sehr 
steiler Anderung ersetzt. Die ,,freie“‘ Ladungsdichte 


oe) — a + en im Inneren des Bereiches I fithrt nun allein noch nicht auf 


die Gesamtladung eines richtig abgegrenzten, physikalischen K6rpers. Denkt man 
sich namlich die Materie des Bereiches I entfernt, wobei sich aber an der 


Abb. 2 


Ladungsverteilung sowohl der entfernten als auch der zuriickbleibenden Materie 
nichts andern soll, so zeigt sich, daB auBer den im Inneren der Bereiche I 
und IT auftretenden raumlich verteilten ,,freien‘‘ Ladungen auch an der Grenzflache, 
makroskopisch gesehen flaichenhaft verteilte Polarisationsladungen auftreten. Diese 
waren vor der Entfernung der Materie nicht feststellbar, da sich die beiden flichen- 
haften Ladungsschichten in der Grenzfliche gedeckt und wegen des gleichen Betrages, 
aber der verschiedenen Polaritét makroskopisch aufgehoben hatten. Die beiden 
sich deckenden Flachenladungen 1 und 2 (oder 1’ und 2’) gehéren aber natiirlich zu 
verschiedenen Materiebereichen, und es ergibt sich daher: Die Ladung eines ,,kérper- 
lichen“ Materiebereiches ist nicht nur durch die makroskopisch in Erscheinung 
tretende réiumliche ,,freie‘‘ Ladungsdichte ©, also durch die ,,freie‘ Ladung Q’@) 


Qed =| ay = (9 +o) dV, (10) 


19 §. die FuBnoten 8 bis 11, S. 287. 
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allen gegeben, sondern die Gesamtladung eines k6rperlichen, physikalisch richtig 
abgegrenzten Bereiches ist vom Standpunkt der Kontinuumsauffassung aus gesehen 
durch die in seinem Inneren réumlich verteilte ,,freie‘ Ladung der Dichte 0’) und 
die an der Bereichsgrenze latent vorhandene, jeweils immer von der Materie des 
angrenzenden Bereiches kompensierte, flaichenhaft verteilte Polarisationsladung der 
Dichte of bestimmt: 


QE? = | oD AV + § odd =Q' + 6 of dA. (11) 
4 A a 


Diese Tatsache wurde erstmals bei magnetisch polarisierten Kérpern erkannt. Bei 
einem permanenten Magneten z. B. treten nach Trennung an einer beliebigen Stelle 
sofort diese beiden vorher nicht feststellbaren Flichenladungen entgegengesetzten 
Vorzeichens ungeindert in Erscheinung: 


ee 


Abb. 3 


Die GroBe dieser noch zum korperlichen Bereich gehérigen, in der Grenzflache 
liegenden Flachenladungsdichte ist entsprechend den Definitionen der Polarisations- 
vektoren im elektrischen Feld durch 


of) =Bn (12) 
(1... 4uBerer HKinheitsvektor der Flachennormalen der Bereichsgrenzfliche), im 
magnetischen Feld, das man auf magnetische Ladungen (Mengen) zuriickfiihrt, durch 
Canasta (13) 


gegeben. 

Andern sich an der Grenzfliche zweier Medien die Materialeigenschaften sprung- 
haft, so kompensieren sich die beiden in der Grenzflache liegenden Flachenladungen 
der aneinandergrenzenden Medien 1 und 2 nicht mehr. Ks tritt jetzt die von Null 
verschiedene Differenz der Betrage der Flachenladungsdichten makroskopisch in 
Erscheinung. Diese Differenz, der Sprung der Normalkomponente, ist dann die 
(negative!®) Flachendivergenz des Polarisationsvektors (i = n,): 


— Div8= Pn, + Fit = Pon — Pin= (Pay. (14) 
Von der (negativen) Flaichendivergenz bedeutet also der eine Term die in der Grenzfliche 
liegende Flachenladung des einen, der andere Term die ebenfalls in der Grenzflache 
liegende Flaichenladung des anderen Materiebereiches! 

An Hand der Gln. (10) und (11) 1a8t sich tibrigens leicht nachweisen, daf die 
makroskopisch resultierend in Erscheinung tretende, raumlich verteilte ,,freie“‘ Ladung 
allein nicht die Ladung eines abgeschlossenen, k6rperlichen Materiebereiches sein 
kann. Die resultierende Polarisationsladung eines abgeschlossenen K6érpers mui 
bekanntlich stets Null, die Gesamtladung also gleich der durch Leitungsstrom zu- 
gefiihrten ,,wahren“ Ladung sein. Dies gibt der fiir einen kérperlichen Bereich giiltige 


194 Das negative Vorzeichen hingt mit der Festlegung zusammen, da der Kinheitsvektor 1 
vom jeweiligen Bereich wegweisen soll. ; 
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Ausdruck Gl. (11), nicht aber die ,,freie‘‘ Ladung Gl. (10) richtig wieder: Gl. (11) 
liefert namlich unter Beriicksichtigung von Gl. (9) zunachst 
Qe — | AV + G odd =| PAV + | WAV + ho9dA = 
A Vv A 


Vv Vv 
= Qo) — | divBdV + 6 Bnd. (15) 
Vv A 

Wegen des GauBschen Integralsatzes verschwinden die beiden letzten Summanden 
in Gl. (15) und es wird die Gesamtladung Q© eines kérperlichen Materiebereiches 
ge) = Q©, also wie erforderlich gleich der ,,wahren“ Ladung allein. . 
Zur Unterscheidung der beiden besprochenen Begrenzungsarten von Materie 
wurden die folgenden Bezeichnungen eingefiihrt (s.%°): Die Begrenzungsart eines 
Materiebereiches, bei der nur die makroskopisch in Erscheinung tretende Ladungs- 
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Abb. 4a 
m... ,,mathematische“, 
ph ... ,,physikalische‘‘ Grenzflache in der kontinuumsmaBigen Darstellung 


dichte erfaBt wird, wird ,,mathematische* Bereichsabgrenzung genannt. Die ,,mathe- 
matische Bereichsabgrenzung ist also prinzipiell z. B. bei der Bildung der Mittel- 
werte der GroBen des Mikrobereiches anzuwenden. Begrenzt man hingegen einen 
Bereich physikalisch richtig als korperlichen Bereich, werden also beim Kontinuum 
auBer den im Inneren raumlich verteilten auch die in der Bereichsgrenze liegenden 
noch zum Korper gehorigen flachenhaften Ladungsdichten bericksichtigt und mit- 
erfaBt, so wird diese Begrenzungsart ,,physikalische“‘ Begrenzung des Materiebereiches 
genannt. Die flachenhaften Ladungen in der Grenzflache treten dabei makroskopisch, 
mit Ausnahme an Sprungflaichen, nicht in Erscheinung, da sie von den gleich groBen 
Flachenladungen umgekehrten Vorzeichens des angrenzenden Materiebereiches 
kompensiert werden. Es ist dann beispielsweise die makroskopisch allein in Er- 
scheinung tretende ,,freie“‘ Ladung Gl. (10) die Gesamtladung eines ,,mathematisch“ 
begrenzten Bereiches, die Gesamtladung eines kérperlichen, ,,physikalisch‘‘ begrenzten 
Bereiches ist dagegen durch Gl. (11) gegeben. Sie stimmt gemi® Gl. (15) mit der 
, wahren Ladung dieser Materie, wie es sein mu, iiberein. 

Schon Maxwell beniitzte die beiden verschiedenen Begrenzungsarten, natiirlich 
ohne die hier verwendeten Bezeichnungsweisen derselben. In den §§ 385 und 386 
seines bertihmten Werkes*! berechnet er z. B. das skalare Potential eines kérperlichen 


*» H. Hofmann: Zur Frage der Bereichsabgrenzung elektrisch und magnetisch polarisierter 
Materie. Elektrotechn. u. Maschinenb., 73. Jg., H. 10 (1956). 

*1 J. Cl. Maxwell: Lehrbuch der Elektrizitaét und des Magnetismus. Deutsche Ubersetzung 
von Weinstein. Berlin: Springer-Verlag. 1883. 
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(,,physikalisch“‘) begrenzten Bereiches magnetisierter Materie und zeigt — worauf 
der Verfasser in allen seinen diesbeziiglichen Arbeiten, insbesondere in 2, immer 
wieder verwiesen hat! —, daB die Darstellung polarisierter Materie durch die im 
Bereichsinneren raéumlich und an der Bereichsgrenzfliche flachenhaft verteilten 
elektrischen oder magnetischen Mengen iquivalent ist der Darstellungsweise, die 
den k6érperlichen Materiebereich kontinuierlich aus Doppelschichten, also aus in sich 
polarisierten, ,,physikalisch“ begrenzten Volumselementen aufgebaut sieht (Abb. 4a). 
Maxwell verwendet ferner auch die ,,mathematische“ Bereichsabgrenzung, und zwar 
in § 402a, wo er nachweist, daB in einem Volumen prinzipiell nur dann resultierende 
Magnetmengen zustande kommen kénnen, wenn der Bereich so begrenzt ist, daB 
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Abb. 4b 
m... ,,mathematische“, ; mee 
ph ... ,,physikalische‘‘ Grenzflache in der kontinuumsmaBigen Darstellung 
an der Grenzfliche ,,... Partikel zerschnitten ...““ werden. Es erweist sich, daB 


der Mittelwert der Ladungsdichte eines derartig, naémlich ,,mathematisch“ begrenzten 
Volumselementes auf die ,,freie“‘ Ladungsdichte 0’ = 0, + op fiihrt, im elektrischen 
Feld also auf Gl. (9), bzw. bei Abwesenheit ,,wahrer“‘ Ladungen auf die Polarisations- 
ladungsdichte gp allein, im elektrischen Feld also auf 


ge) == div %. (16) 


Abb. 4b 1a8t erkennen, da8 ein Uberschu8 an Polarisationsladungen in einem ,,mathe- 
matisch*‘ begrenzten Volumselement auch nur dann zustande kommt, wenn sich 
der Polarisationszustand ortlich andert, also Doppelschichten mit verschiedenem 
Dipolmoment aneinandergrenzen. In diesem Fall iiberwiegt im ,,mathematisch™ 
begrenzten Volumselement die Polarisationsladung eines Vorzeichens. Die »physika- 
lischen‘‘ Grenzflachen hiillen hingegen immer ganze Partikel (Doppelschichten) ein 
und ein derartig begrenzter kérperlicher Bereich besitzt als resultierende Ladung, 
wie schon vorhin an Hand von Gl. (15) nachgewiesen wurde, stets nur die allfallig 
vorhandene ,,wahre“‘, aber keine resultierende Polarisationsladung. Die Abbn. 4a 
und 4b zeigen auch, daB im Gegensatz zur ,,physikalischen™ eine ,,mathematische™ 
Grenzfliche niemals zwei aneinandergrenzende Materiebereiche richtig abzugrenzen 
bzw. zu trennen vermag, da sie die an der Grenze liegenden Doppelschichten zer- 


22 §. FuBnote 21, S. 292. 
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schneidet. Daraus wird ersichtlich, daB die nur auf die riumliche ,,freie‘’ Ladungs- 
dichte Bezug nehmenden Ausdriicke fiir ,,mathematische‘ Begrenzung im allgemeinen 
fiir eine praktische Berechnung der Kraft oder des Drehmomentes auf einen K6érper 
keine wesentliche Bedeutung besitzen kénnen. Trotzdem werden sie oftmals, etwa 
in Zusammenhang mit den nur diese réumlichen Dichtewerte enthaltenden Feld- 
gleichungen des elektromagnetischen Feldes bei theoretischen — aber mehr formalen — 
Betrachtungen herangezogen (unter anderem beim Impulssatz, vgl. 78, Kap. IL/1 
und III/1, und *, Abschn. 6). Das Drehmoment hat vollends nur Sinn in bezug auf 
ganze Partikel, und es wird deshalb im folgenden vornehmlich im Zusammenhang 
mit kérperlichen, ,,physikalisch‘‘ begrenzten Bereichen behandelt. Einzig in dem 
Fall der véllig von Vakuum umgebenen Materie lassen sich auch die fiir ,,mathe- 
matische‘ Begrenzung giiltigen Ausdriicke praktisch verwenden. Die im Vakuum 
verlaufende, einen Kérper einhiillende Flache erfaBt jedenfalls alle zu diesem ge- 
horigen Ladungen, auch wenn man sich diesen in einzelne , mathematisch* begrenzte 
Volumselemente. zerlegt denkt. Der etwa vorhandene Uberschu8 der Polarisations- 
ladungen in den ,,mathematisch“ begrenzten Volumselementen des Ko6rperinneren 
und die starken Uberschiisse an Polarisationsladungen in den ,,mathematisch” be- 
grenzten Volumselementen der an das Vakuum grenzenden Zonen, in denen eine steile 
Anderung der Materialeigenschaften und damit des Polarisationszustandes auftritt, 
kompensieren sich gegenseitig. In Summe besitzt ein vom Vakuum umgebenes 
Materievolumen, auch wenn man von der ,,freien‘‘ Ladungsdichte allein ausgeht, 
natiirlich nur ,,wahre“’ Ladung. Die unmittelbare Ubereinstimmung der aus den 
Ausdriicken fiir ,,mathematische und fiir ,,physikalische‘ Bereichsabgrenzung zu 
gewinnenden Resultate fiir den vom Vakuum umgebenen Kérper wurde fiir das 
Potential und die resultierende Kraft in den diesbeziiglichen Arbeiten des Verfassers 
jedesmal nachgewiesen. Sie wird sich in gleicher Weise auch fiir das Drehmoment 
nachweisen lassen. 

AbschlieBend noch eine weitere Bemerkung zur kontinuumsmaBigen Darstellung 
der Materie: In einem makroskopischen Volumselement von Materie befindet sich 
immer noch eine zur Mittelwertsbildung hinreichende Zah] von diskreten Elementar- 
ladungen. Eine derartige Ladungsverteilung besitzt nun im Falle ,,physikalischer” 
Begrenzung nicht nur ein Dipolmoment allein, sondern auch Momente hoherer 
Ordnung (Quadrupolmoment usw.), wobei allerdings die Wirkung dieser héheren 
Momente bei regelloser Anordnung der Partikel viel geringer ist als die des Dipol- 
momentes (s. 2, § 24). Identifiziert man im Sinne der Kontinuumstheorie den Vektor 38 
bzw. 3 mit dem betreffenden Dipolmoment ©? bzw. I™™ je (,,physikalisch* be- 


erenzter) Volumseinheit B—-D-6 am) 
ON es ee 7 

oe ay men) 
een en ie — 18 
vS D Lage (18) 


(Index ph ... ,,physikalische“ Begrenzung), so bedeutet dies, da man die héheren 
Momente vernachlassigt. Die so vereinfachte Darstellung der Materie durch konti- 
nuierlichen Aufbau aus Doppelschichten kann daher bei Problemen, bei denen die 
hdheren Momente bereits eine Rolle spielen, nicht mehr auf richtige Ergebnisse fiihren ; 
auf die Grenzen der Anwendbarkeit der Kontinuumstheorie wurde bereits eingangs 
verwiesen. 


23 §. FuBnote 11, 8. 287. . 

24 HF. Hund: Materie als Feld. Berlin-Gottingen-Heidelberg: Springer-Verlag. 1954. 

2% Becker-Sauter: Theorie der Elektrizitat, Bd. I, 16. Aufl. Stuttgart: Verlag Teubner. 
1957. : 
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C. Der Magnetismus der Materie 


Bekanntlich gibt es zwei Méglichkeiten, den Magnetismus der Materie zu deuten: 
Man kann ihn auf magnetische Mengen oder im Sinne der Elektronentheorie auf 
stationire Elementarstréme zuriickfiihren. Im ersteren Fall soll die Bezeichnung 
», Theorie magnetischer Mengen‘, im letzteren ,,Theorie der Elementarstrome“ ver- 
wendet werden. Die Mengentheorie des Magnetismus wurde bisher praktisch nur in 
der Magnetostatik verwendet, doch kénnen ohne weiteres auch ganz allgemein die 
Maxwellschen Gleichungen eine Auslegung im Sinne dieser Theorie erfahren, so daB 
die Anwendung der Mengentheorie auch im nichtstationdren Fall méglich ist (s. 2), 
Wenn man sich auch heute im allgemeinen fast ausschlieBlich auf die Elektronen- 
theorie bzw. die Elementarstromtheorie stiitzt, so wird bei der praktischen Berechnung 
der Kraft und des Drehmomentes auf magnetisierte Materie doch immer wieder 
gerne auf die Mengentheorie zuriickgegriffen. Die Aquivalenz der beiden Theorien 
in bezug auf das auBerhalb der Materie erzeugte Feld ist allgemein bekannt, ihre 
Aquivalenz hinsichtlich des Kraftangriffes an der Materie wurde in 2’ nachgewiesen 
und sie wird sich auch fiir das Drehmoment nachweisen lassen. Wesentliche Voraus- 
setzung dafiir ist aber, da® stets nur sinngeméiBe Verkniipfungen von GréBen der- 
selben Theorie vorgenommen und nicht Gréfen beider Theorie miteinander vermengt 
werden. Solche sinnwidrige Verkniipfungen finden sich mitunter in der Literatur; 
sie wurden in ?§ (Abschn. V) und 2° (Abschn. 3) besprochen und bilden einen weiteren 
Grund fir die Schwierigkeiten, die im Zusammenhang mit dem Kraftangriff des 
elektromagnetischen Feldes an der Materie bisher bestanden haben. 

Die Mengentheorie beschreibt magnetisierte Materie zunaichst durch die im 
Bereichsinneren raumlich verteilten Magnetmengen, wobei wegen des Fehlens ,,wahrer“‘ 
Magnetmengen die ,,freien‘’ Mengen ausschlieBlich Polarisationsmengen sind. Dieser 
Zusammenhang wird durch 


omen) — div = — div J = os (19) 


dargestellt. An der Bereichsgrenze liegen die noch zum _ ,,physikalisch*‘ begrenzten, 
kérperlichen Bereich gehérigen flachenhaften Magnetmengen der Dichte 


he se (13) 


Die Analyse der Mengentheorie hat gezeigt, daB im Sinne dieser Theorie der Makro- 
vektor § die Bedeutung der Feldstirke des Magnetfeldes erhalt (s.*°, Abschn. 2), 
er ist im Rahmen dieser Theorie fiir den Kraftangriff an den Magnetmengen und 
den Leitungsstrémen heranzuziehen. 

Die Elementarstromtheorie legt dem Magnetismus der Materie Elementarstroéme 
zugrunde, die im Bereichsinneren mit der raumlichen Dichte 


G® = c,rot 3 (20) 
auftreten. An der Bereichsgrenze sind die flichenhaften Stréme 
{2 =e 3 Xn (21) 


(1... iuBerer Einheitsvektor der Flichennormalen der Grenzflache des betrachteten 
Materiebereiches) zu beriicksichtigen. Bei Anwendung der Elektronen- bzw. der 


26 §. FuBnote 17, S. 289. 
27 §, FuBnote 8, S. 287. 
28 §. FuBnote 8, 8. 287. 
29 §, FuBnote 17, S. 289. 
30 §, FuBnote 17, S. 289. 
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Elementarstromtheorie mu weiters beachtet werden, daB hier dem Vektor 8 die 
Bedeutung der magnetischen Makrofeldstiirke zuzumessen ist (s.*, Abschn. 1). 


(tae 


are 


EEE 
y 
Abb. 5 
m... ,mathematische“, 
ph ... ,,physikalische‘‘ Grenz- 


flache in der kontinuums- 


maBigen Darstellung 


Geht man bis auf ein makroskopisches Volums- 
element zuriick, so ist dieses im Rahmen der kontinuums- 
maigen Beschreibung nach der Mengentheorie eine 
Doppelschicht, nach der Elementarstromtheorie ein 
Wirbelring. Wie schon erwahnt, ist der kontinuierliche 
Aufbau der Materie aus Doppelschichten gleichbedeutend 
mit der Darstellung durch die raumlichen und flachen- 
haften Ladungsverteilungen. Dasselbe gilt fiir den kon- 
tinuierlichen Aufbau aus Wirbelringen und die Dar- 
stellung durch die raétumlichen und flachenhaften Stréme, 
was bereits in * gezeigt wurde. Beim Sprung der ma- 
gnetischen Materialeigenschaften an der Grenzflaiche 
zweier Medien kompensieren sich vom Standpunkt der 
, Hlementarstromtheorie‘‘ aus gesehen die beiden in der 
Grenzfliche flieBenden Flachenstrome der aneinander- 


grenzenden Materiebereiche nicht mehr. Ks tritt in emem 
solchen Fall die Differenz der Betrage der Flachenstromdichten als Flachenrotor von ¢) 


(22) 


in der Grenzflache makroskopisch in Erscheinung (n = n,). Aus Verstaéndnisgriinden, 
namlich um die Form der Gl. (22) unmittelbar aus Gl. (21) herzuleiten, wurde in 
Gl. (22) die Bedeutung des Einheitsvektors der Flachennormalen als jeweils nach 
auBen vom Bereich wegweisender Vektor beibehalten. In der Literatur wird der 
Kinheitsvektor an Sprungflachen meist so indiziert, daB er die Bedeutung des in 
den jeweiligen Bereich hinweisenden Vektors besitzt. Dann waren entsprechend 
dieser letzteren Bezeichnungsweise in Gl. (22) die Vektoren § und n in den Vektor- 
produkten des Flachenrotors zu vertauschen. 


Rot C5 3} = Coase X tat Co G1 X My = Conde X MW — Cy nhs X 1 = Cy (YY X NZ 


Ill. Das direkte, korperfeste Drehmoment im elektrostatischen Feld 
A. Ansatz und Umwandlung in ein Volumsintegral 


Das Drehmoment UT), der korperfesten Krafte, die an den im Inneren des ,,physi- 
kalisch“ begrenzten, kérperlichen Bereiches raumlich und an der Bereichsgrenzflache 


flachenhaft verteilten Ladungen angreifen, ist 


Ton = |t x EMV + Hr x of Edd. 
A 
Driickt man gemé8 den Gln. (9) und (12) die Ladungsdichten durch die entsprechenden 
Feldvektoren aus, so folgt 
Cot ST x € (Pn) dA = |r x EdivDaV — |r x Ediv RaV + 


V 4 Vv 


+ Ot x €) (Bn) dA. 


Das Oberflaichenintegral in Gl. (24) wird nun weiter umgeformt, wobei auch der 
GauBsche Integralsatz Anwendung findet: 


$tr x &) (Bu) dAa = $ (ALB) (1 x €) = $ AX {BP; (x x G)} =| dV {Bs (rx x G} = 
4 
=| VB) (v x €),.dV + | (BV) (v x Gav. 


4 V 
31 §. FuBnote 17, S. 289. 
32-8. FuBnote 9, S. 287. 


(23) 


(24) 


(25) 
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Der Index c bedeutet, da der betreffende Ausdruck — hier der Klammerausdruck 
(vr x ©) — beziiglich des vorherstehenden Differentiationsoperators — hier V — konstant 
zu halten ist. Es folgt aus Gl. (25) weiter 


§ (© x 6) (Bm) dA =(t x EdivPaV + [ (BV) r x GaV = 
V Vi 

=r x Ediv PAV + | (BV) x GaV + |v x (BV) EAV = 
V V 


=|r x Ediv $aV — | E x {BV;v} dV +) r x (PB grad) av. (26) 
Vv Vv Vv 
Wegen 
Via] (27) 
(I ... Einheitstensor) ergibt sich schlieBlich 


$( (rx ©) (Pn) dA = |e x Ediv BdV + |r x (Pgrad)EdV —[E x (PI)dV = 
1% V 
= |r x EdivPdV +\ x x (Pgrad) EdV —|E x PaV. (28) 
Vv V V 
Fihrt man Gl. (28) in Gl. (24) ein, so gelangt man zu 
Tey = |r x {EdivD + (P grad) G} dv —(C x PaV = 
V 


x 


=ae x {of E + (P grad) G} dV +1" x EV. (29) 
Vv 


Der in der geschlungenen Klammer des ersten Integrals in Gl. (29) stehende Ausdruck 
ist aber gerade die an einem korperlichen Bereich angreifende Kraftdichte [s. ® 


Gl. (29), (31) und (32)] 


fer = of? E + (P grad) € = div TS, — 7 TE, (30) 
(T ... transponierter Tensor zu T) mit dem Tensor 
Fy =; D—+1=6; DS (DE)] +> (% &) I, (31) 
so daB also aus Gl. (29) folgendes Ergebnis folgt: 
Tom =| x foyav +|B x EdV =[{r x VT) + Bx Gav. (32) 
V Vv Vv 


Geht man schlieBlich auf das Drehmoment je (,,physikalisch‘‘ begrenzter, kérperlicher) 
Volumseinheit, also auf die Momentendichte, tiber, so erhalt man: 
(el) 


axed 
{Vea =t xX itp He Xe (33) 


Die Momentendichte ist somit nicht allein durch. das Moment der Volumskrafte 
gegeben. Es tritt vielmehr noch ein zweiter Summand, das Vektorprodukt aus 
Polarisationsvektor und Feldstirke, auf. Dieser Summand stellt die Momentendichte 
der Kraftepaare dar, die auf die polarisierten Volumselemente einwirken. Diese 
Kriaftepaare liefern einerseits zur resultierenden Kraftdichte keinen Beitrag. Ander- 
seits ist es aber einleuchtend, daB das Moment der resultierenden Volumskrafte allein 
im allgemeinen keineswegs gleich der resultierenden Momentendichte sein kann, da 
eben noch die Momentendichte der Kraftepaare hinzukommen muf. 


Man erkennt diese Zusammenhinge leichter, wenn man das Resultat Gl. (32) 
bzw. (33) nicht, wie eben durchgefiihrt, vom Moment der auf die raumlichen und 


33 §. FuBnote 8, S. 287. 
Ingenieur-Archiv X1/4 ; 21 
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flichenhaften Ladungen ausgeiibten Krafte ausgehend herleitet, sondern wenn man 
die Kraft und das Drehmoment eines einzelnen polarisierten Volumselementes, also einer 
Doppelschicht, betrachtet. Eine Doppelschicht besitzt die beiden Flichenladungen 
dQ,=0,dA=—($hdA (34) 

dQ, = 0,dA = (% f) dA. (35) 
(... Einheitsvektor nach Abb. 6). Herrscht im 


Zentrum O die Feldstirke €, so sind die Feld- 
stiirken in den Ladungsschwerpunkten 1 und 2 


und 


qlee (3 Vv), (36) 
Abb. 6 
6 =€+(FV)E (37) 


(Die weiteren Glieder der Taylor-Entwicklung k6nnen beim differentiellen Fort- 
schreiten um d8/2 als Glieder héherer Kleinheitsordnung exakt vernachlassigt werden.) 
Die resultierende Kraft. auf das durch die Doppelschicht gegebene Volumselement 
ist demnach 


aR), = dQ, G, +.4Q,6,= —(BHdA (e—(FV)ch +@BH{E+(FV)E}. (38) 


Die beiden Krifte + ($ f) €dA heben sich gegenseitig auf und leisten somit keinen 
Beitrag zur Resultierenden. (Sie verursachen jedoch als Kraftepaar ein Drehmoment!) 
Nach Vertauschung der parallelen Vektoren {$ und d3 und unter Beriicksichtigung von 


(d3 t)dA =dsdA =dV (39) 


erhilt man weiterhin 


dX Gm = (BY) Gav. (40) 


Summiert man hiezu noch die Kraft auf die im Volumselement allenfalls vorhandene 
raumlich verteilte ,,wahre“ Ladung, so folgt 


TS Eh = (0 E + (BV) Ch dV = Fe aV (41) 


mit der bereits in Gl. (30) angegebenen Volumskraft der Kontinuumstheorie auf 
korperliche Volumselemente. Wesentlich ist dabei, daB Gl. (40) eben bereits eine 
resultierende Kraft darstellt und daB auBerdem noch ein Kraftepaar an jedem Volums- 
element angreift. Das Kraftepaar ruft ein zusiatzliches Drehmoment hervor, wie 
das bereits in den Gln. (32) und (33) ausgedriickt ist und auch unmittelbar aus dem 
Drehmoment der Doppelschicht abgeleitet werden kann. Dieses ist gegeben durch 


dX Gh) = %1 x U1 + ve x d§.=(r-F) x dQ, 6, +(r +-$)x dQ,6,= 
dg d. 
=(r-4) x {e-(Sv)eli—eyaa +(t +5) x{E+(Frjcl@paa = 
=1 x (BV)E(RNdA +08 x EBHAA. (42) 


Glieder mit d3* verschwinden wieder wegen der hdheren Kleinheitsordnung. Nach 
Vertauschung der parallelen Vektoren 8 und d3 und wegen Gl. (39) folgt weiterhin 


dXon = {t x (BV) E+P x Gav. (43) 


Summiert man noch das Moment der auf die ,,;wahren‘“‘ Ladungen ausgetibten Kriafte 
so ergibt sich i 


dXiohy = {t X (ova E + (BV)E) +B x dV = {x x f24+BxGQav (44) 


Uber das Drehmoment im elektromagnetischen Feld 299 


in Ubereinstimmung mit dem bereits in Gl. (33) erhaltenen Resultat fiir die Momenten- 
dichte 

toy =f X fe =— Dg a On (33) 
Wie schon erwahnt, ist der erste Summand dieser Beziehung das Moment der resul- 
tierenden Volumskrifte, der zweite bedeutet die Momentendichte der an der Materie 
angreifenden Kraftepaare. Es erweist sich somit, da die Darstellung kontinuierlicher 
Materie mit Hilfe der im Bereichsinneren raumlich und an der Kérpergrenze flachenhaft 
verteilten Ladungen, wie erwartet, auch hinsichtlich des Drehmomentes aquivalent 
ist der Darstellung durch kontinuierlich verteilte ,,wahre‘‘ Ladungen und Doppel- 
schichten. 


In Abschn. IT/B wurde erwihnt, da8 man im Spezialfall des véllig vom Vakuum 
umgebenen Koérpers auch mit Hilfe des fiir ,,mathematisch“‘ begrenzte Bereiche 
giiltigen, nur die raumliche ,,freie‘‘ Ladungsdichte beriicksichtigenden Ausdruckes 
auf dasselbe, richtige Resultat gelangen miisse wie mit Hilfe des fiir kérperliche, 
, physikalisch“ begrenzte Bereiche giiltigen Ansatzes. Das Drehmoment der an den 
raumlich verteilten ,,freien‘‘ Ladungen allein angreifenden Krafte ist 


Tye ko. Caserta; (45) 
Vv 4 

(Index m ... ,,mathematische‘‘ Bereichsabgrenzung). Dabei war die fiir ,, mathematisch‘‘ 

begrenzte Bereiche giiltige Volumskraft [s. 34, Gl. (18) bzw. (23)] 
fo oe 6 = Cdiv € = div T3 =v TS (46) 

mit dem symmetrischen Tensor 
1 ~ 

Tae ees Clare). (47) 


Es 14Bt sich daher das Drehmoment eines k6rperlichen Materiebereiches Gl. (23) 
wegen Gl. (45) im allgemeinen Fall auch 


oy = top er oe Cds TS ne tx (Bn) EdA (48) 
A 

schreiben. Man gelangt ja auch bekanntlich immer zu dem jeweiligen Ausdruck fir 
, Physikalische‘ Begrenzung, indem man zusatzlich zu dem Ausdruck fiir ,,mathe- 
matische“‘ Abgrenzung, der sich ausschlieBlich auf die raumliche ,,freie“ Ladungsdichte 
bezieht, noch den Ausdruck fiir die an der Bereichsgrenze latent vorhandene Flachen- 
ladung hinzufiigt. Verlauft nun die Hiillflache zur Giénze im Vakuum, so verschwindet 
das Oberflachenintegral in Gl. (48) wegen $8, = 0 und es wird in diesem einzigen 
Fall Dee ne: da eben bei einem vom Vakuum umgebenen Korper auch die 
,freie* Ladungsdichte allein in Summe alle seine Ladungen erfabt. 


B. Umwandlung in ein Oberflachenintegral 


Es ist im vorangegangenen Abschnitt III/A gelungen, das direkte, kérperfeste 
Drehmoment der auf die riumlichen und flichenhaften Ladungen eines korperlichen 
Materiebereiches ausgetibten Krafte Gl. (23) in das Volumsintegral Gl. (29) bzw. (32) 
umzuwandeln, also 


(23), (29), (32) TEP, =| rx o'-EdV + Or x op EdA = 
V A 


= | (vr x {o? E + (Pgrad) ©} + PB x E)dV =| {ce x fy + Bx EsavV. (49) 


V V 


34 S$. FuBnote 8, S. 287. 
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Der Integrand dieses Volumsintegrals lie& sich dabei zwanglos als Momentendichte 
der auf die polarisierten, kérperlichen Volumselemente ausgeiibten Krafte deuten. 
Da dabei die Volumskrafte {?) bereits eine resultierende Kraftdichte darstellen, 
mute zusitzlich zu deren Moment zwangsliufig noch ein weiterer Summand, ném- 
lich $3 x €, auftreten, der die Momentendichte der an den polarisierten Volums- 
elementen angreifenden Kraftepaare bedeutet. Fiir eine praktische Anwendung 
wiire es nun auferordentlich wichtig, daB die Umwandlung der GI. (49) in ein Ober- 
flichenintegral gelange, wie dies allgemein bei der Bildung der resultierenden Krafte 
moglich war (s. *, Abschn. VIII): Sowohl die kérperfeste als auch die ponderomotorische 
Kraftdichte konnten als Divergenz eines Tensors dargestellt werden, so da man 
die Resultierende % als Volumsintegral iiber eine Kraftdichte { in einfacher Weise 
mit Hilfe des GauBschen Integralsatzes in ein Oberflachenintegral tiber Spannungen p 
umformen konnte: 

%=|fdV=|VTdV=$TnddA = 6 pd. (50) 

Vv 4 A A 

Der Vorteil einer derartigen Umformung liegt darin, da man von den Verhaltnissen 
des Bereichsinneren unabhingig wird und nur die Feldwerte an der AuBenseite der 
Hillfliche bekannt sein miissen. Obwohl die kérperfesten Volumskrafte Gl. (30) 
selbst aus einem unsymmetrischen Tensor Gl. (31) abgeleitet werden, gelingt trotzdem 
eine derartige Umformung in ein Oberflachenintegral auch beim kérperfesten Dreh- 
moment, was auf die zusatzliche Momentendichte der Kraftepaare zuriickzufiihren 
ist. Man geht hiezu am besten wieder vom urspriinglichen Ansatz Gl. (23) fir das 
korperfeste Drehmoment aus: Fiihrt man in diese Gleichung die Beziehungen (9), 
(12) und (46) ein, so ergibt sich zunachst 


(23), (9), (12), (46) SiS, =| x eh EdV + br x of EdA = 
V A 


= |r x EdivEdV + $r x E(Pn dA = 
Us A 
=|[rxVT@dV + $r x (EP) n}d4=J,+$rx{(E;Pnpdd. (51) 
a A 


Vv 
Das zweite Integral ist bereits ein Oberflachenintegral, so da8B nur mehr der erste 
Ausdruck J, umgewandelt werden muB: 


J,=|rx VT) av = —| VT) x dV = —| WT) x rdV + 
V V V 


als Tel) x td) = —\dV (V1) or (Tey) tay. (52) 
Vv Vv 

Zu beachten ist dabei, daB alle hinter V stehenden Grd8en zu differenzieren sind, 

falls sie von den mit V zusammenhangenden Koordinaten abhingen, es sei denn, 

da durch den Index c angedeutet ist, da sie beziiglich des vorherstehenden 

Differentiationsoperators konstant zu halten sind. Es erweist sich hier wieder die 

formale Schreibweise des GauBschen Satzes als Operatorenregel 


{\dVV= $dAn (53) 
Vv A 
als sehr praktisch, da das erste Integral der letzten Form von Gl. (52), in dem sowohl 


~(e] . , ¢ : < : 
T(} als auch r zu differenzieren sind, in einfacher Weise umgeformt werden kann, 


so da weiterhin folgt: 
J,=— p dA (n Te) x ¢ +) (TRAY) xrdV = Gr xT ndA + | (TEV) x rdV. (54) 
d A 


eee A V, V 
35 §. FuBnote 8, S. 287. 
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Fihrt man J, wieder in Gl. (51) ein und verwendet die fiir den Tensor EN giltige 
Beziehung Gl. (47), so ergibt sich 


Tom = $x x ((GE +658 — FEI) nl ad + | BY) x vay. i» 


FaBt man die beiden Dyaden €;€ und €;% zu €;D zusammen, so gelangt man 
schlieBlich zu 


Tm = hr x {(G2—z@1) al GA + | (Tm V) x ra. (56) 


Da nun der Tensor Te gemaB Gl. (47) symmetrisch ist, verschwindet der letzte 
Summand in Gl. (56). [Der Nachweis hiefiir kann etwa bei °°, III. Teil, Abschn. 1b, 

In. (2) und (3) auf S. 230, anliBlich der Behandlung des symmetrischen Spannungs- 
tensors der Mechanik nachgelesen werden.] Die Dyadensumme des ersten Integrals 
in Gl. (56) entspricht aber dem Tensor Tey nach Gl. (31), so daB also endgiiltig er- 
halten wird: 


Uo) = he x (Tg 1) dd = gr x Pom dA . (57) 


Das Resultat Gl. (57) besagt, daB das kérperfeste Drehmoment allein aus dem Moment 
von Spannungen berechnet werden kann. Es handelt sich dabei gerade um diejenigen 
Spannungen, die tiber die Oberflache integriert auch die richtige resultierende, kérper- 
feste Kraft liefern, da sie durch denselben Tensor T(}, bestimmt sind wie die kérper- 
festen Volumskrafte Gl. (30). ‘ 

Zusammenfassend kann also die tiberaus wichtige Feststellung gemacht werden, 
daB einerseits bei der Darstellung der korperfesten Kraft und des kérperfesten Dreh- 
momentes mit Hilfe eines Systems von Spannungen, die an der Oberfliche eines 
k6rperlichen Materiebereiches angreifen, 


Pom = Tmt (58) 
sowohl die resultierende Kraft als auch das resultierende Drehmoment durch die 
beiden anschaulichen Formeln gefunden werden: 


Soy = $ Pom 44 (59) 
cen = Se nied (60) 


Die Kraft und das Drehmoment erscheinen hier direkt als Resultierende der Spannungen 
bzw. als Summe der Momente derselben Spannungen. Bei der Darstellung mit Hilfe 
der Volumskrafte " 
| (30) fipm = div Tim = Tom 
gelangt man anderseits zu 
Beam = | fom 47 (61) 

und entsprechend Gl. (32) zu 

Ton = J {t x fom + B x €}4P. (62) 


Die resultierende, kérperfeste Kraft ist wohl wieder die Resultierende der Volums- 
krafte, das resultierende, kérperfeste Drehmoment setzt sich aber nicht nur aus den 
Momenten der Volumskrifte zusammen, sondern es treten noch die Momente der 


36 —. Lohr: Vektor- und Dyadenrechnung fur Physiker und hiariede 2. Aufl. Berlin: 
Verlag de Gruyter u. Co. 1950. 
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an den polarisierten Volumselementen angreifenden Kriftepaare, die ihrerseits zur 
resultierenden Kraft Gl. (61) naturgemaB nichts beitragen, hinzu. 


C. Der Spezialfall der Isotropie 


Die Verhaltnisse vereinfachen sich, wenn die Materie isotrop ist und daher Pro- 
portionalitat zwischen D und © bzw. $ und © besteht: 


Dee 2 WY) (63) 
¥=D—-CEC=(ec—1)/E=—7E (64) 
(c ... Dielektrizitiitskonstante, zy ... Suszeptibilitét). Da in diesem Fall Polari- 
sations- und Feldvektor parallel sind, verschwindet der Einflu&8 der Kraftepaare 
beim Drehmoment, der Tensor T{*}, der kérperfesten Krafte, mit denen das Feld 
an einem korperlichen, ,,physikalisch’’ begrenzten Materiebereich angreift, wird 


symmetrisch: 
1 ~ 
(31), (63) Up nero, a eG; c— 5 el = en isotr. * (65) 
Die Symmetrie des Tensors bewirkt, da sich das Drehmoment der Volumskrafte 
allein bereits in das Drehmoment der Spannungen umwandeln 1a8t: 


Tein = ux ea isotr. dV eels x (VV Teh) orn) dV — YX (hae) isotr. 11) dA a 
+ | CTR) isotr. V) x tdV a YX (Tee) isotr. 1) dA ay pr x ee isotr. dA (66) 
A 


[vgl. die Erlauterungen zu Gl. (56)]. Man erkennt dies auch unmittelbar an dem 
Vergleich der Gln. (60) und (62), wenn man dabei beriicksichtigt, daB jetzt die Momenten- 
dichte $8 x © der Kraftepaare wegen {ff gleich Null wird. 


IV. Das ponderomotorische Drehmoment im elektrostatischen Feld 


Die Tatsache, da der klassische Maxwell-Tensor oa aus dem die pondero- 
motorischen Volumskrafte gemaB 


pond == diy pend — Toe (67) 
bzw. die ponderomotorischen Oberflachenkrafte (Maxwellschen Spannungen) gemiB 
pond = T ona ™ (68) 
abgeleitet werden, unsymmetrisch ist, 
é 1 
TOa=E;D—F(@OHL (69) 


bewirkt, daB das Moment der ponderomotorischen Volumskrifte allein nicht in das 
Moment der ponderomotorischen Spannungen umgewandelt werden kann, da 


pao fpona VV =| x V Tyona) 4V = $ x (Tn dA +h(T27) xedV = 
a V 


; pond 
=r X pOadd +f (T247) x rdV 
_ Poona vy ( pond ) xv (70) 
wird und der letzte Summand | (T®,V) x rdV wegen der Unsymmetrie des 
4 
Maxwellschen Spannungstensors nicht verschwindet. Es kann also nur entweder 
das Moment der Maxwellschen Volumskrifte oder das Moment der Maxwellschen 


Spannungen Allgemeingiiltigkeit fiir sich beanspruchen. In der Literatur wird der 
Maxwellsche Spannungstensor stets fiir den Fall isotroper Materie abgeleitet, wo er 
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natirlich wieder symmetrisch ist und wo daher die unmittelbare Umformung des 
Momentes der Volums- in das der Oberflichenkrifte méglich wird. Im allgemeinen 
Fall des beliebigen Zusammenhanges D = € + $$ wird dann vielfach die Giiltigkeit 
des klassischen Maxwell-Tensors angezweifelt. So stellen Sommerfeld und Bopp?” 
das weiter nicht begriindete Postulat auf, da® sich die Kraftdichte im allgemeinen 
Fall ebenfalls aus einem symmetrischen Spannungstensor ableiten lassen miisse und 
der demnach mit dem klassischen Maxwell-Tensor nicht iibereinzustimmen brauche, 
damit auch jetzt die einfache Umformung 


ot <td =Or'x pda (71) 
V A 


mdéglich wird. Die beiden Autoren schlagen unter anderem einen ,,symmetrisierten“ 
Maxwellschen Spannungstensor vor, also eine gewaltsam symmetrisch gemachte 
Abwandlung des klassischen Maxwell-Tensors. Die Erérterungen des vorhergehenden 
Kapitels I haben gezeigt, dafi wegen der Polaritét der Materie eine Symmetrie 
des Spannungstensors und damit die einfache Umformung Gl. (71) nicht erwartet 
werden kann. Es mu eben wegen der an der polarisierten Materie angreifenden 
Kraftepaare zu dem Moment der Volumskrafte noch die Momentendichte dieser 
Kraftepaare hinzutreten. Nicht die Umformung Gl. (71) ist bei polarisierter Materie 
,anschaulich*, wie die beiden vorhin genannten Autoren meinen, sondern vielmehr 
die im Falle des korperfesten Drehmomentes nachgewiesene Gleichheit der Beziehungen 
Gl. (60) und (62). Auch fiir das ponderomotorische Drehmoment wird sich zeigen 
lassen, daB zum Moment der Volumskrafte der entsprechende Ausdruck fiir die 
Kraftepaare hinzuzutreten hat. 

Analog zur Berechnung der ponderomotorischen Kraft entsprechend Gl. (5) wird 
auch das ponderomotorische Drehmoment nach der folgenden, allgemeingiiltigen und 
in der Mechanik schon immer angewendeten Methode berechnet: Das ponderomotorische 
Drehmoment mu8 gleich dem Moment der an einem Korper (also einem kérperlichen, 
»physikalisch*‘ begrenzten Materiebereich!) direkt (korperfest) angreifenden Krafte 
sein, zuztiglich des Drehmomentes, das die durch die umgebende Materie an seiner 
Oberfliche angreifenden Auftriebskrafte hervorrufen: 

Zpona = Lon) + Lauttr.- (72) 
Nun wurde schon in *§ (Kap. VIII/1) erwahnt, da die auf anderem Wege, namlich 
aus dem Energieprinzip gewonnenen Maxwellschen ponderomotorischen Krafte zur 
Voraussetzung haben, daB man die bei einer virtuellen Verschiebung auftretende 
Energieinderung in ein Volums- oder ein Oberflachenintegral tiber den Korper allein 
umzuwandeln vermag. Dies bedeutet eine empfindliche EHinschrankung bei der An- 
wendung der Maxwellschen Ausdriicke, und zwar auf die Fille, in denen sich die 
Materie unabhingig von der umgebenden verschieben 1i8t, also im wesentlichen 
auf einen in homogener (fliissiger oder gasformiger), von ,,wahren“ Ladungen freier 
Materie befindlichen Korper. Die Ubereinstimmung des nach der allgemeinen Be- 
ziehung Gl. (72) berechneten ponderomotorischen Drehmomentes mit dem der Maxwell- 
schen Theorie wird sich also prinzipiell nur in diesem Fall zeigen lassen. In der um- 
gebenden Materie rufen die dort korperfest angreifenden elektrischen Feldkrafte 
einen zur Aufrechterhaltung des Gleichgewichtes nétigen mechanischen Spannungs- 
zustand hervor. Dieser lit sich aus einem (hydrostatischen, aerostatischen) Druck p 
ableiten, so daB die Gleichgewichtsbedingung fiir die umgebende Materie 


ch + pm = Gh — grad 0 (73) 


37 §. FuBnote 1, S. 286. 
38 §. FuBnote 8, S. 287. 
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lautet. Fiir die im ladungsfreien, homogenen Medium herrschenden Volumskrafte (ee) 


ergibt sich wegen 0°? = 0 und « = konst. und unter Beriicksichtigung der im elektro- 
statischen Feld giiltigen Bedingung 

rotG=V x= 0 (74) 
nach Gl. (30) 
fed, = (BV) E=(e = 1) EV) E=(e— IV (EE) —Ex V xO} =(€- NV EG)= 


2 ee 2 2 & & 
(ype ee ee ee eee 


2 2 
7 (75) 
Aus den Gln. (73) und (75) folgt weiterhin 
grad = =e Pite (76) 
und daraus nach Integration der in dem umgebenden Medium herrschende Druck 
€ 
p= +6, (77) 
so daB das Drehmoment der Auftriebskrafte 
€ 
Tau = Gt x (—n){S- +0} a4 (78) 
A 


wird. 

Die vom Orte unabhaingige Integrationskonstante C ware aus den Randbedin- 
gungen zu bestimmen. Sie hat aber auf das resultierende Drehmoment keinen Hinflu8, 
da der durch sie bestimmte, an der Oberflaiche des Korpers angreifende konstante 
Druckanteil — Cn wegen 


roby VY Kit (79) 
kein resultierendes Drehmoment hervorruft: 
dr x(—Cn)dA=ChdAnxr=C\dVV xr=0. (80) 
A A 4 


Das ponderomotorische Drehmoment Gl. (72) wird demnach im elektrischen Feld 
unter Benititzung des Oberflaichenintegrals Gl. (57) fiir das kérperfeste Drehmoment 
und unter Verwendung der Resultate Gl. (78) und Gl. (80) 


e€ e e & 
Tha = LH) + Lame. = Fr x (TER n)dd + $e x (—w {Ay +O} da = 


=$rx{(6D—37 Gl) njad +$r x (—n BP)aa = 
=$r x{(ED-FeI-F@B 6) I)n\aa = §rx ((€:D —>@OI)nlaa. (81) 


Der in der letzten Form von Gl. (81) auftretende Tensor ist nun der klassische Maxwell- 
sche Spannungstensor Gl. (69), das aus der allgemeinen Beziehung Gl. (72) abge- 
leitete bewegende Drehmoment stimmt also mit dem Moment der Maxwellschen 


ponderomotorischen Spannungen p(®, tiberein: 


Tpona = : UX (Thonatt) dA = br x pada. (82) 
A 


Damit ist auch geklirt, daB das Moment der Maxwellschen ponderomotorischen 
Volumskrafte f%), wegen der Unsymmetrie des klassischen Maxwell-Tensors allein 
nicht das richtige ponderomotorische Drehmoment liefern kann. Vielmehr laBt sich 
das ponderomotorische Drehmoment unter Verwendung des Volumsintegrals Gl. (32) 
fiir das kérperfeste Drehmoment und durch Umformung des Momentes der Auftriebs- 
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krafte Gl. (78) mit Hilfe des GauBschen Integralsatzes auf folgende Weise in ein 
Volumsintegral umwandeln: 


(72), (32), (78), (79): ae ee “PS ater. aa 
=| {rx fg +B x Qav + $x x (—m {4 +chda = 
A 
eat x fem + 8 x dV + §ddn BO) xt= 
A 


=|{rx VTQ) +8 x Gav +) avr B® x r= 
V 


Va 


=|(t x {7(2;6-5 1) +B x lav —\r x7 Sav + | Bor x rar. (ss) 
Vv 4 Vv 
Unter Beriicksichtigung der Gl. (79) folgt mit Gl. (67) 


W2a=|(r x {7(DsE-FSI-T POI) +B x E)av= 


Vi 


= |{t x (7 Tyona) +B X EFAV = J (eX fyona +B x ©) AV. (84) 
v 

Die Umwandlung des ponderomotorischen Drehmomentes in ein Volumsintegral 
fiihrt also auf die gleiche Form wie die der Gl. (32) fiir das kérperfeste Drehmoment: 
Ks ist das ponderomotorische Drehmoment wieder durch das Moment der pondero- 
motorischen Volumskrafte zuziiglich des Momentes der bei anisotroper Materie im 
Korperinneren auftretenden Kraftepaare gegeben. Dieses Resultat war auf Grund 
der vorangegangenen Ausftihrungen tiber das korperfeste Drehmoment zu erwarten 
und tiberrascht daher jetzt nicht weiter; es klart aber die Frage, welche der beiden 
Darstellungen, diejenige mit Hilfe der Maxwellschen ponderomotorischen Volums- 
krafte oder die andere mit Hilfe der Maxwellschen Spannungen direkt durch Inte- 
gration auf das richtige Resultat fiihrt: Es ist dies also die Darstellung mit Hilfe 
der Maxwellschen Spannungen Gl. (68), die direkt durch Integrale tiber die Ober- 
fliche des betrachteten Korpers die resultierende ponderomotorische Kraft 


pond x Phona dA 7 Tyrona 1 dA (85) 
und das resultierende, ponderomotorische Drehmoment 
(82) Tyona = Ht X Poona = Ft X (Thora) dA (86) 
A A 


liefert. Die ponderomotorischen -Volumskrifte Gl. (67) fiihren wohl auch durch 
Integration iiber das Volumen direkt auf die bewegende Kraft 
Sond | pone dV = | V Tina). (87) 
Vv Vv 
sie liefern aber allein nicht das richtige ponderomotorische Drehmoment, das viel- 
mehr durch J 
(84) Fora =] (tx fha +B Xx G AV =[{v x VW Tyena) +B x E}AV (88) 
Vv Vv 
gegeben ist. Es lést sich damit schon hier das vorwiegend im Magnetfeld behandelte 
und bisher in allgemeiner Form nicht geklarte Problem, ob und wie weit die aus dem 
klassischen Maxwell-Tensor abgeleiteten Ausdriicke im Falle anisotroper Materie 
ihre Giiltigkeit behalten (vgl.*): Die Anwendung der Maxwellschen Spannungen 


39 §. FuBnote 1, S. 286. 
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ist auch beim anisotropen Korper erlaubt und fiihrt unmittelbar, im Gegensatz vAU 
den Maxwellschen ponderomotorischen Volumskraften, auf das richtige pondero- 
motorische Drehmoment. 

AbschlieBend zu diesem Kapitel soll noch, mehr zur Kontrolle, der Nachweis 
erbracht werden, da& das ponderomotorische Drehmoment fur den vom Vakuum 
umgebenen Kérper direkt gleich dem kérperfesten Drehmoment wird. Diese Uber- 
einstimmung ist erforderlich, da ein Moment von Auftriebskraften mangels solcher 
nicht vorhanden sein kann und daher aus Gl. (72) in diesem Fall 


1 (el) 
Piet 00> pn) 0 (89) 


folet. Die Giiltigkeit der Beziehung Gl. (89) erweist sich sofort, wenn man die beiden 
Drehmomente in der Darstellung der Oberflichenintegrale Gl. (57) und Gl. (82) 
miteinander vergleicht und wegen der im Vakuum verlaufenden Hiillflache in den 
Tensoren die Beziehung 


Do = C, (90) 
berticksichtigt : 
) 1 
Tpona,o = Pt X (Thora, ot) dd = Pr x {(€o; De — F (Do G) 1) n} da = 


= $x x {(G;D)— 5 Gl) n} dA = $e x (THhon) dd = To (91) 
A 


Hier fiihrt nattirlich auch das Moment pitae, der allein an den raumlich verteilten 
, freien’ Ladungen angreifenden Krafte auf dasselbe Resultat, da die im Vakuum 
verlaufende Hiillflache diesmal auch bei ,,mathematischer“‘ Bereichsabgrenzung alle 
Ladungen des Korpers richtig erfaBt. Dieser Nachweis wurde bereits in Kap. III/A 
erbracht. 

Im II. Teil dieser Arbeit werden die analogen Untersuchungen im stationaren 
Magnetfeld durchgefiihrt und schlieBlich wird das allgemeine, nichtstationare elektro- 
magnetische Feld behandelt. 

(EHingegangen am 14. Dezember 1957) 


Stationére Laminarstrémung durch Kreis- und Kreisringrohre 
mit fester bzw. beweglicher Innenwandung und bei statischer 
Druckverteilung* ** 


Von K. Karas, Darmstadt 
Mit 5 Textabbildungen 


Zusammenfassung. Die Geschwindigkeitsverteilung und die sekundliche Durchflu8menge 
bei Laminarstrémungen durch Kreis- und Kreisringrohre mit irgendwie gegebenem Druck- 
gefalle kénnen mittels des angegebenen Membrangleichnisses durch die Auswélbung einer Membran 
dargestellt werden, die gleichformig durch irgendeine Spannung S vorgespannt ist und durch 
einen dem Druckgefalle analog entsprechenden Flaichendruck belastet ist. Man erhalt dann 
das Geschwindigkeitsrelief der Laminarstrémung bzw. die sékundliche DurchfluBmenge der- 
selben, wenn man die Auslenkung bzw. das Volumen des Wolbhiigels mit der Konstanten S/s - 1 
multipliziert, wobei x die dynamische Zihigkeit und J die Rohrlange des laminar durchstrémten 
Rohres bedeuten. Ist die Auslenkung in dimensionsloser Weise dargestellt, so braucht man nur 


die Konstante « durch # zu ersetzen, um statt der Isohypsen des Membranwilbhiigels die 
Isotachen des Geschwindigkeitsreliefs zu erhalten. 


* Dem Gedenken an Herrn Dr. Ing. Camillo Korner, weiland o. 6. Professor fiir Strémungs- 
maschinen an der ehemaligen Deutschen Technischen Hochschule in Prag, gewidmet. 


** Bei der Abfassung der Zeichnungen wurde Verf. von Herrn Dipl.-Ing. K. E. Meier- 
Dérnberg in dankenswerter Weise unterstiitzt. 
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Dieses Gleichnis, das in analoger Weise wohl auch fiir andere Querschnitte Gultigkeit hat, 
bietet neben seiner groBen Anschaulichkeit einen bedeutsamen Vorteil in experimenteller Hinsicht, 
namlich bei Ausmessung des Geschwindigkeitsreliefs einer Laminarstro6mung. Es ersetzt diese 
experimentell schwierige kinematische Ausmessung durch die viel einfacher durchfiihrbare, rein 
statische Ausmessung des Membranwélbhiigels. Umgekehrt kann durch die relativ leicht bestimm- 
bare sekundliche DurchfluBmenge der Laminarstré6mung die im allgemeinen viel zeitraubendere 
Kubatur des Wolbhiigels ersetzt werden. 


Da. Tragheitsglieder nicht beriicksichtigt worden sind, gilt das Gleichnis natiirlich nur fiir 


die sog. ,,schleichenden Bewegungen‘‘, die durch verhaltnismaBig niedrige Reynoldssche Zahlen 
gekennzeichnet sind. 


I. Laminarstrémung und Membrangleichnis 


Auf das in Abb. 1 dargestellte Fliissigkeitsprisma von der Linge 1 wirkt auf seine 
linke Stirnflache der Gesamtdruck p(r, y): rd dr, wahrend die rechte Stirnfliche 
keine oder nur eine vernachlissigbar kleine Druckkraft empfangen soll. Die vier 


ROMIChSE & 


Abb. 1. Ein Flissigkeitsprisma mit dem Einheitsdruck p(r, y) an seiner linken Stirnflache und 
den eingezeichneten Schubspannungen t an seinen vier Mantelflachen. Die Strémung erfolgt 
parallel zur horizontalen Rohrachse a 


Mantelflichen empfangen Schubspannungen t, bzw. t,, wobei der Zeiger r auf die beiden 
Flachen hinweisen soll, deren Normalen die Richtung r besitzen, und analog der 
Zeiger y auf jene Flachen, deren Normalen auf den zugehorigen 7-Richtungen senkrecht 
stehen. Dann gilt nach dem Ansatze von Newton mit x als dynamischer Zahigkeit : 


a) ov 


UT AEE 0 Ti gape 


(Ia; “b) 


Hierin bedeutet v(r,p) die gesuchte Geschwindigkeit parallel /, also parallel zur 
Rohrachse a. 


Die Richtungen von 1, und t, auf den durch r und @ gekennzeichneten Mantel- 
flachen sind wie tiblich in der Richtung —v, ihr Zuwachs auf den durch r + dr bzw. 


308 K. Karas: 


gy + dg gekennzeichneten Mantelflichen, jedoch in Richtung +v angenommen 
worden. 

Setzt man voraus, daB das Rohr horizontal ist (die Verallgemeinerung bei schiefem 
Rohr soll spiter gezeigt werden), so ergibt das Gleichgewicht des Flissigkeitsprismas 
wegen der vorausgesetzten stationiren Strémung: 


—1,rdpl+(r, ne dr) (r + dr) dg l— ty drt + (t» + Ser dy) dri + prdpar =0. 


Nach Ausfiihrung der angezeigten Multiplikationen sowie nach Teilung durch den 
nicht verschwindenden Faktor ldrdg erhalt man: 


(4) ona 
ny to tie 4 Be y Br a0, (2) 


Hierin bedeutet »y das Einheitsgewicht der Flissigkeit. Fiuhrt man (la, b) in (2) 
ein, so erhalt man mit 


J(r, 9) = PSP, (3) 


als ortsverinderlichem Druckgefalle bzw. kinematischer Zahigkeit! x/o die Gleichung: 


" (2 oy 1 oy 


o \ or bt a r a po (4) 


Durch (4) wird v(7, vy) bestimmt. Nennt man v, die etwa vorgegebene oder erst zu 
bestimmende Geschwindigkeit der inneren Rohrwandung mit dem Halbmesser 7,, 
so ergibt sich aus (4) die sekundliche DurchfluBmenge Q zu: 


Ta 


Q=2\ | urdpdr + v.72 a. (5) 


oe 


Li 


Oo 


1 Man vergleiche hierzu etwa das ,,Handbuch der physikalischen und technischen Mechanik“‘, 
herausgegeben von F. Auerbach und W. Hort, Bd. 3, 1. Teil, 1927, Artikel: ,,Elastische 
Schwingungen und Wellen“ von F. Auerbach, S. 283—386, insbesondere S. 355, wo fir die 
Schwingungsauslenkung kreisférmiger Membranen [Gl. (145) ebendort] die Differentialgleichung 


angegeben ist: a2 ae Bw Buy 
ae? ae r)? | ag | 
mit a? = rae wo ¢ die raumliche Dichte und 6 die Dicke der Membran bedeutet. Da pe == 7 as 
ist, so folgt hieraus F f <i “aa 
oa Siete Bri tnas =a 
Ersetzt man die linke Seite dieser Gleichung, naémlich die Tragheitskraft o 6 fe durch — p, 


a ame eee : eae 
da @ Oa die Richtung —w, p jedoch in (6) die Richtung + w besitzt, so erhalt man nach 
Erweiterung mit r sofort Gl. (6). 
1 @ 
Gl. (4) kann wegen J = — 7 ~ gema8 (3) und nach Multiplikation mit 2 auch umgeformt 
werden in 


Op 1 a av 1 dy 
es or te or? % ie = 
und stimmt dann iiberein mit Gl. (3, 33¢) in H. Schlichting, Grenzschicht-Theorie, 1951, 
S. 47, wenn man dort wegen des Fehlens von Raumkraften in Richtung z der Rohrachse 
kK, = 0 setzt und auBerdem die linke Seite dieser Gleichung wegen des hier vorausgesetzten 
laminaren und stationaren Charakters der Strémung streicht. Man vergleiche auch W. Miiller, 
Kinfuhrung in die Theorie der zahen Flissigkeiten, 1932, 8. 29, Gln. (7) mit v, statt ® und 
2 


(9), letzte Gleichung mit k, = 0, wobei in (7) ebendort natiirlich der Term se weefallt. 
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Liegt ein Vollkreisquerschnitt vor, so ist in (5) r; = 0 zu setzen: 
Q=2\ \urdgar. (5a) 
00 
Anderseits genitigt die Auslenkung w(r, y) einer mit der Spannkraft S kp/em homogen 


gespannten Kreis- oder Kreisringmembran, die durch den Flachendruck pr, p) be- 
lastet ist, der Differentialgleichung!: 


ow Fw 1 @w 
8(G+rSe+eS) pr) 0, (6) 
wahrend das Volumen ihres Wélbhiigels mit w, als etwa vorgegebener oder erst zu 
bestimmender Auslenkung des Innenrandes der Kreisringmembran mit dem Halb- 
messer 7; durch 
V=2\ | wrdpdr+w,ren (7) 
r¢, 0 
bestimmt wird, das fiir Vollkreismembranen wegen r; = 0 sich zu 
yey | wr dg dr (7a) 
00 
vereinfacht. Der Vergleich der Gl. (4) und (6) einerseits und der Gln. (5) bzw. (5a) 
mit (7) bzw. (7a) fihrt zu dem folgenden 


Membrangleichnis: 


Liegt eine laminare Flissigkeitsstromung durch ein horizontales Rohr der Linge / 
unter dem Druckgefalle J = p(r, ~)/y 1 vor, so wird das Geschwindigkeitsrelief v(r, q) 
bzw. die sekundliche DurchfluBmenge Q derselben anschaulich durch die Auswélbung 
w(r, p) bzw. das Volumen V des Wolbhiigels einer durch den Flaichendruck p(r, ¢) 
belasteten und mittels S homogen gespannten Membran dargestellt, und zwar erhiilt 
man v(r, y) aus w(7r, gy) bzw. Q aus V nach Multiplikation mit der Konstanten S/x I: 
S 


“nl 


vr, p) = wr, oa: Cea) a |= Leas (8a, b) 


Denn ist in (4) und (6) der Faktor der Klammerausdriicke jeweils gleich 1, so er- 
kennt man, da8 der Faktor von r in (4), nimlich wegen (3) der Faktor: 


OT) Se cM g p 


x PL ok x1 


: é oS : 
aus dem von (6), naémlich = durch Multiplikation mit >; erhalten wird’. 


II. Laminarstrémungen mit statischer Druckverteilung 


Im folgenden sollen die gestellten Stromungsaufgaben mittels des anschaulichen 
Bildes ausgelenkter Membranen gelést werden’. Liegt insbesondere hydrostatische 
Druckverteilung p(r, y) vor, so kann letztere in einen konstanten, im Membranmittel- 


2 Das hier ausgesprochene Membrangleichnis, das in sehr anschaulicher Weise die Stro6mungs- 
geschwindigkeiten v direkt den Verschiebungen w der Membran zuordnet, ist verwandt mit der 
Tatsache, da8 man das Stokessche Reibungsgesetz aus dem Hookeschen Elastizitaétsgesetz durch 
Ersatz der Formanderungen durch die Formanderungsgeschwindigkeiten erhalt. Man vgl. hiezu: 
H. Schlichting: Grenzschicht-Theorie, S. 44. 1951. 

3 Hierbei kann man sich in weitem Mafe der Ergebnisse einer friitheren Arbeit des Verfassers: 
Die Auswélbungen der Kreis- und Kreisringmembranen unter hydrostatischem Druck, erste 
Mitteilung: Homogener Spannungszustand, Ingenieur-Arch. 25, 359—380 (1957), bedienen. Diese 
Arbeit soll kurz durch ,,I“ zitiert werden. 
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punkt O wirkenden Druck p’ = yh mit der hydrostatischen Druckhéhe h und in 
einen antimetrischen Druckanteil p’’ = yr cos y (Abb. 1) aufgespalten werden, denen 
beziiglich die Membranauslenkungen w’(r) und w'’(r,y) entsprechen. w’ bzw. w” 
geniigen den beiden Differentialgleichungen : 


d du’ yh : rd) ow”’ Powe y COs De C : 
CE (4), (8): aed are) Pg staee On carl ecommerce 
(9a, b) 


Den Gl. (9a, b) entsprechen die Losungen 
D 
LE 2D i CE ue We yaa At y2-+O,Inr+G,; vw’ = (— es rt Dir +=*) cos Q. 
(10a, b) 


Die Integrationskonstanten C und D miissen aus den Bedingungen des jeweiligen 
Strémungsproblems berechnet werden. Wegen der Linearitaét der Auslenkungs- 
gleichungen (9a, b) ergibt sich dann die Gesamtauslenkung w zu: 
w(r, y) = w'(r) + w'"(r, ¢). (11) 

Zum Volumen V des Wolbhiigels tragt der antimetrische Teil w’’ nach (10b) nichts 
bei, so daB sich mit (9a) das Volumen V nach (7) fiir Kreisringmembranen bzw. 
nach (7a) fiir Vollmembranen berechnet, wenn man jetzt nachtraglich dort w(r, 7) 
durch w’(r) ersetzt. 

AnschlieBend sollen vier Falle von Laminarstromungen infolge statischer Druck- 
verteilung berechnet werden, von denen sich der erste auf ein Kreisrohr, die drei 
anderen auf Kreisringrohre beziehen. Die Losungen der beiden ersten Falle konnten 
unmittelbar der Arbeit (I) entnommen werden, wahrend die der beiden letzten Falle 
sich wenigstens teilweise auf die Ergebnisse jener Arbeit sttitzen und daher nur geringe 
zusatzliche Rechenarbeit erfordern. 


Ill. Vier Zahlenbeispiele 
a) Stromung durch ein Kreisrohr 
Es ergab sich [I (6), (13)] 


nf 


h 
w = t= P); = are 7) COS (12a, b) 
und somit w nach (11) [I, (14)] 


ae (7,7 — r*) (2h +7 cos ¢). (13) 
Fir konstante Auslenkung w= ww, ergeben sich (Abb. 2) zirkulare Kurven 3. Ord- 
nung [I (18)] als Gleichungen fiir die Héhenlinien 


(t+ ¥)(y +24) — ry + — Bhd =0, (14) 


von denen nur die Ovale innerhalb der Membranfliche in Frage kommen. Der 
Membranmittelpunkt ist auerordentlicher Brennpunkt. Aus (12a) erhalt man 
nach (7a): a yhrg! 

[I (22a)) Vato (15) 


Aus (12a) und (15) kann man nach (8a, b) die Geschwindigkeitsverteilung v’(r) und 
die sekundliche DurchfluBmenge Q sofort durch die der Strémung entsprechenden 
GroBen ausdriicken und erhialt fiir den achsensymmetrischen Anteil 


= 


I 
| 
I 


r ' J 
v'(r) Lummi hrar een ei: 


(16a, b) 


a 
as 
= 
a 
2 
rss 
a 
~ 
i 
=e 
2 
- 


d= ; mR mit Jat. 


ie.) 
x 
~ 


Stationére Laminarstr6mung durch Kreis- und Kreisringrohre 311 


Ist namlich p = y h konstant, so ergibt sich aus (3) sofort J — h/l.. Es soll weiterhin 
die erstere Schreibweise beibehalten werden. Im allgemeinen empfiehlt sich in An- 
lehnung an die anderen Fille die Beibehaltung von Polarkoordinaten r, y und die 
Einfiihrung dimensionsloser GréBen gemas: 


ae Pres Ties ip h w Ss ws a 
= ee 5 —— = a et 
e tp f ; r eM Pe a Vilhes as ans 3° 
pails DLS Son 
ee OE ae og (ay byes devexf) 


Dann erhalt man aus (12a) bzw. (13) 


3 
w’ = v(1— @) 4 oder = oO); 


48 
: ops (18a, b) 
w = (1 =?) (2 + @ cos g) gq oder « = (1 — ¢?)(2»+ ecos¢Q). 
Analog hat man nach (17f) und (16a): 
in y(1— ot) 278, oder r= y(1 — ¢*); 
(19a, b) 


3 
v = (1 — @) (2+ 0 cos ¢) an oder £6 = (1 — @?)(2¥+ ocosg). 


Erteilt man den Gréen w’ bzw. w in (18) oder w’ bzw. v in (19) bestimmte konstante 
Werte w,’ bzw. w, und v,’ bzw. v, oder den GréBen « und f bestimmte konstante 
Werte, so sind durch (18) und (19) die Schichtenlinien der Auslenkungen bzw. Ge- 
schwindigkeitsverteilungen in dimensionslosen Veradnderlichen oder Groen dar- 
gestellt. Letztere (19a, b) werden auch Isotachen genannt. 

Analog erhalt man aus (15) und (8b): 

ae 5 5 
Vasu C= F945 = 0303 mit » <1. (20a, b) 

In Abb. 2 sind die Hohen- 
linien nach (18b) fiir »y = 1 und 
die dort angegebenen Werte « 
sowie auch ein vertikaler Meri- 
dianschnitt durch die Wolbflache 
angegeben. Man vergleiche in I 
die Tabelle 1, S. 365. Man braucht 
bloB die GroBe « durch f zu er- 
setzen, um Abb. 2 sogleich als 
Darstellung des Geschwindigkeits- 
reliefs durch Schichtenlinien auf- 
fassen zu k6nnen. Auch alle 
weiteren Folgerungen tiber den 
Ort des Wélbmaximums m und 
Druckmittelpunktes M behalten 


ere sun dae Paronianseprobict! ~ ‘a, °2) Dis’ Hphenliien “der Walbflacte: dex. Vollkreis- 
ihre sinngemaBe Giiltigkeit. In membran unter hydrostatischem Druck fir die durch 
(16a) erkennt man das bekann- die Werte « gegebenen Héhenlagen der zur Basis der 
te Geschwindigkeitsverteilungs- Wolbflache parallelen Schnittebenen 

gesetz und in (16b) das Hagen- 

Poiseuillesche Gesetz fiir laminare Kreisrohrstromung bei konstantem Druck yh wieder. 
Fiir hydrostatisch verteilten Druck ist jedoch durch (19b) ein neues Geschwindigkeits- 
verteilungsgesetz angegeben, wihrend das Hagen-Poiseuillesche Gesetz auch fiir diesen 


Fall ungeindert bleibt. 
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B. Strémung durch ein Kreisringrohr 
Es ergab sich [I (27), (32)]: 


r TL 
iD aS ( 2In2* +rein +r2in“), 
a a z f 
ey (21a, b) 
rere 
0 b= al Pt (r2 +r?7)r— Paes |cos g, 


Abb. 3. Die Hohenlinien der W6lbflache der Kreisrmgmembran mit fester Innenkreisplatte 
unter hydrostatischem Druck fiir die durch die Werte « gegebenen Héhenlagen der zur Basis 
der WoOlbflache parallelen Schnittebenen 


woraus nach (11) sofort w = w’ + w” folgt. Aus (21a) folgt im Verein mit (7), worin 
aber wegen des festen Innenmantels w, = 0 zu setzen ist [1 (40)]: 


h 
Va ee (22) 
n— 


v; 


Aus (21a) ergibt sich nach Multiplikation mit S/x/ das Geschwindigkeitsver- 
teilungsgesetz und aus (22) ebenso die sekundliche DurchfluBmenge, die zu konstantem 
Druck yh und festen Rohrwandungen gehoren. Bei hydrostatisch verteiltem Druck 
p=y(h-+rcosg) andert sich jedoch das Geschwindigkeitsverteilungsgesetz, das 
mittels (17) sogleich in dimensionslosen GréBen angegeben werden soll. Man erhalt: 


2 1 1 1 
B= tag ene + Inlom) —srlno|+[—e + (1 +72) e—segleose. (28) 
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Dieses Ergebnis folgt aus I (33) sofort mit Beachtung von (8a). Ebenso erhailt man 
aus (22) mit (8b) mit »— 1 und m =—2: 


5 


(24) 


In Abb. 3 sind mittels der Héhenlagen « des Membranwolbhiigels wieder die Héhen- 
linien (23) des Geschwindigkeitsprofils (Isotachen) dargestellt, die fiir die — auch 
in den folgenden Beispielen gewahlten — Gréfen » = 1, ~ = 2 berechnet worden 
sind (man vergleiche in I die Tabellen 2, S. 369, und 3, 8. 370). Daselbst sind auch 
die Meridianschnitte des Wélbhiigels durch die é- und 7-Achse angegeben. Man braucht 


wieder blo8 die GréBe « durch £ zu ersetzen, um sogleich die Abb. 3 als Geschwindig- 
keitsrelief auffassen zu diirfent. 


14 v : 
Q= 3 ing Wet — 1) In — (u? — 1] FF = 0-0405 2 
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C. Strémung durch ein Kreisringrohr mit beweglichem Innenmantel, 
der eine vorgegebene Geschwindigkeit v, besitzt 


Auch hier gilt die Lésung (10a), sofern auch hier wieder zunichst das Membran- 
gleichnis zugrunde gelegt wird. Der Geschwindigkeit v, entspricht hier eine vorgegebene 
Auslenkung w,’ des Innenmantels. Die Grenzbedingungen lauten also hier: 

Wr =r) =9; Weary = Ws. (25) 


Hiermit findet man fiir die Integrationskonstanten C, und (,: 


h ; h , 
i (r,3 te) Ws y h eS (74? r;*) Ws 
C= Cg ae Ts Inr, 
1 r ie 4 r 
In In —% 
qT, qT; 
und erhalt schlieBlich aus (10a): 2 
In -& 
h ) ; r 
Wier aa te lad egal = 74) a0 (26) 
n 


w’’ kann unmittelbar von Gl. (21b) tbernommen werden. Benutzt man wieder die 
in (17) eingefiihrten dimensionslosen GréBen, so erhalt man bei Beachtung von (11), 
wenn man noch die zu «, f analogen dimensionslosen GroBen 

CEO'S) Os Pe ae Iie 


= 3; ee ois (27a, b) 
a 


hinzunimmt: a 
ee : eee ae ts °+(14 a — | cos » (28) 
a= 2\r(1 o) +(1 a)” me *s In uw | Ora Te Q 12 0 Y 
Wahlt man wieder » = 1 und wu = 2, so erhalt man fiir die Ausdriicke in den 
eckigen Klammern von (28): 


5 1 
2(1—@) +-apyme=—Ze) baw. —e+Ge-qe=— Me). 
Mit diesen Abkiirzungen ergibt nun die Auflosung von (28) nach cos 
a + Zo) +7 Ine 
COs Y = Wie) (29) 


Fiir «, wurde — 0°3 gewihlt, also eine Auslenkung w,’ (bzw. Geschwindigkeit Vs), 
die der durch den hydrostatischen Druck zu erwartenden Auslenkung (Geschwindig- 
keit) entgegengesetzt gerichtet ist. In der folgenden Tabelle 1 sind die beiden letzten 

4 Man vergleiche hierzu R. v. Mises: Elemente der technischen Hydromechanik, 1. Teil, S. 45, 


Aufgabe 13, 1914. Lésung ebendort S. 208/09. Die Aufgabe wird dort allerdings nicht in voller 
Allgemeinheit durchgefiihrt, daher ist (21a) der vorliegenden Arbeit nicht angegeben. 
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Summanden des Zihlers von (29) und ihre Summe sowie auch der Nenner — N (oe) 
fiir die in der ersten und letzten Spalte angegebenen Werte von o berechnet worden, 
wobei die Funktionen Z(g) und N(e) der Arbeit I entnommen werden konnten’. 
Nach (29) wurden sodann fiir die oben angegebenen «-Werte die Werte cos gy be- 
rechnet, die in Abb. 4 der Konstruktion der Schichtenlinien dienten. In den letzten 
drei Spalten wurden sodann fiir cosy = +1 (y-Achse, Schnitt 4—B) und fir 


Abb. 4. Die Hohenlinien der Wolbflache der Kreisringmembran unter hydrostatischem Druck, 
deren Innenkreisplatte eine vorgegebene Auslenkung besitzt, fiir die durch die Werte « gegebenen 
Héhenlagen der zur Basis der Wolbflache parallelen Schnittebenen 


cos py = 0 (&-Achse, Schnitt C—D) die zu den einzelnen Groen von e gehérigen 
a«-Werte berechnet und so die Meridianschnitte durch die €- und 7-Achse konstruiert, 
wobei zugleich die Schnittpunkte der Hohenlinien mit den Koordinatenachsen ge- 
wonnen wurden, so da die Hohenlinien noch mehr Konstruktionspunkte erhielten. 
Die Schichtenlinie « = 0 schneidet die negative 7-Achse, fiir die cosy = — 1 
gilt, nicht vor dem Punkte 7 = — 1. Denn fir cosy = — 1 folgt aus (29): 
« = Me) — Ze) —7 4, ne = @ — : e+ a taut @?) += -In ¢8. 
Hieraus folgt: hig 1 Tee ae 
4? Co To” 


se 0-1 
al pa Ord, 
5 Man vergleiche I, Gl. (48a, b, c), S. 368, und ebendort Tab. 2, S. 369. 


somit 


ne 
Q 
R BP] or 
See 
© 
I 
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woraus wegen dg <0 auch da < 0 folgt, das heiBt, die Tangente an den Meridian- 
schnitt A—A im Punkte 0 = 1, p= 2 ist schwach gegen die negativen «-Werte 
zugeneigt®. Man vergleiche Schnitt A—B in Abb. 4. 


SchlieBlich erhélt man aus (7) mit (26): 


h 
Va Ge Chore + 


2 PAS ye 
+ 72 [25 72 — 92) — w/|(E ne — 2") +, re (30) 
In —* : 


v; 


Bei Einfiihrung dimensionsloser GroBen nach (17) und (27) erhalt man: 


3 .) Zing —"P— 1D) og 
v= 2 (eet ee —  — + Sh. (80a) 


Setzt man zwecks Kontrolle in (30a) «, = 0, so liegen die Verhaltnisse des Kreis- 
ringquerschnittes mit festen Kreisrandern vor. Es ergibt sich: 


5 Re 2 
V= na i (eo 1) a (use 1) ” a ie oa \ 
__ Vt % (ut —— line Si 
pages Inu 
Hiermit ist fiir » = 1 der entsprechende Ausdruck in I, S. 370 oben, wiedergefunden. 
Setzt man in (30a) v=1, w= 2, «, = — 0°3, so erhalt man 
V = — 001425 Zee , somit @= — 0:01425 —— 


Es tiberwiegt somit die Stromung gegen den hydrostatischen Druck. 


D. Stro6mung durch ein Kreisringrohr mit beweglichem Innenmantel, 
der eine seiner Geschwindigkeit proportionale Bremsung erleidet 


Diesem Stromungsfall entspricht eine Membran mit elastisch gelagerter Innen- 
kreisplatte, die aber gegen Verdrehen gesichert ist (cz = oo). Dieser Fall ist in I, 
5, 8. 375 bis 379, durchgerechnet, und es kann Gl. (69) daselbst mit cz; = co un- 
mittelbar angewendet werden. Man erhalt: 


2h 
N 


: 2 \ 26h 
w= gig (re —) yea S + (— An +r tint + rein) Ae 4 


(31) 


2 y 2 
al caste Bane SA Ge Se, r2)r —*#7*| cos gl; N =2a8 +¢,n. 


c, ist hierbei die Federkonstante der elastischen Lagerung des Innenrohres gegen 
Verschieben. Fiihrt man noch die weitere dimensionslose GréBe 
tsi ce 
Fano 6 32) 


ein, so erhalt man mit den Einfithrungen (17) aus (31): 


ee LS \inie 
2(1 e) + 20vInu|1—g + (i—-) E , 1 1 
l+olhyz +|-¢ + (1+ a] e sag] eos 9. 
(33) 
Auch hier sollen wieder » = 1, w = 2, ferner o = 3 gewahlt werden. Dann erkennt 
man im ersten Bruch die Funktion ¢,(9) in I (7la), 8S. 377, im Faktor von cos Y 
aber die Funktion — N(g), die schon im vorhergehenden Beispiel verwendet und in 


06 == 


§ Zum Beispiel ergibt sich fiir @ = 0°9 « = —0:01185. 
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Tabelle 1 berechnet worden ist. y,(g) und — N(g) sind auch in den beiden ersten 
Spalten von Tabelle 2 fiir die in den Spalten ganz links und ganz rechts angegebenen 
Werte von @ berechnet worden. Hiermit berechnet sich cos y nach (33) zu: 


x — 9;(@) 
— N(e@) 


Die Werte von cos y nach (34) fiir die oben angegebenen «-Werte sind in den sieben 
weiteren Spalten angegeben. 


cos Y = (34) 


Shit 6-0 


Abb. 5. Die Héhenlinien der Wolbflache der Kreisringmembran unter hydrostatisechem Druck 
mit gegen Verschieben elastisch gelagerter Innenkreisplatte fur die durch die Werte « gegebenen 
Hohenlagen der zur Basis der Wolbflache parallelen Schnittebenen 


Die letzten drei Spalten enthalten dann wieder die nach (34) berechneten Werte 
von « fiir cosy = — 1,0, + 1, die der Konstruktion der Meridianlinien in Abb. 5 
dienten. Diese enthalt auch wieder die Schichtenlinien des Membranwolbhiigels fiir 
die dort angegebenen Werte von «; sie vermitteln wieder ein plastisches Bild des 
Wolbhiigels oder, wenn man « durch f ersetzt, des Geschwindigkeitsreliefs der zu- 
geordneten Laminarstromung. 


Mit », =1+olnu=1-+43In2 = 3:07945 findet man die Verschiebung der 
das Innenrohr reprisentierenden Kreisplatte nach I (74a) zu: 


3 3 
fg Ee Dat 0-06085 +; . », = 0:06085 


(35) 
4S wv, % 
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Den zugehérigen Wert «, erhalt man nach (17e) zu: 
CBECTS) 
ya 
Daher umschlieBen in Abb. 5 die Hohenlinien « = 0°1 bis 0°4 den Innenkreis r = 7,, 
wahrend jene mit « = 0°5, 0°6, 0°65 in sich innerhalb der Membranringflache zurtick- 
kehrende Kurven darstellen. Man entnimmt dem Meridianschnitt A—B in Abb. 5 
ferner Xpax = 0°68, somit ist Umax, = 0°68/8 y 7,3/% 1 = 0°085 y 7,3/x 1. 


— 8 - 0:06085 = 0°487 = f,. (35a) 


KS 


Das Volumen des Wolbhiigels V erhalt man unmittelbar aus I (61): 


{(ret =r) 2208 +| (rat — rt) In — (2 — 729] o4}, 


v% 


yphau 


Van 


(36) 
N=228 +¢,In—* 


oder mittels (17) und (32) in méglichst dimensionsloser Form aus: 


Vitae v 7% + i ‘ 5 3 
a S 8p#(1 + olny) {u 1+ [(@ 1) np — (@ — 1)] 6}. (36a) 


Mit y=1, w= 2, o=3 erhalt man aus (36a) 


5 5 
V=015344, somit @Q=0'153 ae 


x 


als sekundliche DurchfluBmenge’. 


7 Fur elliptische und rechteckige Kanalquerschnitte wurden die Laminarstrémungen unter 
statischer Druckverteilung vom Verf. analog untersucht. Im ersteren Falle lieBen sich geschlossene 
Lésungen angeben, wahrend die gut konvergenten Reihen des letzteren Falles eine rasche zahlen- 
maBige Auswertung zulieBen. 

(EHingegangen am 20. Dezember 1957) 


Uber Kérperfarben im Réntgengebiet 
Von J. Braunbeck, Wien 


Zusammenfassung. Aus der Monotonie der Wellenlangenabhangigkeit des Absorptions- 
koeffizienten zwischen Paarbildungsabsorption und K-Absorptionskante wird eine notwendige 
Bedingung fiir das Vorkommen von Korperfarben hergeleitet. Die Existenz dieser Bedingung 
erklart das Fehlen bestimmter Farben in Réntgen-Farbphotographien. 


Bei der Betrachtung von Réntgen-Farbphotographien fallt haufig das Fehlen 
bestimmter Farben auf. Derartige Photographien werden in der Regel fiir Zwecke 
der Roéntgendiagnostik mit durchfallender Strahlung angefertigt, deren Wellenlainge 
kleiner ist als die der K-Absorptionskante des schwersten in radiographisch wirksamer 
Konzentration im Aufnahmeobjekt enthaltenen Elementes. Im Spektralbereich 
zwischen Paarbildungsabsorption und K-Absorptionskante ist der Absorptions- 
koeffizient eine monoton wachsende Funktion der Wellenlinge. Nachstehend wird 
aus dieser Monotonie eine notwendige Bedingung fiir das Vorkommen einer Farbe 
abgeleitet. 


Verallgemeinert man den Farbbegriff auf Wellenlingen auBerhalb des sichtbaren 
Spektralbereiches, so kann man eine Anordnung von drei wellenlangenabhangigen 
Detektoren mit den Empfindlichkeiten H(A); i = 1, 2, 3 als farbempfindliches System 
bezeichnen. Das Zahlentripel [#,(A), H(A), #;(2)] kann als Koordinatenangabe eines 
Vektors H; betrachtet werden. Halt man den Anfangspunkt des Vektors H, fest und 
variiert den Parameter 4, so erhalt man einen Kegel, den sog. Spektralkegel, der die 
Schar aller reellen Farbvektoren einhiillt. 
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Kiner auf das System auffallenden Strahlung mit der spektralen Intensitats- 
verteilung J(4) wird durch die Beziehung 


A, = | B,(A)-I(a) da (1) 
0 
ein Farbvektor A, zugeordnet. 
Bezogen auf ein bestimmtes farbempfindliches System wird eine Strahlung als 


weif bezeichnet, wenn die Koordinaten des ihr zugeordneten Farbvektors untereinander 
gleich sind. 


© 


| B,(a)- (a) da = | B,(A) -T(a)da = j B4(2) -I(A) da. (2) 
0 0 0 


Durchdringt die weiBe Strahlung ein Objekt, so wird sie verfarbt. Bezeichnet man 
die Dicke des Objektes mit 6 und seinen Absorptionskoeffizienten mit 1(A), so gilt 
fir den der verfarbten Strahlung zugeordneten Farbvektor: 
B, = | BA): 1(A) + e-*™ Pda. (3) 
0 
Nach den Regeln der Integralrechnung existieren drei ausgezeichnete Wellenlangen /,, 
fiir die 


bs) 


B, = e-"4) 0 | BA) (a) da (4) 
0 
gilt. 
Wegen (2) koénnen die Integrale in (4) durch einen Faktor C ersetzt werden. 
B, = C+ eo", (5) 


Nun treffe man die Zuordnung der H(A) und damit der 2; zu den Indizes so, daB 
die Ungleichung 


ey rie (6) 
erfillt ist. Wegen der Monotonie der Funktion 7(A) folgt daraus 
B, = B, = B, (7) 


als einschrinkende Ungleichung fiir die Koordinaten des Farbvektors. Das farben- 
empfindliche System wird also nur Farben aufzeichnen, deren Farbvektoren im 
Schnittvolumen des Spektralkegels mit der aus den Ebenen z = 0, x = y und y = z 
gebildeten Pyramide liegen. 

Mitunter verwendet man farbempfindliche Systeme, die nur zwei wellenlangen- 
abhangige Detektoren enthalten. Das Zahlenpaar H(A), 1 = 1, 2, ist dann als Koordi- 
natenangabe eines Vektors in der Ebene deutbar. Bei Variation von / bestreicht der 
Vektor H,;, dessen Anfangspunkt wieder festgehalten wird, einen Ebenensektor, der 
alle im betreffenden farbempfindlichen System méglichen Farbvektoren enthalt. 

Die Gl. (1) gilt unverindert, wihrend die Definition der weiBen Strahlung nunmehr 

| £,(A)- I(A) da = | B,(A)+ (A) da (8) 
0 0 
lautet. Durch einen Absorber verfirbter Strahlung wird auch hier durch die Gl. (3) 
ein Farbvektor B,; zugeordnet. Aus der Monotonie der Funktion 1(A) folgt sodann 


Rie Be (9) 
Die vorkommenden Farbvektoren liegen also in der Uberlappung des Spektralsektors 
mit dem durch die Geraden y = 0 und w = y begrenzten Ebenensektor. 


Technischer Vorteile wegen verwendet man in der Rontgen-Farbphotographie 
oft einen Detektor mit der spektralen Empfindlichkeitsverteilung H(/), der nach- 
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einander drei verschiedenen Strahlungen mit den spektralen Intensitétsvertei- 
lungen J,(A) ausgesetzt wird. Bei der Variation von 4 beschreibt nun der Vektor J; 
den Spektralkegel, auf dem die Kérperfarben der Monochromatoren liegen. Die 
drei Strahlungen und der Detektor haben die WeiBbedingung 


© © 


[ 2(2) -1,(2) da = | B(A)+ 1,(A) da = | B(A)+ I9(a) da (10) 
0 0 0 


zu erfiillen. 
Einem Absorber von der Dicke b mit dem Absorptionskoeffizienten +(A) wird durch 
B,; = | B(A)- IA) ed (11) 


0 
ein Farbvektor B, zugeordnet. Analog zu (4) existieren drei ausgezeichnete Wellen- 


lingen J, welche die Bedingung 


B, =e". [ H(A) + 1A) da (12) 
i) 
erfiillen. Aus (12) folgt wegen der Monotonie des Absorptionskoeffizienten wieder (7). 
Fiir den Fall eines Detektors, der von zwei Strahlungen verschiedener spektraler 
Energieverteilung getroffen wird, erhalt man in analoger Weise die Gl. (9). 
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Zur Theorie der Gleichrichtung 
Von K. Prachar, Math. Inst. d. Univ. Wien 


Mit 4 Textabbildungen 


Zusammentassung. Es wird die Differentialgleichung eines Gleichrichterkreises mit Wider- 
stand und Kondensator aufgestellt und die Gestalt der zugehdrigen Lésungen untersucht. 


I 


Ein Gleichrichter bestehe aus einer Diode (s. Abb. 1) zusammen mit einem 
Kondensator C und dem Verbrauchswiderstand R, den wir als Ohmschen Wider- 
stand annehmen wollen (in Abb. 1 ist die Heizung der Diode weggelassen). Die 
Kennlinie einer ,,idealen Diode“ zeigt Abb. 2. Diese Kennlinie gibt an, welcher 
Strom 7 durch die Diode flieBt, wenn eine gewisse Spannung e an die Diode gelegt 
wird (wobei e positiv gerechnet wird, wenn die Anode héheres Potential hat als die 
Kathode). Wir kénnen diese Kennlinie in der Form schreiben: 


0 eee 


ee Gh OPE) 


(1) 
und I/tgg =A gesetzt wurde (A > 0). 


Wenn auch diese Kennlinie keineswegs genau den tatsichlichen Verhiltnissen 
entspricht — der Knick bei e = 0 ist nicht so ganz ,,abrupt‘‘ und die Kurve ver- 
lauft fir geniigend groBe Werte von e infolge der beschrankten Elektronenemission 
fast parallel zur Abszissenachse (etwa so, wie die strichlierte Kurve aus Abb. 2) 


as 
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pa doch fiir die Rechnung oft die obige idealisierte Form der Kurve zugrunde 
gelegt. 

Legt man nun an die beiden Eingangsklemmen eine Wechselspannung U, sin w t, 
U, > 0, so flieBt durch den Widerstand R der gleichgerichtete Strom J. Zur Be- 
rechnung von J wird in den Lehrbiichern meist angenommen, daB die Spannung w 
am Widerstand konstant ist, das heiBt, 
da die Kapazitaét C so groB ist, daB die 


Y Sin wt 


Abb. 1 Abb. 2 


Wechselspannungen aller auftretenden Frequenzen kurzgeschlossen werden. Sei w 
diese — vorlaufig unbekannte — konstante Spannung. Dann wird e = U, sinwt — u 
und wenn man 7 = f(e) als Funktion von wt betrachtet, ergibt sich das Bild aus 
Abb. 3, in der zum Vergleich die Kurve 7 = 1 U,sinwé strichliert eingezeichnet 
ist. Es gilt ja 

_ A(U,snwt—u) fir U,sinwt >4, | 


a = 
0 Jor U-sinat< 4. j 


(2) 


Hieraus wird der Gleichstromanteil durch Mittelwertbildung ermittelt: 


7% 


Bie gee 
Dae (es 1 ( ; 58; 
I=5-\idot=3— | A(U,sinwt—u)dot= 
0 3 
$x 
— = (U,sin 8 — 19) = +~—* (sin d — 6 c08 9), (3) 


wobei # durch Uysin (F ob 8) =u definiert ist, das heiBt, U,cos? = u. 
Indem man annimmt, da durch 
den Widerstand nur der Gleichstrom- 
anteil flieBt, kann man # aus der 
Gleichung u=T- R, das heiBbt, 
U,.c08 Os= = U, & (sin 8 — 0 cos 9) 
(4) 
berechnen: tg®?—%0=a2/A Rk, 0< 
<0 <2 = zx. Weiter ergibt sich damit 


I=22i1U, (sin 8 — 0 cos 8) 
und “=U cos: 
Der durch den Kondensator flieBende Strom ist gleich 
Ig=1—T. (5) 


Durch Fourier-Entwicklung nach Vielfachen von wt erhalt man schlieBlich noch 
die Amplitude, mit der ein bestimmtes Vielfaches der Grundfrequenz w in Iq auftritt. 
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Wenn auch die obige Art der Berechnung praktisch vollkommen ausreicht und 
auch theoretisch den Kern der Sache trifft, so ist sie doch prinzipiell insofern ungenau, 
als die Annahme einer konstanten Spannung w= uw am Widerstand nicht zulassig 
ist (es kénnte ja durch den Kondensator dann tiberhaupt kein Strom flieBen und im 
ganzen Kreis daher nur Gleichstrom flieBen). Im folgenden soll die exakte Form der 
Lisung der zugehérigen (nichtlinearen) Differentialgleichung untersucht werden 
(natiirlich wird durch die Annahme einer idealisierten Kennlinie ohnehin eine gewisse 
Ungenauigkeit begangen). 

II 


Mit den Bezeichnungen aus I hat man einerseits i = /f(e) und anderseits 
i= (u/R) +C%u. Wegen e = U,sinwt—u ergibt sich fir w die Differential- 
gleichung 

Ci + =f(Upsin ot —u). 
Nach (1) ergibt sich 


U A 
WES GR Cl 


U,sinwt—u) fir u < U,sinoat, (6) 


. U cr: . 
“=— pO fir. «a> U, sitet. (7) 


(Fiir C co geht diese Gleichung in wu = 0 tiber.) Setzen wir noch 
1 A 
Rate rece eg Be (8) 


und fiihren statt wt eine neue Verinderliche ein, die wir wieder mit ¢ bezeichnen 
wollen, so schreiben sich die obigen Gleichungen in der Form 


a@—=—(«+6)u+6U,sint fir u < Ujsint, (9) 
e=—aw fir u> Upsini. (10) 
(x >0, 6 >0.) Die allgemeine Losung von (9) bzw. (10) lautet: 
w= a ap Mla + 6)sint — cost} + K e—~%*% (11) 
bzw. 
ge Kent (12) 


mit beliebigen Konstanten K, K. Wir wollen nun beweisen, daf die Gln. (9), (10) 
eine periodische Lésung besitzen und daB jede andere Lésung gegen diese konvergiert. 
Zunaichst beweisen wir den 


Hilfssatz: Set t, <<t, und seien u(t), u(t) zwee Lésungen von (9), (10) mit 
Uz(t,) < uUp(t,). Dann gilt 
u(t) < u(t) fir t <t<t, (13) 
und 
Ug(to) — Uy(te) < Ue(t,) — Uy(t,)- (14) 


Beweis: Nach bekannten Satzen iiber die Existenz der Lésungen von Differential- 
gleichungen geht durch jeden Punkt der #, w-Ebene genau eine Lésung von (9), (10). 
Da sich die Integralkurven w,(t), w(t) nicht schneiden kénnen, folgt (13) 

Weiter hat man 


or, tg 
| W o(t) dt <\ w(t) dt, 
1 aT 


denn bei festem ¢ ist % eine monoton abnehmende Funktion von u, und zwar sowohl 
fir w < Usinié als auch fiir w > Uysint. Daraus ergibt sich (14). 
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Wir betrachten nun die Lésungen w,(t), u a(t) vo 


n (9), (10), welche der folgenden 

Anfangsbedingung geniigen: u,(0) = — A, u,(0) = Z- dabei soll A eine so groBe 
positive Konstante sein, da die Kune UWih=Ae**, 0 <t <2, ganz in 
u > U,siné verlauft. In 0 <t < 2m ist dann w,(t) mit dieser Kurve identisch und 
es gilt’ also w.(2 2) <A. Die Lésung (11) von (9), welche u(0) = — A erfiillt, lautet: 
o—— mite ae U(% +- 6) ) sin ¢ — cost} + (— A +4 Te iecayr) err (15) 


Indem man A eventuell noch vergréSert, kann man erreichen, daB diese Kurve fiir 
0 <t <2 ganz in wu < Uysint verliuft und iiberdies, daB 


(22) = ge gR (— Lt ert 2n) — A eo tam 
OU, er bs 
ee eee ar) ie 
groBer als — A = u(0) wird. Diese Kurve ist dann fiir 0 <t <2a mit u,(é) 
identisch und es gilt 
a,(0) = — A, 2,(22) > —'A, 
uO) =A, U (2%) <A. 


Da u,(0) < u,(0), so folgt nach (13) u,(2 2) <u,(22). Es mu8 also einen Wert a 
geben, — A <a <A, so daB fiir die Lésung w,(t) von (9), (10), welche der Be- 
dingung u,(0) =a geniigt, gilt: 

U2) = U,(0) =a. 


Damit ist die Existenz einer periodischen Lésung bewiesen. 


Es kann nur eine solche Loésung geben. Ware nimlich u,*(t) eine zweite, hievon 
verschiedene, so miibte 
Ug*(0) — U9(O) = Uo*(2 %) — Uo(2 m2) 


gelten, was einen Widerspruch zu (14) ergabe. 


Ill 
Wir beweisen nun, da jede Losung von (9), (10) in gewissem Sinn gegen diese 
periodische Lésung konvergiert. 


Sei u(¢) irgendeine Losung von (9), (10) und etwa u(0) <<a = u,(0). Nach dem 
Hilfssatz ist dann 
U (0) — U(0) > U2 2) — u(2 7), 


das heiBt, u(2 2) > u(0); ebenso ergibt sich u(4a) > u(22) usw. Bilden wir nun 
die Funktionen /,,(t) = u,(é + 22) — u(t +2n2) = u,(t) — uf +22) fir n = 0, 
1,2,..., so ist jede dieser Funktionen in 0 <¢ <2a monoton abnehmend und 
positiv. Ferner hat man /,(2 2) = fn4i(0) und f,4,(t) <f,(t) fir 0 <t < 2a. Daher 


existieren die Grenzwerte lim /,(0) und Bas mle z) und haben denselben Wert c, 
n> Co 
> 0. Ferner existiert fiir jedes ¢ aus 0 ie Le 2a der Grenzwert eee ) und ist 


wegen f,(0) > fr(t) > fn(2) ebenfalls gleich c. Es existiert also ay eum t) — u(t + 


n—->Co 
+ 2n)} und ist gleich c. Wir beweisen nun, daB c = 0 sein mub. 


Wegen der eben genannten Monotonieeigenschaften ist die Konvergenz der 
Folge f,(t) gleichmaBig in 0 <¢ <2a. Ks gilt ja 


0 <filt) — fasm(t) < (AQ) wa (0 <t <2z). 
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u(t + 22) konvergiert also fiir n — oo gleichmaBig gegen u(t) —¢ (0 <t < 22). 
Da u(t) den Gl. (9), (10) geniigt, folgt 


Ut +-2nm) = — («a + 6)ult+2nnm)+ 6Uysnt (u< U,sint) 
u(t + 2nn2) = — aut + 2nz) (wu > Uysint). 
Daher konvergiert w(t + 22) fiir m — co gleichmaBig gegen 
— (ox + 6) [u(t) —c] + 6Upsint  [ug(t) —¢ < Upsint] 
— (« + 8) [uo(t) — ¢] [wo(t) —¢ > Uysin#]. 
Da die abgeleitete Folge der Folge u(t + 2) gleichmaBig POR PEE ist, dirfen 


Differentiation und ital vertauscht werden: 


lim w(t + 22) = cou ult + 22) = 7 [uw (t) — cl]. 
Es ergibt sich also, daB u(t) — ¢ shenfalle den Differentialgleichungen (9), (10) geniigt. 
Da u(t) — ¢ eine periodische Funktion ist und nach dem oben Bewiesenen u(t) die 
einzige periodische Lésung von (9), (10) ist, so folgt c= 0. 


IV 
Die periodische Lésung, deren Existenz in II bewiesen wurde, verlauft — wie 
jetzt gezeigt werden soll — ganz in u > 0. Dazu geniigt es zu beweisen, da aus 


u(0) = 0 folgt: u(t) > 0 fiir 0 <¢ < 2a; denn nach dem oben Bewiesenen gilt fiir 
die periodische Lésung u,(t) dann sicherlich w(t) > u(t), da sich fiir u(t) < u(t) 
sofort u(0) < u(0) = 0 und u(22) — u(22) > u(0) — u,(0) ergabe [im Widerspruch 
u (14)]. 

Von u(0) = 0 an ist die Lésungskurve w(t) der Gl. (9), (10) solange mit der ent- 
sprechenden Lésungskurve (11) von (9) identisch, als diese im Bereich wu < Uy sint 
verliuft. Die Gleichung der letzteren lautet nach (15) (A = 0) 


OU, 
Ta(aceaen 
Es gilt u’(0) = 0, w’’(0) = 6U,; fiir geniigend kleine ¢ ist u(t) monoton wachsend 
und < U,sin¢. Daher verlauft die Losungskurve w(t), fiir die w(0) = 0 ist, zunachst 
im taneren des Bereiches 0 <t <a, 0 <u < U,sint. Der Punkt, in dem die 
Losungskurve diesen Bereich zum _ erstenmal verlst, méoge dem Parameterwert 
¢ = ¢* entsprechen. Es kann dann nicht w(é*) = 0 sein, denn die Losungen von (9), (10), 
die fiir irgendein ¢* mit 0 <t* <a verschwinden, verlaufen fiir geniigend nahe 
an t* gelegene Werte von ¢, t < ¢*, auBerhalb dieses Bereiches. Es kann aber auch 
nicht u(t*) = U, sin t*, t* < are tg (— 1/x) sein, denn dort wire u’(t*) = — « u(t*) = 
= — « U,sint* < U, cos t*, so dag wieder fiir geniigend nahe an ¢* gelegenes 
t <t* sicherlich [¢, u(t)] nicht im Bereich 0 <t <a, 0 <u < Up sint liegen wiirde. 
Ks ergibt sich also arc tg (— l/x) <t* <a, u(t*) = U,siné*. Fir #* <t <2a 
ist u(t) durch die Gleichung 


u(t) — u(t*) e— * ¢—t*) te 2 ES 9 zt) 


gegeben; denn w(t) gentigt dann (10) und es ist w(t) > Uy siné fiir t* <t < 22a, da 
die Steigung von U)siné fiir ** <<t <a abnimmt, pancetta die Steigung von u(t) 
zunimmt. Daher folgt schlieBlich — wie bohauptet — u(t) SO “fired 4% = 2ix, 
Daraus ergibt sich u(t) > 0 (0 <t <2) und auBerdem 


u(t) = ) sin t — cost + e— geet 


| Up(t) > 0. 
0 
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Vv 
Zur tatsaichlichen Berechnung des Gleichstromanteils 
1 [molt 
pes [ Ug(t 
as 27 R dt 
0 
bemerken wir, da8 die Kurve u,(¢), wenn ¢ von 0 bis 2 a wachst, zunichst im Bereich 
u > U,sint verlaiuft, sodann, etwa fiir t = ¢,, in den Bereich u < U,sint eintritt, 
: : 1 2 é 
diesen etwa bei t= t,, >” <t, <a, wieder verlaBt und fiir t, <¢ < 22 wieder 


in wu > U,sint liegt. Abb. 4 zeigt eine Skizze des Verlaufes von u,(t) und von zwei 
anderen Integralkurven. 


——— 


Abb. 4 


Unter Beniitzung von (12) erhailt man als Bedingung fir die Stellen ¢, und ¢, 


Us Gis teas SHE telt-ontty Sey Uo sin ty; (1S) 
aus (11) ergibt sich 
Le — pe {(% + 4) sin t, — cos tj} + Ke @794 — Uysin, (17) 


und dieselbe Gleichung fiir ¢, an Stelle von ¢,. Eliminiert man K aus beiden 
Gleichungen, so erhalt man 


eet 94 6(1 + «6 + o?) sind, + dcost,} = 
— ola +9)te {(1 + « 6 + «*) sint, + 6 cos bots (18) 


Aus diesen beiden (transzendenten) Gleichungen kénnen ¢, und ¢, bestimmt werden. 
Das Integral itiber u(t) (0 <¢ <2) setzt sich dann aus den beiden folgenden 
Bestandteilen zusammen: 


te 
I Le aan 52 {(~ + 6) sin t — cos t} + Ke @+0t) dp 
ty : 


6U, (EO p (4+ )tn) = 
= Tria {(% + 4) (cost — eosé,) +sint, —sints} +S 4-5 cules 
= a 5 {(sin ¢, — sin t,) + 6(cos t, — cos ty)} 


[unter Beniitzung von (17) und der analogen Gleichung fiir ¢,]; 


20 + ty U U 
; as Rese ele irae ; 
| U,e—*¢— sin t, dt = — —> sin t, - e** (e oath) — e598) = —* (sin tg = sin 4) 


te 
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[unter Beniitzung von (16)]. SchlieBlich wird 
1 

I = 3% \ wit) dt = 
0 


Ook, £( 
2n2a(a+ 6)R 


sin t, — sin#,) + « (cost, — cost,)}. 


Zusatz bei der Korrektur: Ein ausreichend genaues Naherungsverfahren 
zur Bestimmung der periodischen Losung von (9), (10) wurde von J. Kammer- 
loher (Hochfrequenztechnik III, Leipzig 1942, S. 189f.) angegeben. Ersetzen wir 
zum Vergleich mit den dort erhaltenen Formeln in unserem Resultat U)sint durch 


1 
U,coswt (d. h. ¢t durch wt + =) Mit So (tg —t) =9 und = @ (t, + 4h) = — p er- 
ergibt sich aus (16) 


—— 
tg tg y = Ig ee eg) 


Statt dieser Formel ergibt sich loc. cit. S. 192 [Formel (292 b)] naherungsweise 
tg 0 tg p = (xn — 9)/m RC. Driickt man mittels (*) y durch # aus und setzt dies in 
(18) ein, erhalt man eine etwas komplizierte Gleichung zur Berechnung von #, 
welche loc. cit. S. 191 [Formel (291 a)] durch tg? —% = 2/A R [vgl. oben (4)] ersetzt 
wird. Wenn man J durch # und y ausdriickt, findet man: 


J =214U, cos y (sind — 8 cos @. 


(Eingegangen am 31. Dezember 1957) 


Zum Aufsatz: 


Uber eine Verallgemeinerung des Potentiales 
auf Strémungen mit Drehung 


Von Klaus Oswatitsch, Stockholm 


(Osterr. Ing.-Arch. 10, H. 2/3, 1956, S. 239—241) 


Herr W. D. Hayes machte mich in liebenswiirdiger Weise darauf aufmerksam, 
da ich im ersten Satz meines Aufsatzes die voreilige Behauptung aufstellte, in jedem 
Stromungsfeld existiere eine zur Geschwindigkeitsrichtung orthogonale Flachenschar. 

Wie aus Gl. (2) meines Aufsatzes hervorgeht, folgt aus dem Ansatz (1), daB der 
Geschwindigkeitsvektor auf dem Wirbelvektor senkrecht stehen muB: 


fo rot w= 0. (2a) 


Ks wird in der Differentialgeometrie gezeigt, daB Gl. (2a) nicht nur die notwendige, 
sondern auch die hinreichende Bedingung fiir die Existenz einer orthogonalen Flaichen- 
schar auf dem Vektorfeld ww darstellt. 


Meine Ausfithrungen beziehen sich also nur auf Stro6mungen, welche Gl. (2a) 
gentigen. Sie werden auch als komplexlamellar! bezeichnet. Wie groB die Klasse 
dieser Strémungen ist, erscheint noch unsicher. Sicher gehéren ihr die ebene und 
die achsensymmetrische Strémung an. C. Truesdell? nimmt an, daB es eine Reihe 
weiterer solcher Stro6mungen gibt und leitet fiir diese Wirbelsitze ab. 


1 C. Truesdell: The kinematics of vorticity. Indiana Univ. Publ. Science Ser. No. 19, 94 
(1954). 


* C. Truesdell: Verallgemeinerung und Vereinheitlichung der Wirbelsaitze ebener und 
rotationssymmetrischer Flissigkeitsbewegung. Z. angew. Math. Mechan. 31/3, 65—71. 


A. Basch: Franz Jung + 327 


Franz Jung + 


Am 3. Dezember 1957 verschied der emeritierte ordentliche Professor der All- 
gemeinen und Analytischen Mechanik an der Technischen Hochschule in Wien 
Dr. phil. Dr. techn. h. c. Franz Jung unerwartet an den Folgen einer iiberstandenen 
Grippe. Am 14. Mai 1872 als Sohn eines Richters zu Hohenelbe in Bohmen ge- 
boren, studierte Jung nach Gymnasialjahren in Wien, Mies und Prag an der Deutschen 
Universitat in Prag, legte 1896 die Lehramtspriifung fiir Mathematik und Physik 
ab und wurde 1899 zum Doktor der Philosophie promoviert. Der Titel seiner dem 
Gebiete der mathematischen Physik angehérenden Dissertation lautete: ,,Uber das 
Problem der Verteilung der Elektrizitaét auf zwei leitenden Kugeln‘‘. Von 1898 an 
wirkte Jung unter Professor Franz Stark als Assistent fiir Mechanik an der Deutschen 
Technischen Hochschule in Prag und war gleichzeitig als Supplent am Deutschen 
Gymnasium in der Prager Vorstadt Smichow titig. Als Jung 1899 eine Stelle als 
provisorischer Lehrer am Gymnasium in Briix antreten sollte, verzichtete er auf die 
Fortsetzung der sicheren Mittelschullaufbahn, um an der Deutschen Technischen 
Hochschule in Prag bleiben zu kénnen, an der er sich 1904 fiir Mechanik habilitierte. 
Schon im folgenden Jahre wurde seine Lehrbefugnis an die Technische Hochschule 
in Wien tibertragen, wo er die Assistentenstelle fiir Reine Mechanik unter Josef 
Finger tibernahm und gleich mit Abhaltung von Parallelvorlesungen und Vor- 
lesungen tiber Hydromechanik betraut wurde. 1911 erfolgte unter Bestellung zum 
Vorstande des neugeschaffenen Zweiten Institutes fiir Allgemeine Mechanik seine 
Ernennung zum auBerordentlichen, 1919 die zum ordentlichen Professor. Aufgabe 
seines Institutes war die Einfiihrung der Studierenden des Maschinenwesens und der 
Elektrotechnik in die Mechanik, wozu sich bald auch jene der technischen Physiker 
gesellte. Die Entpflichtung, die 1941 auf Grund deutschen Gesetzes erfolgte, unter- 
brach nur ftir kurze Zeit seine Lehrtatigkeit, da sich Jung unter den schwierigen 
Verhaltnissen, denen die Technische Hochschule 1945 gegentiberstand, ihr sofort 
zur Verfiigung stellte, als Honorarprofessor sein Institut wieder tibernahm und bis 
1947 in vollster Tatigkeit blieb. 


Jungs wissenschaftliche Bedeutung erscheint dadurch gekennzeichnet, da er von 
Anbeginn seines Wirkens die Stellung der zu Anfang dieses Jahrhunderts praktisch 
noch wenig verwendeten Vektorrechnung in der theoretischen und technischen 
Mechanik erkannt und selbst durch mehrere Verdffentlichungen gefordert hat, von 
denen hier ,,Ableitung im raéumlichen GroBenfelde“ (Zeitschr. f. Math. u. Phys., 
Bd. 567) und ,,Die Polarableitung in rechtwinklig krummlinigen Koordinaten“ (Akad., 
Wien, 1908) hervorgehoben sein mégen. Auch zur Tensorrechnung, der Weiterbildung 
der Vektorrechnung, hat Jung durch Verdéffentlichungen in den Jahresberichten der 
Deutschen Mathematikervereinigung und in den Sitzungsberichten der Wiener 
Akademie der Wissenschaften beigetragen. In der Literatur wurde vielfach die 
raumliche Ableitung als ,,Jungsche Ableitung“* und der allgemeine unsymmetrische 
Tensor nach Jung als ,,Affinor“, die Tensorrechnung selbst als ,,Affinorrechnung* 
bezeichnet. Andere Arbeiten Jungs betrafen Dynamik des Kurbeltriebes, Massen- 
ausgleich, Hydrodynamik, die elektromagnetischen Feldgleichungen, die Schwere- 
feldgleichungen der Einsteinschen Relativitatstheorie. Die letzte Verdffentlichung 
Jungs war eine im Almanach der Osterreichischen Akademie der Wissenschaften 
erschienene Wiirdigung seines Prager Lehrers Ernst Mach. 

Jung war seit 1914 mit Martha, geb. Schimmek, der Tochter eines Freundes und 
Kollegen seines Vaters, verheiratet. Der Ehe entspro8 ein Sohn Franz, der sich 
dem Studium der Physik zuwandte, aber 1944 ein Opfer des zweiten Weltkrieges 
geworden ist. Von diesem harten Schlag konnte sich Jung, frither eine heitere gesellige 
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Natur, nie mehr erholen. Das schwer getroffene Ehepaar fiihrte seither ein stilles, 
zuriickgezogenes Leben. Jungs Giite und Bescheidenheit haben ihm viel Freunde 
erworben. Seine stete Hilfsbereitschaft wurde von allen, die ihn kannten, immer 


hochgeschatzt. 

Im Studienjahr 1930/31 war Jung Rektor der Technischen Hochschule in Wien, 
die ihm 1947 anlaBlich der Vollendung seines 75. Lebensjahres und der Beendigung 
seiner Lehrtitigkeit das Doktorat der technischen Wissenschaften ehrenhalber verlieh. 

Das Osterreichische Ingenieur-Archiv betrauert in Franz Jung einen seiner Heraus- 
geber, der ihm seit der in schweren Zeiten 1946 erfolgten Griindung angehdrte. 


A. Basch, Wien 
Buchbesprechungen 


Dreidimensionale Uberschallprobleme der Gasdynamik. Von K.-R. Dorfner. (Ergebnisse der ange- 
wandten Mathematik: Heft 3.) Mit 44 Textabb., VII, 150 8S. Berlin-Géttingen-Heidelberg: 
Springer-Verlag. 1957. DM 26.—. 


Auf dem Gebiet der theoretischen Gasdynamik haben die Mathematiker schon viele ausgezeich- 
nete Arbeiten geliefert. Aus dem Mathematischen Institut des Professors Sauer ist nun vorliegen- 
des Werk, das sich mit schwierigen dreidimensionalen Uberschallproblemen beschaftigt, hervorge- 
gangen. Hinweisend auf die Notwendigkeit, Vernachlassigungen durchzufithren, um zu Naherungs- 
lésungen fiir allgemeine raumliche Uberschallprobleme zu gelangen, bespricht der Verfasser die 
Linearisierung der Potentialgleichung und in Anlehnung an Arbeiten von L. Schwartz ein 
Rechenverfahren, das sich fiir das Rechnen mit gewissen unstetigen Funktionen als niitzlich erweist. 
Es handelt sich um das Rechnen mit U-Funktionen, ,,uneigentlichen Funktionen“, das auch in 
die theoretische Gasdynamik Hingang gefunden hat. Ganz allgemein werden im 2. Kapitel das 
,,.Dickenproblem“ und das ,,Auftriebsproblem“ eines Tragfliigels behandelt. Im folgenden Kapitel 
werden vom Verfasser Stémungsfelder mit konischer Symmetrie besprochen und zum Schlu8 ein 
Kapitel iber homogene Stérung als naheliegende Verallgemeinerung hinzugefigt. 

Das Werk bietet sehr viel, erfordert aber der schwierigen Materie wegen ein eingehendes Studium. 

R. Bruniak, Wien 


Statik und Dynamik der Schalen. Von W. Fliigge. Zweite, neubearbeitete Auflage. Mit 
121 Textabb., VIII, 2868S. Berlin-Géttingen-Heidelberg: Springer-Verlag. 1957. Geb. 
DM 28.50. 


Man freut sich, mit der vorliegenden Neuauflage eines der grundlegenden Werke der 
elastizitatstheoretischen Literatur einen lange gehegten Wunsch erfillt zu sehen. Es ist dem 
Verfasser gelungen, durch Beschrankung auf das Grundsatzliche den Umfang des Buches in 
den alten Grenzen zu halten — ein sehr schwieriges Unternehmen im Hinblick auf die starke 
Entwicklung, die das behandelte Gebiet seit dem Erscheinen der ersten Auflage genommen hat. 

Die Darstellung ist auBerordentlich klar und anschaulich. Es sei hier im besonderen die 
schéne Diskussion des eigentiimlichen Verhaltens der hyperbolischen Paraboloidschalen hervor- 
gehoben, bei denen die Membrantheorie auf Differentialgleichungen vom hyperbolischen Typus 
fihrt. Ebenso mu8 auf die sorgfaltige Darstellung der Stabilitatsprobleme hingewiesen werden. 
Uberall tritt deutlich das Bestreben hervor, den mathematischen Formalismus zugunsten einer 
anschaulichen Deutung zurtickzudrangen. 

Gegentiber der ersten, 1934 erschienenen Auflage ist der Abschnitt tiber Membrantheorie, 
der etwa die Halfte des Buches einnimmt, etwas erweitert. An Stelle der Membrantheorie der 
Translationsschalen ist eine Theorie fiir beliecbige Schalenformen getreten. In der Biegetheorie 
wurde die Kegelschale mit linear verinderlicher Wandstirke neu aufgenommen. Auch der 
Abschnitt uber Stabilitaét wurde erweitert. Die Kapitel tiber Faltwerke und tiber Schwingungen 
sind im wesentlichen unveradndert geblieben. 

Das Buch gehért auf den Tisch jedes Statikers. H. Parkus, Wien 


Theorie schallnaher Strémungen. Von K.G.Guderley. Mit 125 Textabb., XV, 3768S. Berlin- 
G6éttingen-Heidelberg: Springer-Verlag. 1957. Geb. DM 42.—. 


Die Theorie der Strémungen mit Mach-Zahlen nahe an eins gehért zu den schwierigsten 
Gebieten der Gasdynamik. Der Grund hiefiir ist in dem physikalisch und mathematisch grund- 
verschiedenen Verhalten der Unterschall- bzw. Uberschallstrémung zu suchen. Wahrend die 
Differentialgleichungen der reinen Unterschallstrémung vom elliptischen und die der reinen 
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Uberschallstrémung vom hyperbolischen Typus sind, erscheinen bei schallnaher Stroémung im 
Stromungsfeld sowohl Uber- als auch Unterschallgebiete und somit beide Gleichungstypen 
nebeneinander. Die Méglichkeit der Linearisierung besteht im allgemeinen nicht mehr. 

Die Entwicklung ist hier noch iiberall im Flu8 und viele Forschungsergebnisse (insbesondere 
auch die des Verfassers des vorliegenden Buches) sind in schwer zuganglichen und zum Teil 
uberhaupt geheim gehaltenen Spezialberichten niedergelegt. Man mu deshalb dem Verfasser, 
einem der fithrenden Fachleute seines Gebietes, fiir die hier gegebene systematische Darstellung 
besonders dankbar sein. Sie gestattet es nicht nur dem nicht unmittelbar an der Entwicklung 
Beteiligten, sich einen Uberblick iiber den Stand der Forschung zu verschaffen, sondern bietet 
vor allem demjenigen, der sich diesem Gebiet aktiv zuwenden will, eine sichere Ausgangsbasis. 

Zum Studium des Buches, das sich auf die reibungsfreie Strémung ohne Warmeleitung be- 
schrankt, ist, soweit es sich um die aerodynamische Seite handelt, die Kenntnis der wichtigsten 
Tatsachen und Methoden der Gasdynamik etwa in dem Umfang, wie sie an den Technischen 
Hochschulen vorgetragen werden, durchaus ausreichend. Uberdies stellt der Verfasser die Grund- 
lagen noch in einem sehr schén geschriebenen einleitenden Kapitel zusammen. Die mathematischen 
Anforderungen an den Leser sind allerdings nicht gering und erfordern neben einer guten Aus- 
bildung intensives Mitdenken. Das Buch wird sowohl den Aerodynamiker wie den Mathematiker 
ansprechen. Der Verfasser ist angewandter Mathematiker. Er vermeidet es aber durchaus, in 
das rein Formale abzugleiten, eine Gefahr, die besonders bei der Behandlung schallnaher Strémun- 
gen nicht gering ist. 

Der Inhalt sei noch kurz angedeutet. Nach einer sehr sorgfaltigen Diskussion der Vernach- 
lassigungen bei der Vereinfachung der Strémungsgleichungen und einem Abschnitt tiber die 
linearisierte Gleichung werden zwei exakte Lésungen der vereinfachten schallnahen Potential- 
gleichung gegeben. Dann folgt die Behandlung der ebenen Strémung mittels der Hodographen- 
methode, der ja gerade fur schallnahe Stroémungen besondere Bedeutung zukommt. Hierbei 
wird im wesentlichen nicht die exakte Hodographengleichung, sondern die der schallnahen 
Potentialgleichung entsprechende, nach Tricomi benannte Gleichung verwendet. Dieser Teil 
des Buches, der etwa zwei Drittel des Gesamtumfanges einnimmt, ist besonders fesselnd und 
bietet eine Fille interessanter Ergebnisse. Er ist allerdings keineswegs leicht zu lesen. Den 
Abschlu8 bildet em kurzes Kapitel uber achsensymmetrische Str6mungen. H. Parkus, Wien 


Spannbetonbau. II. Teil. Von W. Herberg. Mit 240 Textabb., 65 Tafeln im Text und 21 Tafeln 
im Anhang, VIII, 407 8. Leipzig: B. G. Teubner Verlagsgesellschaft. 1957. Geb. DM 26.—. 


Der vorliegende 2. Band bildet eine notwendige und gerne entgegengenommene Ergaénzung 
des 1. Bandes, der 1956 erschienen, den Spannbetonbau theoretisch behandelte und alle im 
Spannbetonkérper zu erfassenden Spannungszusténde mathematisch darstellte, wobei auch 
auf die Baustoffeigenschaften gebiihrend Riicksicht zu nehmen war. 

Wie und auf welchen verschiedenen Wegen die rechnerischen Voraussetzungen im Spann- 
betonbaukérper tatsachlich erfillt werden kénnen, ist in diesem 2. Band in Wort und Bild aus- 
fiihrlich dargestellt. Es ist ein Vorzug dieses Buches, in die grofe Zahl der verschiedenen, 
patentierten Bauverfahren eine iibersichtliche Systematik gebracht zu haben. Dies gilt gleicher- 
mafen fiir die Gestaltung der Spannglieder und deren Verankerung, der zum Anspannen einge- 
setzten Geraite und Bauverfahren, der dem Kraftflu8B angepaBten Fiihrung der Spannstaébe usw. 

Der Spannbetonbau verwendet die Bauformen des Stahlbetonbaues; er ist aber in seinen 
Querschnittsabmessungen und damit im Eigengewicht der Konstruktion bedeutend sparsamer 
und gibt die Méglichkeit, kithnere Bauwerke zu schaffen. Bisher vom Stahlbeton nicht erreichte 
Spannweiten von Balken, Bogen, Platten und Schalen sind ihm méglich. Die Kthnheit der 
Bauten, ausgedriickt in ihrer groBen Schlankheit, findet ihre Grenze nur in der Empfindlichkeit 
fiir Schwingungen. Und gerade dieses rissefreie Ausschwingen liegt im Wesen des Spannbetons. 
Wenn der Verfasser die Absicht hatte, alle diese Fragen in Konstruktion und Bauausfiihrung 
ausfiihrlich zu beleuchten, dann darf wohl behauptet werden, daB dies vdllig gelungen ist. 

Seine tabellarischen Zusammenstellungen ausgeftihrter StraBen- und Eisenbahnbriicken 
zeigen nicht nur die verschiedensten Gestaltungsméglichkeiten auf, sondern geben auch die 
wichtigsten Abmessungen an, mit einem Literaturhinweis, der den Leser zu einer ausfithrlicheren 
Beschreibung des Bauwerkes fiihrt. Rund 40 Briickenbauwerke beschreibt der Verfasser in ihren 
Einzelheiten; iibersichtliche Planzeichnungen und gelegentlich auch Zeichnungen iiber Details 
sowie Lichtbilder unterstiitzen das Verstandnis. 

Die vielfache Verwendung des Spannbetons im Hoch- und Tiefbau behandelt das néchste 
Kapitel. Decken- und Dachkonstruktionen, Balken- und Rahmentragwerke, Wasserbehalter, 
Schornsteine, Leitungsmaste, Druckrohre usw. werden beschrieben und iibersichtliche Werks- 
zeichnungen und Bilder beigebracht. 
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Ein ausfithrlich gehaltenes Schrifttumverzeichnis weist dem Leser weitere Wege. Das 
Stichwérterverzeichnis soll die Beniitzung des Buches als Nachschlagwerk erleichtern. Dem 
Lernenden ein ausgezeichnetes Unterrichtsbuch, dem Suchenden ein Nachschlagwerk, ist dieser 


zweite Teil von Herbergs Spannbeton ein wertvoller Zuwachs der Spannbetonliteratur. ; 
F. Pongratz, Wien 


Technische Schwingungslehre in ihren Grundziigen. Von £. Hiibner. Mit 208 Textabb., XI, 
322 S. Berlin-Géttingen-Heidelberg: Springer-Verlag. 1957. Geb. DM 29.40. 


Der Verfasser sétzt sich zum Ziel, ein Lehrbuch zu schreiben, das die Liicke zwischen den 
schwer zu erarbeitenden Werken von hohem wissenschaftlichem Niveau und den zu wenig in 
die Tiefe dringenden elementaren Biichern iiber Schwingungslehre auszufiillen vermag. Bei 
einem solchen Unternehmen erscheint es geboten, die Grundlagen besonders klar und eindringlich 
herauszuarbeiten. 

Leider wird der Verfasser gerade dieser Forderung nicht gerecht. Es erscheint zunachst 
zumindest didaktisch anfechtbar, wenn an die Spitze der dynamischen Betrachtung der ,,Gleich- 
gewichtssatz der Punktmechanik“ gestellt wird, anstatt von der dynamischen Grundgleichung 
auszugehen. Wenn aber dariiber hinaus die D’Alembertsche Tragheitskraft als Reaktion der 
Massenkraft (,,jene Kraft, die zum Beschleunigen der Masse aufzuwenden ist‘) angesehen wird, 
so ist dies véllig dazu angetan, die klaren Grundlagen der Newtonschen Mechanik zu verwirren. 
Eine derartige, offensichtlich von der elektrischen Analogie angeregte Auffassung zeigt, daf solche 
Analogien, wenn sie zu weit getrieben werden, den Grundlagen der Mechanik Gewalt antun. 

Abgesehen von dieser Schwiche in den Grundlagen, die heute leider ein sehr verbreitetes 
Ubel darstellt, bringt das Werk eine Fiille von Material und eine groBe Anzahl zweckmaBiger 
Methoden, die dem Schwingungstechniker willkommen sein werden. Die vielen, vollstandig 
durchgerechneten Beispiele bringen eine reiche Illustration des Lehrstoffes und erleichtern das 
Studium fiir den Anfinger. Bei der Behandlung des Fadenpendels auf S. 76 wurde allerdings 
ubersehen, daB es sich hier um ein System mit zwei Freiheitsgraden handelt, so daB die angegebene 
Schwingungsdauer unrichtig ist. 

Der gréBte Teil des Werkes ist der Behandlung linearer Systeme gewidmet. Am Ende finden 
sich je ein Abschnitt tiber Messung von Schwingungen und tiber nichtlineare Schwingungen. 
Dabei beschrankt sich der Verfasser auf Systeme mit einem Freiheitsgrad. Bei der Abfassung 
des Werkes, das keinen Anspruch auf Vollstandigkeit erhebt, wurde darauf Wert gelegt, auch 
die neueste schwingungstechnische Literatur einzuarbeiten. G. Heinrich, Wien 


Thermodynamics of Heat-Power Systems. Von f’. W. Hutchinson. Mit zahlreichen Textabb., XIX, 
490 S. Reading, Mass., USA: Addison-Wesley Publishing Company, Inc. 1957. Geb. $ 8.50. 


Die vorliegende Thermodynamik, dessen Verfasser eine zwanzigjahrige Erfahrung als Lehrer 
besitzt, soll den Studenten in folgenden sieben Abschnitten : Source, Storage, Liberation, Transition, 
Transformation, Transportation, Utilization, von den Quellen der Warme bis zu den Vorrichtungen 
(Maschinen) fiihren, in denen thermodynamische Anderungen erfolgen. Der Verfasser ist der 
Meinung, da eine ausfithrliche Darstellung der qualitativen Seite sehr niitzlich sein kann und 
will in seinem Werk einen Kern der Grundbegriffe und Prinzipe bringen. Wenn er somit die quanti- 
tative Seite im groBen anderen Lehrbiichern iiberlaBt, so finden wir doch auch in seinem Buch 
eine groBe Anzahl von Beispielen mit Lésungen, die sehr geeignet sind, das in den Kapiteln aus- 
fithrlich Besprochene exakt anzuwenden. Um auf die Vielseitigkeit des Werkes hinzuweisen, sei 
erwahnt, da8 der Leser z. B. im zweiten Teil (Energy in Storage) unter anderem in die Grund- 
gedanken der Atomenergie und ihrer Gewinnung eingefiihrt wird. 

Die Ansicht des Verfassers, da® eine weitlaufige Darstellung fiir den Studierenden von Nutzen 
sein kann, ist zu begriiBen, wenn sie so ausfallt, wie in diesem Werk, das viele Vorziige besitzt, 
darunter auch die, geschickte Hinweise auf entsprechende Literatur und Probleme zu geben. 
Studenten und Interessenten kann dieses Werk wiirmstens empfohlen werden. 

R. Bruniak, Wien 


Kinfiihrung in die Physik der magnetischen Werkstoffe. Von K.M. Koch. Unter Mitarbeit von 
W. Jellanghaus. (Einfihrung in die Physik der elektrotechnischen Werkstoffe: ITI. Teil.) 
Mit 176 Abb. i. T. und 7 Tab. i. Anhang, VII, 208 S. Wien: F. Deuticke. 1957. Geb. S 186.—. 


Die Werkstoffe der Technik, darunter insbesondere die magnetischen Werkstoffe, zeigen 
heute eine geradezu stiirmische Entwicklung, der in Theorie und Praxis an Hand der zerstreuten 
Literatur zu folgen nicht leicht ist. Die Herausgabe einer mehrbandigen ,,Kinfiihrung in die 
Physik elektrotechnischer Werkstoffe‘‘, von denen der vorliegende Band nur einen Teil darstellt 
bedeutet daher eine begriiBenswerte Aufgabe. ’ 


Buchbesprechungen 331 


Die Anwendung der magnetischen Werkstoffe erfordert ein entsprechendes Ausmaf an Theorie. 
Innerhalb eines gewissen Bereiches geniigt die Kontinuumstheorie, die jedoch nur oberflachliche 
Angaben gestattet; das Eintreten in die Atomphysik und soweit diese nicht ausreicht, in ad hoc- 
Theorien, wird notwendig, sobald das wirkliche Verhalten der ungezahlten magnetischen Werk- 
stoffe erfaBt werden soll. DemgemaB8 beschiftigt sich das erste Kapitel des vorliegenden Bandes 
mit der Entwicklung der Grundbegriffe, das zweite mit den atomistischen Theorien des Magnetis- 
mus ; sodann werden in den weiteren Kapiteln die Elementarprozesse der Magnetisierung, Probleme 
der realen Magnetisierungskurve, ferromagnetische Werkstoffe im magnetischen Wechselfeld, 
die magnetischen Werkstoffe in der Elektrotechnik und die Messung magnetischer GréBen be- 


handelt. Vier Anhange mit Zahlentafeln und Daten, Literaturangaben, Register und Firmen- 
hinweise beschlieBen den Band. 


Die Ausfithrungen sind didaktisch gehalten und sichtlich dem ,,Dienst am Leser‘‘ gewidmet. 
Es ist nun allerdings eine schwierige Aufgabe, dem Leser gleichsam im Fluge die ndétigen 
kontinuumstheoretischen und atomphysikalischen Unterlagen zu vermitteln und es liegt die 
Frage nahe, ob nicht zumindest die ersteren, die kontinuumstheoretischen, vorauszusetzen wiren. 
Auf Seite 15 wird angegeben, da® die Einfiihrung des Vektors % ,,erfahrungsgemaB“ fiir ver- 
schiedene Probleme von Vorteil ist. Angesichts der Geschlossenheit der Maxwellschen Theorie, 
in der dieser Vektor eine fundamentale Rolle spielt, scheinen derartige Auslassungen iiberfliissig. 
Daftr konnte etwa kurz und klar der Sachverhalt, da8 der Flacheninhalt der Hysteresisschleife 
den entsprechenden Verlust angibt, aus dem Poyntingschen Vektor abgeleitet werden. 


Sobald man in die Konstitution und das Verhalten der speziellen magnetischen Materialien 
eintritt, gibt es eine Fille von Einzeltatsachen und zugehérigen Hypothesen oder Theorien, auf 
die hier nicht eingegangen werden kann. Das Kapitel ,,Die magnetischen Werkstoffe der Elektro- 
technik“ ist verhaltnismaBig kurz; doch wird im Vorwort ein eigener, den magnetischen Werk- 
stoffen gewidmeter Band in Aussicht gestellt. 


Der vorlegende Band stellt eine didaktisch betonte, handliche Zusammenstellung der nétigen 
Unterlagen dar, die bei der Verwendung magnetischer Materialien bendtigt werden und kann 
jedem an diesem Gegenstand Interessierten empfohlen werden. L. Kneissler, Wien 


The Analysis of Engineering Structures. Von A.J. S. Pippard und J. Ff. Baker. Dritte Auflage. 
Mit zahlreichen Textabb., XII, 564 S. London: E. Arnold (Publishers) Ltd. 1957. 60s. 


Das Werk behandelt die Grundlagen der Baumechanik fiir Studierende des Bauingenieur- 
wesens an englischen Universititen. In ihm ist die Festigkeitslehre mit der Baustatik vereinigt, 
wobei der Stoff wie folgt gegliedert wurde: Spannungen in geraden und gekriimmten Balken, 
die allgemeinen, auf elastische Gebilde beziiglichen Theoreme, Formanderungen elastischer 
Gebilde (Balken und Fachwerke). Spannungen in statisch unbestimmten Fachwerken, Knickung, 
Durchlauftrager, Rahmen, Momentenverteilungsverfahren (Crosssche Methode), AbriB der Theorie 
des Stahlbetons, Theorie der Bogentraiger, Ringe und Hangebriicken, EinfluBlinien statisch 
bestimmter und statisch unbestimmter Systeme, Gewédlbe (Naherungsberechnung), gekriimmte 
Trager, die senkrecht zu ihrer Ebene belastet sind, Stahlskelettbauten, Theorie der Talsperren, 
Bodenpressung und Erddruck, Experimentelle Methoden zur Spannungsbestimmung und Plasti- 
zitatstheorie. 

Diese Zusammenfassung einerseits und die Auswahl des Stoffes anderseits, die mehr darauf 
abgestellt ist, von allem etwas zu bringen, entspricht den Ausbildungsverhaltnissen in England, 
wo viele nur die Unterstufe des Hochschulstudiums durchmachen, dabei aber eine abgerundete 
fachliche Bildung erhalten sollen. 

Unter diesem Gesichtswinkel betrachtet, erfiillt das Werk, aus welchem die einfache Dar- 
stellungsweise besonders hervorzuheben ist, seinen Zweck. EL. Czitary, Wien 


Atomphysik. Grundlagen und Anwendungsbeispiele fiir Metallkundler und Ingenieure. Von 
F. Regler. (Aluminium Ranshofen Mitteilungen: Sonderheft 4.) Mit 39 Textabb., 7 Tab., 43 S. 
Salzburg-GroBgmain: Technische und Fachverlagshandlung F. Wessiak. 1956. 8 42.—. 


Der Autor vermittelt in vorliegendem Werk mit viel padagogischem Geschick in méglichst 
allgemein verstandlicher Form die theoretischen Grundlagen der Atomphysik und die Technik 
der Kernspaltung einschlieBlich ihrer vielseitigen Anwendung. 

Im einzelnen werden nach einem kurzen Riickblick auf die Entwicklung der Atomistik die 
Elemente der Atomhiille mit ihrer Beziehung zu den Wellenlaingen des elektromagnetischen 
Spektrums und zum periodischen System behandelt. AnschlieBend werden der Spaltungsmechanis- 
mus des Atomkernes, das Prinzip der Teilchenbeschleuniger (Zyklotron und Betatron) und die 
kiinstliche Atomumwandlung beschrieben. In einem besonderen Kapitel wird dann auf die 
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Kernspaltung und die Grundtypen der Atomreaktionen und deren Wirkungsweise eingegangen. 
Der letzte Abschnitt befaGBt sich mit der friedlichen Anwendung der Atomenergie, die insbesondere 
dem Metallkundler und Ingenieur neue Priif- und Forschungsbehelfe bietet, wie etwa die Klarung 
der Diffusion fester Kérper ineinander, die Markierung von Atomen der Zusatzelemente, die 
Verbindung der Metalle bei Kalt- und Warmbearbeitungsprozessen, die Strukturanalyse durch 
Neutronenstrahlenbeugung, die Beeinflussung der technologischen Eigenschaften durch Be- 
schieBen mit Neutronen, die Messung der Wirkungsquerschnitte der Metalle fiir Neutronen ver- 
schiedener Geschwindigkeiten, die Verwendung von Radioisotopen u.a.m. Damit vermittelt 
dieses Werk ein abgerundetes Bild vom heutigen Stand der Atomphysik und -technik und ver- 
leiht dem Leser auBerdem die grundlegenden Kenntnisse fiir die eingehendere Beschaftigung 
mit diesem Wissensgebiet. C. Kammerer, Wien 


Raum und Zahl. Von K. Reidemeister. Mit 31 'Textabb., VII, 1518. Berlin-Géttingen-Heidel- 
berg: Springer-Verlag. 1957. DM 19.80. 


Das vorliegende Werk eines sehr kritischen Denkers besteht aus zehn weitgehend unab- 
hangig voneinander lesbaren Essays, und zwar der Reihe nach betreffend den Ursprung des 
geometrischen Denkens, eine Anzahl von Kegelschnitts- und Winkelteilungsmechanismen (aus 
dem NachlaB von C. Otter, 1598—1660), die Grundlagen der analytischen Geometrie, den 
Unterschied der ,,Gegenden‘“‘ im Raum (Orientierungsfragen), Anschauung und Begriff (erdrtert 
an Hand der transfiniten Ordnungszahlen und der Elemente der kombinatorischen Topologie), 
Geometrie und Logik, eine bemerkenswerte neue Begriindung der Infinitesimalrechnung, eine 
Note iiber Gau8, Geometrie und Zahlentheorie (Konstruktierbarkeitsprobleme) und schlieBlich 
,,Prolegomena einer kritischen Philosophie“. Der Verfasser will mit seinen Darlegungen das 
Wesen mathematischen Denkens und insbesondere die Rolle der Anschaung dabei beleuchten. 
Wer Sinn fiir erkenntnistheoretische Fragen und die Grundlagen der Mathematik hat, wird das 
Buch mit Interesse zur Hand nehmen. W. Wunderlich, Wien 


Zusammenstellung eines MefSwagens fiir Bau- und Raumakustik. Von K.-F. Rittinghaus. 
(Forschungsberichte des Wirtschafts- und Verkehrsministeriums Nordrhein- Westfalen: Nr. 213.) 
Mit 17 Textabb., 87S. Ké6ln und Opladen: Westdeutscher Verlag. 1957. DM 19.80. 


Die fortschreitende Industrialisierung und Motorisierung auf der einen und Leichtbauweise 
auf der anderen Seite machen das Aufsuchen und die Bekaimpfung von st6renden Gerauschen 
bzw. Larmen zu einer dringend notwendigen Aufgabe. Der Autor unternimmt nun den Versuch, 
die Zusammenstellung eines MeBwagens fiir die Bau- und Raumakustik zu beschreiben. 

Kinleitend werden in diesem Heft allgemeine Erérterungen tiber Gré8en und Grundbegriffe 
der Akustik angestellt. Dabei ware zu bemerken, da8 siimtliche Begriffe auf DIN-Blatter bezogen 
sind. Die eigentlichen bau- und raumakustischen Messungen basieren meist auf den Méglichkeiten 
der Elektroakustik. Detailliert werden nun die Messungen von Lautstarke, Nachhallzeit, Luft 
und Trittschalld4mmung (besonders fiir Wohn- und Krankenhauser) bzw. Kérperschall und 
auBerdem die Messung der Diffusion und Deutlichkeit besprochen. Fiir diese Messungen sind 
nicht nur die Schaltbilder angegeben, sondern auch die vom Institut fiir elektrische Nachrichten- 
technik der Rhein.-Westf. Techn. Hochschule Aachen entworfenen und verwendeten Vordrucke 
fiir die Auswertung der MeBresultate. Die fiir die eben angegebenen Messungen notwendigen 
Instrumente sind zu einer Einheit in einem Kraftwagen zusammengebaut, um leicht beweglich 
uberall eingesetzt werden zu kénnen. 

Hine ausfiihrliche Literaturangabe schlieBt diesen Band ab. Die deutliche und klare Dar- 
stellungsweise machen dieses Heft sicher zu einem wertvollen Hilfsmittel fiir jeden mit solehen 


Problemen Beschaftigten. FE, Pléckinger, Wien 
Higentiimer, Herausgeber und Verleger: Springer-Verlag in Wien I, Mélkerbastei 5. — Fiir den Inhalt verantwortlich: Prof. 
Dr. Franz Magyar, Wien IV, Technische Hochschule, Karlsplatz 13. — Druck: Manzsche Buchdruckerei, Wien IX, Lust- 
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Wasserstrahlpumpen zur Forderung von Fliissigkeiten. Von Dr. techn. Ferdinand 
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Assistent an der Technischen Hochschule, Wien. Mit 60 Textabbildungen. VI, 75 Seiten. 
Gr.-8°, 1958. S 90.—, DM 15.—, sfr. 15.40, $ 3.55 

Das Anwendungsgebiet der Strahlapparate in vielen Industriezweigen ist sehr groB. Trotzdem sind in der technischen 

Literatur, wenn man von einigen ausfiihrlichen Aufsitzen absieht, nur spirliche Angaben zu finden, nach denen es mdglich 


ist, Wasserstrahlpumpen zur Férderung von Fliissigkeiten richtig zu dimensionieren. Diese Liicke in der Literatur zu schlieBen, 
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Technische Elektrodynamik. Von Dr.-Ing. Dipl- Ing. Franz Ollendortt, 
Member I. E. E., Member A. I. E. E., Senior Member I. R. E., Research- Professor der _ 
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In Vorbereitung befinden sich folgende Bande: ‘i 
Band II: | 


-Dritter Teil: Kinetische Elektronentheorie. Badsehes (N oise) ei Gent: 
mechanik der Elektronik. 


Vierter ails Kristallelektronik.: 
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Die Welt der Vektoren. Einfithrung in Theorie und Anwendung der Vaktoten™ " 
Tensoren und Operatoren. Von Dr.-Ing. Dipl.-Ing. Franz Ollendorff, Professor der Elektro- 
technik und Vorstand des Elektrotechnischen Laboratoriums, The Hebrew Technical 


i Haifa (Israel). eh 68 Textabbildungen. VIII, 470 Seiten. Gr.-8°. 1950. 
S 228.—, DM 37.50, sfr.39.—, $9.— 


Ganzleinen S 240.—, DM 40.—, sfr.42.—, $9.50. 


Vom Ursprung und Ende der Metaphysik. tine Studie zur Welt- 
anschauungskritik. Von Ernst Topitsch, Professor an der Universitat Wien. IV, 320 Seiten. 


Gr.-8°. 1958. i Ganzleinen S 156.—, DM 26. —, sir. 26.60, $6. 20° : 


. Angesichts der néehet prekaren geistigen Situation der Gdeuawart: anBesiaite der bedenk- 
iss Stagnation des philosophischen Denkens unserer Zeit infolge seiner Aufspaltung in eine 
gleichermaBen sterile wie vorlaute existentialistische ,Katastrophenphilosophie* und neumittel- 
alterliche ,Restaurationsphilosophie‘ kann man gar nicht erfreut genug sein iiber die vorliegende 


Publikation, deren Grundgedanken schon infolge ihrer unbestreitbaren heuristischen Frucht-- 
barkeit ohne Zweifel einmal einen ,Eckstein der philosophischen Erkenntnistheorie« bilden : 
werden... Mit diesem Werk ist ein Grundstein gelegt, auf dem in Zukunft mit Aussicht auf 


groBten Gewinn weitergebaut werden kann. Hs bleibt nur zu wiinschen, daS auch tatsichlich 


in dieser Richtung weitergedacht werde, daB die so notwendige kritische Selbstbestimmung 


der Gegenwartsphilosophie mehr und mehr um sich greife!. . .“ Neue Wege 
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